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Effect of Surface Condition of a Sphere on Its Water-Entry Cavity 


ALBERT May 
Naval Ordnance Laboratory, White Oak, Silver Spring, Maryland 
(Received April 3, 1951) 


An investigation of factors which influence the formation of a cavity when a sphere enters water ver- 
tically is described. It was found that cavities are formed at much lower entry speeds for contaminated 
spheres than for clean ones, but that there is little dependence on the type of the contaminant. In regard 
to the scaling of this effect, the tendency to form a cavity varies little with change of size or specific gravity 
of the sphere. Water-entry cavities are observed to form at velocities much lower than those required for 
incipient cavitation in the water tunnel; consequently, a relation between the two effects appears improb- 
able. Some dependence of cavity formation on viscosity might be inferred from the observation that cavities 
were formed at lower entry speeds by spheres coated with liquids of high viscosity than by spheres without 


this coating. 





INTRODUCTION 


HE dependence of water-entry behavior on the 
surface condition of the missile was first investi- 
gated by Worthington and was reported by him! about 
fifty years ago. He studied the phenomena accompany- 
ing the vertical water entry of spheres and reported a 
strong dependence on the surface condition. He found 
that a clean sphere dropped into water with very low 
velocity would produce a “smooth” splash consisting 
of an upward jet unaccompanied by a cavity. With an 
increase of height from which the sphere was dropped 
it was found that the splash became “‘rough”’; a cavity 
was formed and a “basket-shaped” splash was pro- 
duced. The transition from smooth to rough splash 
was gradual but strongly surface dependent ; prolonged 
and careful polishing resulted in increased heights from 
which the sphere could be dropped while still obtaining 
a smooth splash. 

Initially a test was carried out to verify the conclu- 
sion suggested by Worthington’s work, that the mere 
handling of spheres previously cleaned in a solvent was 
sufficient to alter the nature of the water entry. This 
verification was easily obtained. Figure 1 shows por- 
tions of two high speed motion pictures taken of such 
water entries. Smooth steel spheres one inch in di- 
ameter were dropped 6.2 feet, so that they struck the 

‘A.M. Worthington and R. S. Cole, Trans. Roy. Soc. (London) 


IMA, 175 (1900). A. M. Worthington, A Study of Splashes 
(Longmans, Green, and Company, London, 1908). 


water with a speed of 20 ft/sec. Before they were 
dropped the spheres were washed in alcohol and dried. 
One of them was then dropped into the water from a 
clean tissue to prevent contamination of the surface, 
while the second was purposely contaminated by a 
slightly exaggerated handling such as might naturally 
result when the spheres were dropped by hand into the 
tank. The photographs of Fig. 1 were taken through a 
glass-sided tank at about 1900 frames/sec with an 


Eastman Type III high speed camera with synchronized 
flash lighting.” 
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Fic. 1. High speed motion pictures of water entry of clean and 
handled spheres at 20 ft/sec. Upper pictures: one-inch steel sphere 


cleaned in alcohol. Lower pictures: clean sphere with subsequent 
handling. 


2R. A. Anderson and W. T. Whelan, J. Soc. Motion Picture 
Engrs. 50, 199-208 (1948). 
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SPHERE DIAMETER (INCHES) 


Fic. 2. “Critical” velocities which mark the transition between 
entries of steel spheres with and without cavities. 


Similar pairs of spheres were dropped into water at 
30 and 45 ft/sec. At 30 ft/sec it was found that a cavity 
was always produced but that the cavity was generally 
smaller when the sphere was clean. At 45 ft/sec no 
difference in the entries was observable. 

In a further investigation of this effect observation 
was directed primarily to the question, whether or not 
a cavity was formed on the entry of the sphere. It was 
found, however, that this distinction is not always 
definite ; there are intermediate cases where the cavities 
are very small or consist only of a small cloud of tiny 
bubbles. 

& Records were obtained photographically at 64 
frames/sec by means of a modified Bell and Howell 
camera. The camera was modified by replacing the usual 
rotating shutter having an aperture of 135 degrees, 
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with one of about 3 degrees. This shortened the ex. 
posure to about 1/7800 second. 


EFFECT OF VARIOUS SURFACE CONTAMINANTS 


Since the cavity formed by relatively high speeq 
water entry is practically independent of the surface 
condition, it would appear that the cavity produced by 
contaminated spheres at low speeds should be regarded 
as an extrapolation of the high speed results, rather 
than the cavityless behavior of clean spheres. It js 
obvious, however, that handling is not a very satis. 
factory way of achieving the end, especially since the 
skin condition of the technician may vary greatly, 
Accordingly tests were carried out to determine whether 
some other surface preparation might be used. One-inch 
and $-inch steel spheres were dropped into water from 
various heights with the following surface contaminants: 
handling, petroleum grease, whale oil, lanolin, hexa- 
decane, hexadecane with 10 percent oleic acid, a baked 
silicone coating, and handled silicone. 

The petroleum grease was a grade provided for tech- 
nical purposes. Hexadecane is an oil that will not 
spread on water unless contaminated, for example, 
with oleic acid. These greases and oils were applied to 
the spheres in a very thin layer by wiping with a cloth 
impregnated with the substance. The silicone coating 
was prepared by placing spheres in polymethylsiloxane, 
DC-500, 70 centistokes, for two days at 150° to 160°C. 
A very thin but rather durable yellowish coating results, 

Handled or greased spheres produced cavities at 
lower speeds than did clean silicone-coated spheres. 
Little difference exists apparently between the effects 
of the various other contaminants in producing cavities. 
None was found to give a more hydrophobic coating 
than handling; in fact there was some evidence that 
handled spheres were most likely to produce cavities. 


SCALING OF THE WATER ENTRY OF STEEL SPHERES 


Even when high speed motion pictures were used it 
was difficult to decide when an entry should be classified 
as “rough.” In addition, the entries did not repeat well. 
Two entries under presumably identical conditions 
would result in no cavity in one trial, and in a small or 
large one when the test was repeated. To facilitate the 
location of the transition velocity the maximum vd- 
umes of the cavities formed for various entry speeds 
were measured on the motion pictures and were plotted 
against speed separately for the different sphere sizés. 

Various steel balls were used in the test, with d 
ameters between } and 3 inches. All were commercial 
ball bearings. Their surfaces were reasonably perfect, 
but there was evidence of stain marks on some whid 
could conceivably effect the experimental results. 

Clean and greased spheres were used. The forme 
were cleaned in ether and were used reasonably sow 
after cleaning to prevent contamination or accumul- 
tion of dust. They were dropped into the water fro 
clean forceps. The greased spheres were wiped with! 
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CONDITION OF A SPHERE ON 


cheese-cloth rag impregnated with petroleum grease 
and were held with this rag prior to dropping. The 
dropping heights were so chosen that each change in 
height increased the entry velocity by 1 ft/sec. 

The volumes obtained for the entries of clean and 
greased spheres are listed in Tables I and II, respec- 
tively. Two or three readings are given for many condi- 
tions and the agreement is frequently poor. These 
volumes were plotted against entry velocity and graphs 
were drawn for each sphere size. The velocities corre- 
sponding to the intersections of these graphs with the 
velocity axis have been plotted against sphere diameter 
in Fig. 2. These velocities may be regarded as “‘critical”’ 
velocities of the transition between entries with and 
without cavity. 

The data of Tables I and II show that the critical 
velocity is higher for clean spheres than for greased 
ones. Worthington! obtained a transition velocity of 
nearly 19 ft/sec after careful polishing of a $-inch steel 
sphere. Figure 2 would indicate a value slightly greater 
than 20 ft/sec. This is satisfactory agreement in view 
of the different methods used in preparing the surface. 
Worthington suggested dust as the cause of early sepa- 
ration of water from the sphere, and regarded the transi- 
tion velocity of more than 23 ft/sec obtained by drop- 
ping the sphere through a flame, as substantiating this. 
The behavior might also be explained by the presence 
of gas nuclei on the surface. 

It should be noted that in no part of our investiga- 
tion were extreme precautions taken to obtain per- 
fectly clean surfaces. The solvents used in cleaning the 
surface were probably not entirely grease free. The 
precautions required to produce perfectly clean surfaces 
are inconsistent with the conditions under which weapon 
models are tested. 

From the relation between the critical velocity and 
the sphere size it should be possible to determine the 
similarity laws which govern this transition. For water- 
entry behavior, since a free surface is involved, one 
may expect dependence on the Froude number, F 
=1"/lg, where, for the problem at hand, » is the entry 
velocity, /, the diameter of the sphere, and g, the ac- 
celeration of gravity. Furthermore, the factor that 


TaBLE I. Maximum volume of cavity (in sphere volumes) 
formed by clean steel spheres. 











Entry speed Sphere diameter in inches 

(ft/sec) F 4 1 3/2 2 3 
25 21 
24 37 13 il 
23 22 13 17 22 20, 20 
22 13,2 32, 16 1,0 14 21 24 
21 } ye HS ey 4 21, 20 19,8 
20 4 Zz 28,0 7,0 8, 13,0 
19 4,0 0,0 17,0 12,0,0 25, 18, 6 4 
18 1,0 0 2,0 1,0 Y 0 
17 T,0 1,0 T,0O , 0 
16 0,0 T,0O T,0O 0 0 
15 1 T 0 
14 0 0, 0,0 0 
13 T 
12 0 








T =trace (cavity less than 0.2 sphere volume). 


ITS WATER-ENTRY CAVITY 1221 
determines whether a cavity will be formed is whether 
the water clings to the surface of the sphere as the 
sphere passes into the water. This will apparently result 
in a dependence on surface-tension forces and hence on 
Weber’s number, W=plv?/T, where p is the density of 
the water and T is its surface tension. In the expres- 
sions for these two quantities only the velocity and 
diameter were varied. 

The behavior shown in Fig. 2 for greased spheres 
appears to be approximately that of the Froude cri- 
terion for diameters up to about ? inch, but for larger 
diameters the transition velocity is not correspondingly 
high. 

At high entry velocities cavities may always be ex- 
pected to accompany water entry; the initiation of the 
cavity when the missile first enters the water may or 
may not be related to the initiation of cavitation bub- 
bles as observed in steady-state motion in a water 
tunnel. However, it is of interest to compare the limit- 
ing conditions under which the water-entry cavity and 


TaBLeE II. Maximum volume of cavity (in sphere volumes) 
formed by greased steel spheres. 











— Sphere diameter in inches 

(ft/sec) } a J 1 3/2 2 3 
21 55 18 
20 48 
19 33 38 22 17 
18 48 33 20 
17 34 30 24 17,6 
16 47 20, 18 25, 11 24, 23,8 14,6 7 
15 145 9 37 2, 3. 7 11, 6,1 53, § 12,6 
14 38 3 41,6,1 19,7, T 16, 11, 4 31,1 13,8 
13 86 51, 41 6, 4 4,2,2 4 5,7 
12 27 5 | y 7,1 3,0 0 Yi 
11 41,18 10 27,2 x T,0 T 
10 14 17 T,0O T 
9 32 18 17 
g 23 12 y Oe 
7 15, 12 0 
6 2 6 
5 5 0 
4 1 0 
3 1 1 








cavitation bubbles are produced. The velocity at which 
cavitation might be expected to occur may be esti- 
mated from experimental data obtained from water- 
tunnel investigations on cylinders with clean hemi- 
spherical noses.’ Such data indicate that, for atmos- 
pheric pressure, cavitation first becomes visible at 
about 52 ft/sec. Figure 2 indicates that the critical ve- 
locities for greased spheres do not exceed one-fourth of 
this value and it appears improbable therefore that the 
limiting conditions for initiation of the water-entry 
cavity are the same as those for the production of cavi- 
tation bubbles. 

There is no clear evidence of Weber scaling in the 
data of Fig. 2. In going from }-inch to 2-inch spheres 
the Weber number would remain constant if the ve- 
locity were reduced to half its initial value. The only 
velocity decrease is found in the data for clean spheres 

3H. Rouse and J. S. McNown, Cavitation and Pressure Distribu- 


tion (Head Forms at Zero Angle of Yaw) (State University of Iowa 
Studies in Engineering, Iowa City, Iowa), Bulletin No. 32. 














1222 


where the change of value in going from }- to 2-inch 
spheres is less than fifteen percent. 


EFFECT OF SPHERE DENSITY 


An investigation was made of the effect of the de- 
celeration of the sphere, on cavity formation, by drop- 
ping greased spheres of various densities into water. 
The substances used and their specific gravities are as 
follows: tungsten, 16.8; steel, 7.8; Duralumin, 2.8; 
magnesium, 1.8; polystyrene, 1.1. 

Data obtained for the entry of four sizes of Duralumin 
spheres are given in Table III. All transition velocities 


appear to be very close to 12 ft/sec. For the five ma- . 


terials of different densities, 1-inch spheres gave transi- 
tion velocities between 12 and 13 ft/sec, and }3-inch 
spheres between 10 and 12 ft/sec. 


EFFECT OF LIQUID COATING 


Worthington! and Bell‘ investigated the entry of a 
sphere into a liquid after the sphere had been dipped 
into the same or another liquid. They found that the 


TasLe III. Maximum volume of cavity (in sphere volumes) 
formed by greased Duralumin spheres. 








Entry speed Sphere diameter in inches 





(ft/sec) 4 1 3/2 3 
20 60 3 
19 
18 23 22 27 
17 57 25 
16 45 23 28 11 
15 61 14 8 
14 1 4 10 10 
13 + 8 8 3 
12 0 0 0 0 
11 1 2 0 
10 0 0 








transition between entries with and without cavities 
occurred at much lower entry speeds than when the 
spheres were clean. 

Table IV lists five liquids of varied physical proper- 
ties into which spheres were dipped in investigating 
this behavior. After immersion in one of these liquids, 
1-inch steel spheres were dropped into water at speeds 
of 5, 10, and 20 ft/sec. The maximum sizes of the cavi- 
ties produced were determined from high speed motion 
pictures. 

For entries at 20 ft/sec large cavities of almost the 
same size (50 sphere volumes) were produced with the 


*G. E. Bell, Phil. Mag. 48, 753 (1924). 
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TaBLeE IV. Liquids into which sphere was dipped. 

















Kinematic Mass ad. 

atk ee a 
Liquid pte a (dynes/cm) _ stokes) = 
Water 1.00 72 1 0.15 
Methyl “Carbitol’”* 1.02 41 4 0.17 
Aerosol OT» 1.03 27 23 0.24 
Propylene glycol 1.04 40 45 0.39 
Glycerol 1.26 63 322 0.72 








* Diethylene glycol monomethyl ether (Carbide and Carbon Chemicals 
Corporation). 

b Aerosol OT Clear 25 percent. A wetting agent produced by the Ameri. 
can Cyamid and Chemical Corporation. 


three most viscous liquids, while cavities of about ten 
sphere volumes were produced with liquids of low vis- 
cosity. When the entry speed was reduced to 10 ft/sec 
no cavities at all were formed when the low viscosity 
liquids were used and cavities of one to two sphere 
volumes were formed with the three of highest vis- 
cosity. Even at 5 ft/sec small cavities were observed 
for the aerosol and glycerol. 

Since the formation of a cavity at low speeds after 
dipping in a liquid may be caused by the shearability 
of the liquid layer on the sphere surface as suggested 
by Worthington,! an investigation was made as to the 
amount of liquid that might be expected to adhere to 
the sphere in the case of the various liquids. One-inch 
steel spheres suspended from a pointed electromagnet 
were immersed by bringing a container of liquid up 
from below. Then the container of liquid was removed, 
the sphere was held in place for 5 seconds to allow any 
easily detachable drop to fall off, and then the sphere 
was dropped on the pan of a scale. 

The amount of liquid which was conveyed to the 
scale pan with the 1-inch sphere is given in Table IV 
for the various liquids. It is seen that these masses are 
in good correlation with the viscosities of the liquids. 
The correlation with specific gravity is probably for- 
tuitous. 

The author wishes to express his indebtedness to 
several members of the Hydroballistics Division; to 
Dr. J. H. McMillen for valuable discussions and sug- 
gestions throughout the progress of the work; and to 
Miss J. C. Woodhull, Mr. D. J. Milano, Mr. N. R. 
Corngold, and Mr. P. M. Aronson for their contribu- 
tions to this project. He wishes also to thank Professor 
G. Birkhoff of Harvard University and Dr. W. A. 
Zisman of the Naval Research Laboratory for helpful 
discussions. 
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Numerical Computation of Neutron Distribution and Critical Size* 


A. S. THomPson 
Atomic Energy Research Department, North American Aviation, Inc., Downey, California 


(Received December 4, 1950) 


A numerical procedure is presented for the calculation, by diffusion theory, of the critical size of a nuclear 
reactor, and of the time-dependent distribution of neutrons in the reactor as a function of position (one- 
dimensional space) and neutron energy. The partial differential equations for diffusion theory are approxi- 
mated by a system of ordinary differential equations, for which a stepwise convergent solution for neutron 
distribution is obtained. The critical size of the reactor is found from a comparison of the magnitude of the 
slowing down densities of the thermal and fission energy limits. The procedure is applied, by way of example, 
to a reflected sphere. 

This method for calculating neutron distribution is thought to be new in that an integration technique is 
developed and used for the solution instead of differentiation. Primarily because of this, the convergence to 
the correct solution is rapid. In a one-dimensional space all physical quantities can be considered as either 
continuous or discontinuous functions of both position and eriergy. Such items as fission and nonproductive 
capture in the slowing down range can be included as part of a standardized procedure. Certain transient 


a, he... oe 


i ieee 





and steady-state perturbation like motion of control rods can be studied. 





I, INTRODUCTION 


EVERAL methods are available for the calculation 

of the distribution of neutrons and the critical size 
of nuclear reactors."* The simplest of these assumes 
that all neutrons have the same energy. Since a simple 
diffusion equation is applicable, analytical solutions can 
be obtained for a wide variety of problems in a simple 
and straightforward manner. Since, however, the slow- 
ing down of neutrons is ignored, the results are un- 
realistic. A more sophisticated approach to the problem 
is made with ‘‘age”’ theory. For certain idealized cases, 
namely, bare homogeneous reactors with no absorption 
during slowing down of the neutrons, exact analytical 
solutions are easily obtainable. For most cases of re- 
actor design the idealizations of age theory also are 
unrealistic. In multigroup theory the neutrons are 
assumed to exist at discrete energy levels, the reaction 
among neutrons in various energy levels being related 
by boundary conditions. The larger the number of 
groups involved, the better is the accuracy of the 
calculations and, unfortunately, the greater is their 
complexity. In order to get analytic solutions it is 
necessary in multigroup theory to consider the re- 
actors homogeneous and uniform over the region of 
space under consideration. Calculations of critical size 
and neutron density are, in general, complex and 
tedious, the accuracy of the calculations being limited 
by the complexity of the calculations. 

In this paper a numerical procedure is presented for 
calculation, in accordance with diffusion theory, of the 
critical size of a nuclear reactor, and of the time-de- 
pendent neutron distribution in the reactor as a func- 


* This report is based upon studies conducted for the AEC under 
Contract AT-11-1-Gen.-8. 
‘Clark Goodman, editor, The Science and Engineering of Nuclear 


in (Addison-Wesley Press, Inc., Cambridge, Massachusetts, 
7). 


Harry Soodak and Edward C. Campbell, “Elementary pile 
i AECD-2201, Clinton National Laboratories, August 4, 


tion of position (one-dimensional space) and energy. 
The procedure is applied to a reflected sphere, which 
is symmetrical about its center, in which the reaction 
of the neutron with its environment is assumed to 
change both with radial position in the sphere and with 
neutron energy. The partial differential equations of 
diffusion theory are approximated by a system of 
ordinary differential equations, for which a stepwise 
convergent solution for neutron distribution is ob- 
tained. The critical size of the reactor is found from a 
comparison of the magnitudes of the slowing down 
densities at the thermal and fission energy limits for 
the steady-state solution. 

Certain advantages are inherent in the numerical 
procedure, as opposed to the various analytical pro- 
cedures which are customarily used. For a problem of 
the distribution of neutrons in one-dimensional space, 
all the physical variables can be considered as either 
continuous or discontinuous functions of both position 
and energy with little increase in complexity. Such 
items as the fast fission effect and absorption in the 
slowing down range can be included as part of a stand- 
ardized procedure. The numerical solution can be used 
to study the effects of certain transient and steady- 
state perturbations, for instance, motion of control 
rods, on the over-all distribution. The one disadvantage, 
as compared with analytical methods where these are 
applicable, is that all calculations must be repeated 
each time a new type of reactor is considered, since 
parametric relationships are, in general, obvious only 
for geometrically similar reactors. The numerical pro- 
cedure presented here is limited to those cases of neu- 
tron distribution to which diffusion theory is applicable.’ 


Il. MATHEMATICAL DEVELOPMENT 


The mathematical development of the numerical 
method presented can be divided into three general 


3 Robert E. Marshak, Revs. Modern Phys. 19, 185-238 (1947). 
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parts. First, the diffusion equations for the thermal and 
slowing down ranges are stated in a dimensionless form 
suitable for numerical computations. Second, the par- 
tial differential equation for the slowing down range is 
converted to a system of ordinary differential equations. 
Third, the ordinary differential equations for both the 
slowing down and thermal ranges are written in a 
standardized form and integrated. The integrated equa- 
tions obtained are the basis for the numerical com- 
putations. 


A. Statement of Diffusion Equations 


It is assumed that the slowing down process can be 
described by the following partial differential equation, 
adapted from gaseous diffusion theory,'? valid for 
thermal reactors: 


1 1 O(nvt >>.) nv 
-V-—Vnv— nv > == —-——_+—_,, 


(1)* 
wr dlogE vd 


where >>,=transport macroscopic cross section for 
neutrons at energy E, n=density of neutrons per unit 
logarithmic energy range at energy E, v=velocity of 
neutrons at energy E, >°.= total macroscopic cross sec- 
tion of matrix for absorption of neutrons at energy E, 
>..= macroscopic cross section of matrix for scattering 
of neutrons at energy E, =mean fractional loss of 
neutron energy per collision with scattering matrix, 
and ¢= time. 

Equation (1) is subject to two spatial boundary con- 
ditions which are determined by the geometry of the 
reactor. It is necessary that the neutron densities at 
the fission and thermal energy limits be tied together 
by the following two simultaneous equations for the 
thermal range: 


1 1 O(nv) th 
—~V-—V (nv) n— (22) th Dat (nv):& = » & 
3 > vol 





t 


and 


(nv) o& Le -_ 


Bas 


n(n) Xyd(logE), (3) 


Eth 


n(nv) en yt 


where n»=number of neutrons produced per fission, 
(nv)yn= thermal neutron flux, (v)o>=neutron flux at 
fission energy boundary, (v),=neutron flux at thermal 
energy boundary of slowing down range, and }°, 
=macroscopic cross section for fission at energy E. 
Equations (1) through (3) can be rewritten in dimen- 
sionless form as follows: 


V-aVo(a, B, t) —b¢(a, B, t) 
OL tcdh(a, B,t)] AALEcdH(a, B, t) ] 
= + 


0g > nkcdvdt 


* By a generalization of the terms 2, and 2,, Eq. (1) can be 
applied to reactors which are not thermal without added numerical 
complication. 





(4) 
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V: aV¥(a, t) - dy (a, t) +ridd(a, 1, t) 
=N(1/Cas) Lave, 1)/at], (5 





(a, 0, t)&d=mp(a, t)e+ J ee dB, (6) 
where 
a=Don/ds, 
b=)06/Lin, 
c= 1/log(Eriss/Etn), 
d=) ./Xiu, 
e=2)s/Xi4, 
A=3(R Yn)’, 
log(Etise/E) 
log(Etiss/Etn) 


R=critical dimension of reactor, 


>= transport macroscopic cross section for refer- 
ence material at thermal energy, 


$(a, B, t)=nv/(nv)thmax, 
¥(q, t) = (nv) tn/(12)thmax, 
(nv)thmax= Maximum thermal flux, 
a= dimensionless space coordinate ; for a sphere, 
for instance, da=dx/ax*, where x=r/R,1 


being the variable radius, and R the critical 
radius. 


B. Conversion of Partial Differential Equation to 
System of Ordinary Equations 


The two terms on the right of Eq. (4) represent the 
total derivative of neutron density with respect to 
change in energy level during a transient slowing down 
process. This can be demonstrated as follows. The total 
derivative of neutron density with respect to energy 
is given by 


dLécdo(a, B, t)] aL &cdo(a, B, t) ] : dL Ecdo(a, B, t)] dt 
dB aB ™ at dB 





The rate of slowing down is given by 
dB/dt=> nédcv. 
Hence, the total derivative is given by 
dLécdo(a, 8, t)] ALEcdo(a, B,t)] AL Ecdo(a, B, t)] 
aa ’ 
dg ag Y ntcdvat 


which is equal to the right side of Eq. (4) divided by) 

The partial differential with respect to 6 on the 
right of Eq. (4) can be written as the limit of a di 
ference, 


dLécdo(a, B, t)] L Ecd(a, B, t) — tcdo(a, B—A§, t) 
op AB-+0 Ag 











’ 





In 


“fer. 


rere, 
R, 1 
tical 


n to 


t the 
t to 
lown 
total 
ergy 


| dt 


dB 
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and approximated for the interval A8=k by 
dLtcdo(a, B,t)] tcdo(a, B, t)—Ecdo(a, B—k, t) 
ap k 


The partial differential with respect to ¢ in Eq. (4) can 
be written 


ddLEcdo(a, B, t)] pdrEcdp(a, B—k, 1) 
> nécdvdt k 





(7) 





, (8) 


where u is a measure of the transient variation in neu- 
tron density during slowing down, given by 


k 1 AL Ecd(a, B, t)] 
""Snteds tedola,B—k,t) at 


Substituting Eqs. (7) and (8) into (4) gives an approxi- 
mate expression for the neutron density in the slowing 
down range at 8= m in terms of the density at B= m—k: 





écd 
V-aVo(a, m, )-r( 14+ )o(a, m, t) 


Ecd 
= 1+ 1) molar m—k,t). (9) 


In Eq. (5) the time for absorption can be described 
> nvAt=> > s/>- a= 1/0. 


The change in the quantity ¥(a) during the interval 
At can be written 


Ay(a, t)= ai vido(a, 1, t)/b, 
v= (—b/td)[Ay(a, t)/$(a, 1, t)]. 


Equation (5) can hence be written 


V-aVy(a, t)—dbp (a, t)= —X(1+ r)tdd(a, 1, t). 


where 


(10) 


C. Integration of Equations 


An integrated form will now be developed for Eqs. 
(9), (5), and (6) for one-dimensional problems, that is, 
where all space variations can be expressed as functions 
of one coordinate. For such one-dimensional problems 


V-aV=(1/x2)(d/dx)ax%(d/dx), 


where the exponent g has the following values for 
various coordinate systems: 


Cartesian 0, 
Cylindrical polar 1, 
Spherical polar 2, 


and x is the coordinate distance divided by the critical 
dimension R. Multiplying both sides of Eqs. (9) and 
(5) by ax2¢ gives for Eq. (9) 


{DP?—g,*(a)} d(a, mM, t)=h,(a)(a, m—k, t), (11) 
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where 
D=d/da=ax%(d/dx), 
gi(a) = Lax?*r(b+ Ecd/k) }}, 
hy(a) = —ax**(1+p)écd/k, 
and for Eq. (5) 
{ DP— g(a) }y(a, t)=he(a) (a, 1, 4), (12) 


where 
g2(a)=[ax?nb |}, 


h(a) = — ax?9)(1+- v)&d. 


Equations (11) and (12) are of the type, for a given 
time, ¢, 


{ D’— g*(a)} ¥(a)= X(a). (13) 


It can be shown by direct substitution that Eq. (13) is 
factored‘ into 


{D+ f(a)}{D—f(a)} ¥(a)=X(a), (14) 
where f(a) must satisfy the nonlinear equation 
Df(a)+f?(«)— g’(a)=0. (15) 
Let 
{D—f(«)} ¥(a)= U(a). (16) 
Then 
[D+f(e)}U(e)=X(a). (17) 


Equation (16) is integrated to give 


V(a)=exp(+ f fdas) 
| f e0(- f ites urardarta, (18) 


and from Eq. (17) 


v(a)=exp(— f jaar) 
x| J ‘ex (+ J. fdas )X (ca)dart-c. (19) 


Combining Eqs. (18) and (19) to eliminate U(a), 


v(a)=exp(+ f'fi)} f exo(—2f sao) 
x| J “exp (+ J “fdax) X (az2)dacg 
| +e1)tarter. (20) 


‘E. L. Ince, Ordinary Differential Equations (Dover Publica- 
tions, New York, 1944), pp. 23, 128. 
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Equation (20), subject to Eq. (15), is the basis for a 
numerical solution to Eq. (10) for the slowing down 
range and to Eq. (12) for the thermal range of neutron 
distribution. These equations, in conjunction with Eqs. 
(6) and (15) are sufficient for a complete description 
of the neutron flux density as a function of position, 
time and energy. Numerical forms are needed for the 
solution of Eqs. (15) and (20). It is to be noted that an 
analytical solution to these equations is available only 
for a few very special cases where all physical quan- 
tities are constant throughout both space and energy 
variation. For a general case Eq. (15) can be written 
for numerical solutions 


fla)=g0)— [ tf@)-e@)\da, (21) 
0 
Equation (21) is to be solved by stepwise integration. 


Ill. PROPERTIES OF EQUATIONS 
A. Transient Solution 


Equations (9) and (10) indicate that the effect of 
time variation on neutron density is to multiply the 
solution for each step on the slowing down range by a 
factor (1+) and the solution in the thermal range 
by a factor (1+¥v), where uw and » represent the per- 
centage change in neutron density due to the transient. 
The numerical procedure yields the transient solution 
when values different from zero are assigned to » and 
v. A steady-state solution is obtained as a special case 
when yu and » are assigned the value zero. The total 
percentage change in neutron density during one neu- 
tron cycle, consisting of m steps in the slowing down 
range and then the thermal range, is given by the 
product, 


14+2= (1+) TI (1+1)m. (22) 


The time for one neutron cycle is given by 


n+1 1 n k 
" S, wn 5 se 5 alk = 


B. Discontinuities 








It is. often expedient to treat some of the physical 
variables in a reactor as discontinuous; for instance, at 
a junction of two different materials. For these cases 
the function g(a) is discontinuous, and at the discon- 
tinuities the integrals of Eqs. (20) and (21) are indeter- 
minate. These quantities are to be determined from the 
conditions that both the flow of neutrons and their 
density must be continuous across the boundary. 
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These conditions are stated 


Y =continuous, 


DY =continuous. 
From Eq. (16), 


DY (a)=U(a)+f(a)¥ (2). 


It can be seen that if f(a) is continuous, both Y and Dy 
are also continuous, as is also the solution to Eq. (21), 
the indeterminacies disappearing in this case. 


C. Properties of the Nonlinear Equation 
(df/de)+f?—g’=0 


This equation is a special form of the Riccati equa- 
tion for which the complete solution can be obtained 
in terms of a known particular solution.‘ So far as the 
author knows, the particular solution cannot be ob- 
tained analytically if g is completely arbitrary. How- 
ever, some useful information concerning the properties 
of the equation can be obtained for simple forms of g 
for which an analytical solution is obtainable. The 
simplest case is g equal to a constant, for which the 
particular solution is 


fi= (24) 
The complete solution to the equation is of the form 
f=gi—1/z. (25) 
Substituting this value for f into the original equation 
(dz/da)—2gz+1=0, (26) 
for which the solution is 


z= (1/2g:)[1+exp2g:(a+<c) ]. (27) 


The complete solution for f for the case where g=g, is 
hence 


f/gi=exp[2gi(a+c)—1]/exp[2gi(a+c)+1] 
=tanh [gi(a+c)]. (28) 


The value of f/g; is asymptotic to unity for large values 
of gi(a+c). The slope of the curve is unity for gi(a+<¢) 
equal to zero and is asyniptotic to zero for large values 
of gi(a+c). This equation enables us to study the 
effects of a discontinuity in the value of g. Assume that 
f/g has the value unity for values — © <a<0 and that 
g has the value go for values — © <a<0O and 4g; for 





DY 
1), 


jua- 
ned 
the 


Ow- 
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of ¢ 
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(27) 


£1 1S 


(28) 


alues 
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values 0<a<+. Then a plot of f(a) is as shown 
schematically in Fig. 1, where the value of c has been 
so chosen that f(a) is continuous at zero. 

As can be seen by inspection of the original differ- 
ential equation, it is characteristic of the equation that 
the greater the difference between g’ and /?, the larger 
is the slope of the function f in such a direction as to 
eradicate the difference between the two functions. 


IV. PROCEDURE 


The numerical procedure for finding the neutron dis- 
tribution and critical size is as outlined below: 

1. Guess a curve for the thermal neutron distribu- 
tion ¥:(a), and at the same time a value for the constant 
\1, which is a measure of the critical size of the reactor. 

2. Solve Eq. (6) for the distribution of fast neutrons 
¢(a, 0) (neglecting for the first trial the fast effect as 
represented by the integral). 

3. Solve Eq. (21) for f(a) for all energy levels of 
interest, including thermal. 

4. Solve Eqs. (10), (20) in succession for the neutron 
distribution at the energy levels of interest, ¢(a, 0.1)— 
(a, 1.0), using ¢(a, 0) from Step 2 as a starting point 
and getting finally (a, 1.0). 

5. Solve Eqs. (12) and (20) for the thermal neutron 
distribution (a). 

6. Compare curves y¥2(a) and ¥(a), obtaining an 
estimated value of the constant a from the equation 


1+2= J Walo)da / J ‘Yala)da. (29) 


7. Estimate a new value of the critical constant » 
from the equation 


do/A= 1+ 2. 














Fic. 2. Reactor configuration. 
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Fic. 3. Thermal neutron distribution for various trial \-values. 


8. Repeat Steps 1 to 7, using Y2 and dA, as a starting 
point. 

9. Plot (1+) as a function of \ for these trials, 
extrapolating to find the value of \ for which (1+2)=1. 
Continue the process until a value of \ is found which 
checks the value (1+2)=1. This is steady-state solu- 
tion, since 2 measures the percentage change in neutron 
density during one neutron cycle. 


V. NUMERICAL EXAMPLE 
A. General 


For demonstration purposes a simple numerical ex- 
ample for a spherical reactor will be worked.t The 
hypothetical reactor considered is primarily composed 
of graphite as a structural and moderating material. 
It is assumed that the fuel material is separated 
uranium-235 uniformly distributed (1.5 mg per cm’) 
throughout a central portion, or core, of the reactor. 
The core for structural purposes is encased in an 
aluminum can, which for cooling has an air gap at its 
outside. A plug in the center of the core is composed of 
graphite which carries no fissionable material and a 


t The values of constants used are not the best available but 
are those cited for illustrative purposes in reference 1. 
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TABLE I. Computation layout for the critical solution. 


® ® © 
Bx? f; a 


3obb 


a 


0.9745 


0.10 | 0.9338 10% 


x 107% 
x 1072 
0.25 6.25 x 10 


0.8126 6.363 


x 10~ 





x 10° 


138 x 10° 


942 x 10% 


536 x 10° 


x 10° 


576 x 109 





8 


931x107? 


075x10~ 


238x10~ 


766x107? 


. 588x107? 


872x10 


5232x107 
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TABLE I.—(Continued). 


@®|lelelelfletleloe ® le lea 
" . $(x, 0.2) “ - $(x, 0.4) 
©©|®e@|Lel9@|e| os |oalfelee|l/sa| oe 


0 ° 0 ry 0 0 5.64710" |5, 64721072 
' 


0.1522 


0.81! 0 -2 


0. 


1.634 |3.618x10~ 
65721072 
1. 


50€x1073 


-489x10-8 


569x10~ |4.9242107% 


-227x10 


30510719 Io, 62421075 
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TABLE I.—(Continued). 


@5) @) @ ® @ 
d(x, d(x, 0.8) 
f* 


a 4°@1G 


5 5. 3181072 


-795x107)9 Io, 534x107 
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® 


&) 


0.1919 


0.11292 


®) 
_ 


7162x1072 
1.602x107~* 


1x1073 


625x10° 

-273x10°8 
913x10~° 
632x107? 


-870x10~ 


i) 


.869x10~3 
.73@x1079 
.282-1073 
.3281073 
-115x107~4 


2026x1074 


TABLE I.—( Continued). 


@) 


Fx? 


0 
-2 
1 
0 
0.4074 
0.6620 
0.9166 
1.27% 
1.630 
2.088 
2.546 
3.106 
3.667 
4.329 
4.91 


6.518 
7. 
8.250 


9.218 
2185 


9 


49.88 
ie) 


+ 
’ 
0 


1.694 
1.773 
1.942 
2.110 
2.294 
2.477 
2.67. 

2.872 
3.084 
3.297 
3.524 

7 


G8 


f2 


0 
3.66x10~% 
1.30x1072 


4e48x107* 


09 


i) 
0.146 
0.723 
1.30 
1. 
1.9 
2.27 
2.55 
2.80 
3.06 
3.27 
3.48 
3.66 
3.83 


@) 


Var) 


0 


3.651073 


3.981072 


2.236 


2.326 


2.413 


2.498 


1, 


2 542n1U™ 
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1.17810 


8.885x10~ 
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TABLE I.—( Continued). 


&) @) @) @) 
e 1©@@/@@|f*@ 


C) ") ry r) 
32x 107? 
0.12640 | 3.16x10~3 
0.3218 
0.5173 }1.925x10~2 
0.8807 
1.2441 |6.329x1072 
1.828 
2.412 |1.547x1072 
3.27% 
4.135 |3.184x1072 
5 


7.962 
9,821x1071 


@ 
@® © |4*@ 


r) -955x1072 
0.6270 
1.2539 -642x107) 
1,515 
1.776 |8,88%x1072 
1, 
-890x1072 


735x107) 


5530x1072 


3.351x107) 


2.497x1072 
1, 
1.827107 
0.9 
‘ 1.330x1072 
0.1738 
0.1714 
0 


t 


0 
0.7 1. 


0.6914 

0.6342 {1.148107 
0.5592 

0.4843 6891072 
0.4332 

0.3820 523x107 
0. 

0.3026 8121072 
0.2708 

0.2390 |3.45810 
0.2140 

0.1890 |2.3¢8x107* 
0.1681 

0.1472 |1.547x1072 
0.1306 

0.1141 94x1073 
0.1010 

-805x1072 |3.890x1073 
-780x1072 
756x10~2 





@ @) @) 


Y> (x) Resulting 


¥, (x) Y Shape 
@) v2 (x) 


0.9955 1.004 
0.9680 
0.9248 
0.8916 
0.8520 


0.8046 
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TABLE I.—(Continued). 























































































































REACTOR CONSTANTS 
MEDIUM PLUG AND CORE ALUMINUM | AIR GAP 
REFLECTOR 
Range tt Seas 0.05<x 0.50<x 0.505 sx 
e x<=l] x<40.50 x £0,505 x 50,525 
d 6.00 x 10? — 
k A 5 SS ence enee 
c 2. 2 ee 
|) 
t 3.68 x 1072 [3.68 x 10-11 7.80 x 107 0 
a = 2t)/2¢ 1 1 472 o 
*s,  .60x10%p.50 x 1073|1.20x 10%] 0 
Z, _|3.90 x 10-1 3,90 x 10-1 | 8.00 x 10% 
"Zr 0 2.10 x 1073 0 
> = D,/ Lah.77 x 10 6.79 x 1073 | 3.26 x 10% 
azle/ Zyl 1.6 1.06 2.17 x 107 
te = Ly/ Dp 0 5.70 x 1073 0 
7 0 2.30 0 
€ n.58 x 1072 hse x 1072 | 7.00 x 10% 
D C) 7.23 x 1072 0 
5 b.76 x 107 £.76 x 107 | 2.50 x 102 
F 5.8 .07 x 102 | 1.96 x 107 
G n.005 x 107.005 x 1031 9.14 x 101 
DEFINING EQUATIONS 
= ne/Ed [See Eq. (6)] 





-hy(x)/ax? = g)*(x)/ax® = (X€ ed/k)x? 





See Ea. (11) 








= go?(x)/ax* 


= br x? 


G S-ho(x)/ax? = d€ dx? 


See Eq. (12) 





Critical Trial Solution Given By \ = 6000 


"I, =b= 2,=e=0 for the slowing down range since 

absorption and fission are neglec- 
ted for this range. 
**The method of calculating "f" is given by Eq. (38). 
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reflector is provided at the outside of the core which also 
is composed of graphite with no fissionable material. 
This reactor is shown schematically in Fig. 2. 

It is desired to find the critical size, and the resulting 
distributions of neutrons at various energy levels, for 
a spherical reactor geometrically similar to that of 
Fig. 2. The values of various physical parameters 
needed for the solution are taken from reference 1, and 
are given in Table I. 

In Fig. 3 are shown the distributions of thermal neu- 
trons for three trial values of \. The value \=6000, 
representing a reactor of diameter 95.6 in., gives the 
steady-state distribution of neutrons. The other curves 
for A=5000 and A=3000 give transient solutions for 
neutron distribution in subcritical reactors of the same 
geometrical configuration. Figure 4 gives the steady 
state distribution (A=6000) of neutrons at six energy 
levels in the slowing down range. 

From Fig. 3 it can be seen that the distribution of 
thermal neutrons is realtively independent of the as- 
sumed value of A. The convergence to an acceptable 
distribution is very rapid even when the originally 
assumed distribution fails radically to satisfy boundary 
conditions. This is demonstrated in Fig. 5 where, for 
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Fic. 4. Steady-state slowing down range neutron distribution 
at 6 energy levels (A= 6000). 
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Fic. 5. Thermal neutron distribution convergence through 
numerical procedure application. 


comparison, the curves are shown for a purposely bad 
initial assumption, the distribution resulting from one 
application of the numerical procedure to this initial 
distribution, and the correct distribution as obtained 
in Fig. 3. , 

In Fig. 6 is shown a plot of the number (1+) for 
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Fic. 6. Departure of neutron distribution from steady state for 
trial values of A(1+2=yW2/y;). 
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the three values of \ used in Fig. 3. This plot was useq 
in obtaining the correct value of A for a steady state 
distribution of neutrons. As has been pointed oyt 
previously, values of (1+) different from unity repre. 
sent transient states of neutron distribution. 

Time for neutron cycle: The time for one neutron 
cycle can be calculated from Eq. (23), which is 


i=E a = HE Ee cl 


1 k ofl 
= + -| f tant ), 
5 =| Sk a * | (for constant ) 


Taking the thermal and fission energy levels at 1/40 ey 
and 2 Mev, respectively, and solving for the slowing 
down energy levels by the formula, 


Ey Ey 
E,=—————_= — E=next higher energy level 
[E;/Eu]*? 38 











(which arises from the definitions of 8 and k), we may 
solve for the other energy levels and all associated 
belocities which are given below: 


E (ev) v (cm/sec) 
2 X10°(f) 1.39 10° 
5.26X 108 © 2.24 108 
1.39 108 3.65 X 107 
3.65 X 10! 5.91X 108 
9.61107 9.60X 10° 
2.5 K10-*(th) 1.55 105. 


Substituting into Eq. (23) these velocities and the 
values of the nuclear constants given in Table I, we get: 


1 
5 v 
n k - - 
a 5 | =7.35X10- sec; 


t=2.59X 10-*+-7.35 X 10-°= 2.66 X 10- sec. 





=2.59X 10 sec; 
th 





B. Discussion of Calculations 


For the present example, using 22 intervals of in- 
tegration and seven energy levels (six in the slowing 
down range and one for the thermal range), it was 
found that a computer, using an automatic calculating 
machine, required almost two eight-hour days to com- 
plete a set of calculations for one neutron cycle. The 
three trials described here, which should in general be 
sufficient for a complete steady-state calculation of @ 
given reactor, were hence done in a six-day period. It 
is to be noted that doubling either the number of energy 
intervals or the number of intervals of integration wil 
approximately double the amount of work for a com- 
plete calculation. 





Th 


fx 


ste 


qui 
spe 
to ¢ 





ed 
ite 
ut 


ron 


dev 
ving 


may 
ated 


1 the 
e get: 


of in- 
lowing 
it was 
ilating 
> com- 
e. The 
eral be 
n of 2 
iod. It 
energy 
on will 
a com 





NUMERICAL COMPUTATION OF NEUTRON 


The headings of the columns and the defining equa- 
tions of Table I, together with the following section on 
the computation of “‘f(x),” describe how Table I is 
obtained in accordance with the procedure outlined in 
this report. For tabular purposes, x is substituted for a 
according to 


da=dx/ax*. 


The trapezoidal rule is used for integration. The 
boundary conditions on Eq. (21) and Y’(x)=0 at x=0 
and Y(x)=0 at x=1. Therefore, co=0 and 


l *1 fdx2 
a= [en(-2f —) 
0 90 ax? 
#2 fdx; dx» 
| f vexn(+f —)x(o2)— fan, 
ax;? ax? 


Thus, for tabular purposes, Eq. (21) becomes 


women fi) fon(-2f' 2) 
«| f exp(+ fi = ) x(e)— fash (30) 


Computation of “f(x)”: The solution of Eq. (21) for 








f(x) is an initial value problem and is performed by 


stepwise integration. For computational purposes, this 
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equation must be written in terms of x by recalling that 
da=dx/ax*. 


Writing Eq. (20) in terms of differences in order to 
facilitate integration, we have 


z+dz 9?(x)dx i 








f(x+Ax) =se)+ J 


? 














ax? 4 ax 
as g(0)= 
g(xt+Ax) g*(x) Ax 
=f(0)+| +E 
ao(x+Ax)? a,x? J 2 
[——— ae 
ao(x+Ax)? a,x?J 2 
where a;=a(x), d2=a(x+Ax). 
Solving explicitly for f(~+ Ax), 
f(x+Ax) = —a2(x+Ax)?/Ax 
f?(x) 2f(x) a2(x+Ax)* 
+(0)e+-ax)( + 
ax(x+Ax)? Ax (Ax)? 
g°(x+- Ax) £ (x) 1 1 7\! 
ao(x+ Ax)? qua ao(x+Ax)? a,x? 


or approximately (using two terms of the binomial 
expansion of the radical), 


g(x) 





2(x+ As 
fa)+| Abi 


2a2(x-+ Ax)? 
flx-+Ax)= 


Pp 1 ! 
+ [+= |—_- — as 
2a;(x)? 2 Lao(x+Ax)? a,x? 





1+[f(x)Ax/a2(x+Ax)?] 


In Table I, f; and f2 stand for the sets of values re- 
quired in the slowing down and thermal ranges, re- 
spectively. An example of the use of data from the table 
to calculate a value for /; is as follows: 


f,(0.10)=9.53X10----, column (5), row (5), 





Ax=5X10°; 
Ip g(x+Ax) g?(x) 
“ |=s.ass -+, column (4), row (6), 
YNa(xt+Ax)? a,x? 


ao(x+Ax)?=2.25X 107: - 


neglecting 
fr x)f 1 1 


-— [as 
2 Lao(x-+Ax)? a,x? 


-, column (3), row (7), 





because of “a” generally being constant and Ax being 
small, we td 





9.53X 10-°+-4.485 X 5X 10 
9.53X 10° X 5X10 
2.25 X 10-* 
=2.63X107'- - 


f(0.15) = 





i+ 





-, column (5), row (7). 
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High accuracy of measurement in the electrolytic tank can be achieved by a suitable ratio of impedance 
as between tank and measuring circuit. With this requirement fulfilled, the axial distribution of the sym- 
metrical lens was explored and found to satisfy a relation Vze= Vo+Ae~*®), where Vo is the potential of 
the outers with respect to the middle electrode, ‘z,’ the distance from the center of the lens, while A and } are 
constants from any one geometry. 

The effect of varying the separation between the elements was specially investigated. It was found that 
the maximum value of the gradient between electrodes increases as the separation between them increases, 
reaches a flat maximum between 0.8 and 0.4 D (D, bore diameter) and then decreases with further increase 
in spacing. This means that a decrease in separation produces a weaker field and an increased focal length. 

The potential at the saddle point is very sensitive to changes in separation in the region up to about 0.5 D; 
it was also found to vary linearly with bore up to 3 D. 

The relaxation method of field plotting was also applied to map the distribution; the results are com- 





pared with those given by the tank. 





I. INTRODUCTION 


HE symmetrical electron lens has recently come 

to have an increased importance through its use 

in the electrostatic electron microscope. The search for 

an instrument which embodies the simplification thus 

introduced, for example, the less stringent voltage re- 

quirement, has focused attention on the need for a 

fuller understanding of the effects of a variation in 

parameters on the field distribution and the focal 
properties of such lenses. 

In its simplest form this type of lens consists of three 
apertured diaphragms—the two outers (usually con- 
nected together) are symmetrically disposed with re- 
spect to the middle one, which may be at cathode 
potential. It has the merit that under these conditions 
the focal length remains invariable and is independent 
of any of the usual supply variations to the outer elec- 
trodes, and, in fact, this is the way in which the system 
is used in the E.S. electron microscope. 

The published work on this system is mainly confined 
to the general pattern of field forms—Johannson and 
Scherzer' and to the extent of the aberrations, par- 
ticularly the spherical. 

An investigation by Bruck and Romani’ has provided 
information on the variations in axial field form as the 
dimensions of the central electrode were varied, keeping 
the outers and separation constant. The present work 
supplements this in a fuller way by seeking to deter- 
mine the structure of the whole field, and was directed 
also towards elucidating the effect and importance of 
variations in other parameters, e.g., separation. Atten- 
tion was also paid to the conditions at the “saddle 
point”—little has appeared in the literature about the 
sensitivity to change of this interesting singular point. 

The method adopted on the purely experimental 
side consisted in exploring the field, using the electro- 


* Now at University of Liverpool, Liverpool, England. 
1H. Johannson and O. Scherzer, Z. Phys. 80, 193 (1933). 
? H. Bruck and L. Romani, Cahiers phys. (October, 1944). 


lytic tank technique. Confirmation of these results was 
sought on the numerical side by applying the Southwell 
“relaxation” method to this study of axially symmetrical 
fields. 


Il. THE ELECTROLYTIC TANK METHOD 


The method employed in using the electrolytic tank 
to map out the potential field for conventional electron 
lenses has been described in detail by Jacob.’ 


The potential on or off the axis is picked up by a | 


probe and read directly, or it can be fed to an amplifier 
whose output in turn is applied to a cathode-ray oscillo- 
graph ; the time base of the latter is synchronized with 
the movement of the probe. The resulting trace seen 
on the fluorescent screen is the actual form of the 
potential distribution and can be calibrated in terms 
of the potential difference between the electrodes in 
the tank. 

The first alternative, however, was used. The meas- 
uring instrument was a vacuum tube voltmeter of high 
impedance (approx. 4MQ) which was connected to the 
probe, and read the potential directly, point by point, 
at constant intervals of 1 mm throughout the field. 
Many exploratory experiments were performed on 
calculable systems in order to establish measuring con- 
ditions which would result in high accuracy. 

It was finally found that faithful reproduction of a 
potential distribution required an increase in the im- 
pedance of the vacuum tube voltmeter to about 34MQ, 
by insertion of a 30MQ resistance in series with the 
probe lead. This corresponds to values of impedance 
ratio of at least 110/1 between the impedance of the 
measuring instrument and that of the circuit across 
which it is placed. With this high ratio it was then 
possible to reach the highest accuracy which could be 
read on the instrument. This requirement, while highly 
desirable, was, however, not a very convenient one, for, 


3L. Jacob, Phil. Mag. 26, 570 (1938), and Proc. Phys. So. 
(London) 63B, 75 (1950). 
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Fic. 1. Field plot of symmetrical lens. 


with this high value of series impedance, it was neces- 
sary to allow the probe to remain at each selected point 
for several seconds to insure full charging to the correct 
value. 

To expedite the time required to take a full set of 
readings it was decided to lower the series resistance to 
20MQ. Further tests showed that all readings of voltage 
would be uniformly low by 1.3 percent on this account, 


but this could readily be allowed for on the calibration 
curves. 


The over-all accuracy of measurement under these 
conditions can hence be stated as about 1.8 percent, of 
which 1.3 percent arose from the impedance variation, 
while an instrumentation error of +0.5 percent was 
ascribed to the vacuum tube voltmeter. 

In order to make a comparison between the tank and 
relaxation method of field plotting, it was decided to 
make a complete field plot for one particular case in 
which the electrodes had the following dimensions in 
the tank: outer diameter 120 mm, inner bore 20 mm, 
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thickness 1.65 mm, separation 16.5 mm. Because of the 
mirror symmetry about the axis and midplane, it was 
only necessary to explore one quarter of the whole field. 


Ill. THE POTENTIAL DISTRIBUTION 


The complete form of the field obtained in the tank 
is shown in Fig. 1. It is seen to be non-uniform and with 
the gradient smaller in the vicinity of the middle elec- 
trode than at the outer. For the “einzel” condition, the 
electrons are continually decelerated up to the “saddle 
point” and then accelerated beyond this on the way 
out of the lens. Glaser and Henneberg* have shown 
theoretically that the limiting equipotentials cross at 
the saddle point making an angle of 54°44’ with the 
axis. In the present series of measurements this angle 
was measured to be 55°, which agrees very well within 
the accuracy of measurement with the theoretical value. 
When the saddle point potential is slightly greater 
than that of the cathode, the potential on the mid- 
plane at some radial distance from the axis may become 
negative ; under this condition, the electrons which pass 
through the region near the center of the lens will be 
focussed, whereas those on the edge of the beam may 
be reflected back as in a mirror. 

An examination of the axial distribution shows that, 
in conformity with the suggestion of Scherzer,® it may 
be expressed in the form V,=Vo+Ae~*®, where Vo 
is the potential of the outer electrodes with respect to 
the middle one, “z” is the axial distance from the center 
of the lens, and A and 6 are constants and have the 
values — 10.56, and 0.005 for the geometry chosen. This 
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Fic. 2. Radial gradients at midplane and outer electrode. 


* W. Glaser and K. Henneberg, Z. tech. Phys. 16, 222 (1935). 
*O. Scherzer, Z. Physik 101, 23, 593 (1936). 
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form of function, though simple to use when substitut. 
ing directly, is difficult to integrate; hence the deriva. 
tives which are required to determine the focal length 
were obtained from the following formulas given by 
Comrie,® involving the calculus of differences: 


1 1 
Vs) =AVs—-A*V 2+—AV s+---, (1) 
6 30 
J V : 
r I = A2V 2——A'V:+---, ) 


where A"V, is the nth difference of potential. Since the 
interval dz between successive steps was small (1 mm), 
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Fic. 3. Off-axis field distribution. (Bore diam as unit.) 


it was found justifiable to neglect differences of order 
higher than A‘V,. 

To get a more complete picture of the force variation 
in the field, the radial gradients were calculated; they 
are shown in Fig. 2 for two different positions—(l 
across the plane of the outer electrode—here dV /dris 
positive and the electrons are accelerated away from 
the axis; (2) across the midplane—here 0V/dr is nege- 
tive and the acceleration is towards the axis with the 
converging effect predominating. 

The off-axis field distributions shown in Fig. 3 2 
dicate that the gradient is smaller on than off the axis, 
and that it progressively increases with distance off the 
axis. It has a maximum value which occurs at a poi! 
situated at a distance about } of the separation, awaj 


6L. J. Comrie, Interpolation and Allied Tables (His Majesty’ 
Stationery Office, London, 1936). 
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Fic. 4. Variation of axial peak gradient, and focal length, with 
separation. Figures on points give distance from outer electrodes 
at which peak occurs. 


from the middle electrode, and in nearly all cases falls 
in the vicinity of both electrodes. The slope of the 
curves increases with radial height and the fields tend 
to become more uniform between the electrodes. This 
fact is helpful when applying the “relaxation process,” 
since for the first round, the assumption of a uniform 
field now appears to have a justifiable basis. The slight 
droop at the larger off-axis distances in regions outside 
the outer electrode is to be attributed to the absence of 
a shield in this space. Actually, this would not occur in 
an electron microscope in which this system was used 
for focusing. Its effect on the focal properties of the 
lens can be shown to be quite small. 


IV. THE EFFECT OF VARIATION IN SPACING 


The axial potential distribution was determined in 
the tank for eight different spacings for the “einzel” 
condition, with a constant PD between the same set 
of electrodes as was used to obtain the complete field 
plot. Analysis of the sets of curves shows that the larger 
the separation S, the closer is the potential at the center 
of each electrode to the actual electrode potential, 
while for any fixed separation the gradient shows a flat 
maximum between the electrodes falling steeply in the 
vicinity of the central one. The maximum value of the 
gradient increases with increase in separation reaching 
a flat maximum between 0.8 and 0.4D (where D is the 
bore diameter as a unit) and then decreases with further 
increase in spacing. Figure 4 shows this variation in 
peak gradient with spacing as well as the corresponding 
effect on the focal length. The latter has been calculated 
from the weak lens approximation which is known to be 
in error by some 40 percent, but has been used through- 
out to give an indication of the general direction of the 
change when variations are imposed on the system. 
The exceedingly steep rise for the shorter spacings is 
very pronounced; but for values of S>D there is little 
change with spacing, which implies that the integrated 
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effect of the field as the lens is widened out remains 
almost constant. 

The variation of potential at the saddle point with 
variation in S is shown in Fig. 5. Here, as in the 
previous curve for the focal length, the potential is 
quite sensitive to changes in spacing at the smaller 
spacings and comparatively insensitive at values of 
S>D. This agrees with observations already made— 
that for the larger spacings little change occurs at axial 
points close to the midplane; the influence of the 
potential of the outer electrode is very marked indeed 
as the spacing is decreased. 

In an actual electron microscope using symmetrical 
lenses as elements, a shield is provided to screen the 
beam from external electrostatic effects, and also to 
provide a field-free space through which the electrons 
pass before reaching the object. To simulate the effect 
of this shield on the lens field, an open ended hemi- 
cylinder was placed in the tank on one side of the sys- 
tem only, and in contact with the outer electrode. Its 
dimensions were such that the ratio of diameters of 
shield to bore was 8/1. The plots then showed no effect 
on the form or value of the field in the inter-electrode 
space, but a rise in potential over a short distance on 
the side at which the screen had been introduced. The 
effect on the focal length when integrated over the 
whole field was small and could be neglected. 

A substantial error may arise, however, from the 
effect of the position of the object which is usually 
chosen to be in field-free space. In the present instance 
the object took the form of a semicircular copper disk 
60 mm in diameter and silver plated. It was placed at 
the following distances in front of the outer electrode 
(a) 0.2D, (b) 0.1D, (c) in contact. It was soon found 
that the potential distribution is considerably affected 
by the position of the object for small electrode sepa- 
rations. The potential at the saddle point increases as 
the object is brought nearer to the outer electrode, 
while the point at minimum potential shifts away from 
the midplane on the side remote from the object. 
For still closer positions, the gradient on the object side 
of the midplane increases rapidly and appears to become 
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Fic. 5. Saddle point potential for various separations. 





1240 J. R. 





ry) 200 
7 1 
™ ~} 
e i 
j- : 
i “| 
> ~{ 
|. =f 








on Led on ead wae one ad oa? 


os oe ome 
———} Goss Oiunarsn ( SEPARATION AS au) 


Fic. 6. Effect of bore diameter on saddle point and focal length. 


almost linear in the interval before contact is made. 
The focal length for small electrode separations (0.2D) 
and for an object distance in the range up to 0.1D has 
a value about 90 percent of that in the absence of the 
object. Outside this range the value of focal length in- 
creases to about 98 percent of its value when the object 
is at infinity, and the position of the object then hardly 
affects the focal properties of the system. 


V. VARIATIONS IN BORE 


There were two limitations in the experimental set- 
up which prevented a direct field plot for different bore 
diameters. The first involved the depth of the probe 
below the surface; this required that the bore diameter 
should not be reduced below 10 mm for high accuracy. 


POTENTIAL DISTRIBUTION: RELAXATION METHOD 
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The second arose from the size of the electrolytic tank 
for the preliminary tests had shown that in order to 
avoid leakage effects, the minimum ratio of outer elec. 
trode diameter to bore diameter should be at least 6/1, 
Since the maximum diameter which could be accom- 
modated was 120 mm, it was thus possible only to use 
a maximum bore diameter of 20 mm. To overcome 
these difficulties use was made of the “principle of 
dynamical similarity.” This, in effect, states that for the 
same configuration of electrodes the electron paths 
become scaled up by a factor of m when all the dimen- 
sions of the lens system are scaled up by the same 
factor. An application of this principle provided a 
means for introducing variations in bore diameter jn 
the five sets of electrodes chosen for the experiment 
whose dimensions are shown in Table I. 

The last column indicates the variation in bore when 
all sets are transformed to have a common thickness of 
1.65 mm, and the same separation of 13 mm. Figure 6 
shows that the potential at the saddle point increases 
rapidly and linearly as the bore is opened out up to 
values of about 3S; above this point it tends towards 
a constant value independent of further increase in 
aperture diameter. The potential gradient, however, 
decreases with increased aperture and the integrated 
effect over the field results in an increased focal length. 
When the bore diameter is of the order of the separa- 
tion, the slope of the focal length/bore diam curve is 
small and agrees as to the order of variation with that 
found by Bruck and Romani up to their limiting ratio 
of 1.85. Beyond this point and when it has increased 
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TaBLE I. Electrode dimensions for studying variations 
in bore diameter. 











Bore diam Thickness Separation Corresponding 
No. mm mm mm bore mm 
1 10 1.65 13 10 
2 20 1.65 13 20 
3 20 0.965 8 32.5 
4 20 0.635 5 52 
5 20 0.406 3 86.7 








to some 2.5.S, the present work shows that the focal 
length rises steeply with a slope increased by a factor 
of four. A feature of the curve is the abrupt change in 
slope at the ratio 2.5. It is interesting to note that the 
alteration of focal length with bore diameter and sepa- 
ration are in opposite directions, the former effect 
being much the large at large values of D. 


VI. THE RELAXATION METHOD 


The principles of this method and its application to 
the case of an immersion lens have been described in 
detail by Motz and Worthy’ and in another paper by 
Motz and Klanfer.* The unit of mesh length in the 
present case was taken as the thickness / of an elec- 
trode, and a coarse net of squares constructed for the 
whole field, in which the bore diameter occupied twelve 
such units to within 1 percent, while the separation was 
kept at the value 10. 

Because of the symmetry about the axis and mid- 
plane,‘a quarter section with the specified boundary 
conditions was taken and subdivided into three parts: 
(1) a fine net of mesh length 3/ near the outer electrode, 
(2) a fine mesh of length 3/ near the midplane, (3) a 
medium net of mesh length ¢ for the rest of the space. 
These are shown in Fig. 7. The line which separates the 
fine net from the medium net was considered the com- 
mon boundary for both nets (shown thick in the figure). 
The potentials at the nodal points on the common 
boundary line for the fine net region were determined 
by interpolation, keeping the values of the potentials 
on the common boundary line fixed. By suitable choice, 
points in the fine net region were included in those for 
the medium net region and their values recalculated. 
This changed the values of potential on the common 
boundary at the nodal points of the medium net. The 
potentials at all the nodal points for the fine net region 
on the boundary were redetermined by interpolation, 
and assuming these values constant, the values of 
potential at the nodal points in the fine net were re- 
calculated. This process was repeated till constant 
values were attained at the common boundary for both 
types of net; the accuracy was taken far enough to be 
correct to two places of decimals for comparison with 
the tank experiments. It was necessary, starting from 
the initial assumption of a uniform distribution between 
* 7H. Motz and W. D. Worthy, J. Inst. Elec. Engrs. (London) 
92, 522 (1945). 


1945 Motz and L. Klanfer, Proc. Phys. Soc. (London) 58, 30 


outer and middle electrodes, to carry out some fifteen 
to twenty operations before achieving the final result 
for the distribution. The ES field was constructed 
from this by drawing equipotentials at }-volt intervals 
throughout the field. 


VII. RELAXATION AND TANK: A COMPARISON 


The two methods of field plotting were compared 
by reading off the potential on the axis at 1 mm inter- 
vals beginning from the midplane. The percentage 
variation in potential is plotted against axial distance 
z in Fig. 8, where it can be seen that, as compared with 
the tank, the relaxation results for the axial distribu- 
tion vary from about —1.5 percent to +1.5 percent; 
the relaxation results are lower near the middle elec- 
trode, but higher near the outer one. This has been 
readily accounted for by the 1 percent variation in bore 
diameter. Since the distributions furnished by the tank 
method are uniformly low by 1.3 percent at all points 
because of the ratio of impedances involved in the 
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Fic. 8. Tank and relaxation: a comparison. 


circuit, the numerical method must then be almost 
correct for axial points in the vicinity of the outer 
electrodes and it follows that it must be correct to the 
same extent for all other axial points. Analysis of the 
experimental results indicates that the 3 percent varia- 
tion at the saddle point is to be expected for a 1 percent 
variation in bore. Comparison over the whole field reveals 

(a) the equipotentials cross at the saddle point at 
an angle which within the limits of experimental error 
agrees with the theoretical value of 54°44’; 

(b) the slope of the equipotentials due to the edges 
of the electrodes have an increased value in the nu- 
merical method; this could not be detected in the tank 
without special amplifying circuits; 

(c) in the region between the electrodes, there is 
almost perfect agreement in form as between the two 
distributions. 

In conclusion, the authors wish to express their 
thanks to Dr. Motz for a helpful discussion on the 
relaxation method. 
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A solution is given for the problem of the scattering of plane electromagnetic waves from a sphere with 
a concentric spherical shell. The solution is general, and under appropriate conditions is reduced to the 
well-known solution for scattering from a single sphere. 





INTRODUCTION 


HE theory of scattering of plane electromagnetic 

waves from an isotropic sphere has been worked 

out by Mie! and is concisely presented by Stratton.’ 

This paper will present the solution for scattering from a 

sphere with a concentric spherical shell. The application 

of this solution to certain problems in radar meteorology 
will be considered in a subsequent communication. 


FORMULATION OF THE PROBLEM 


Figure 1 shows the orientation of the incident plane 
wave and the scattering configuration with respect to a 
rectangular coordinate system with its origin at the 
center of the sphere. The inner sphere, the shell, and the 
surrounding medium are called regions 1, 2, and 3, 
respectively. The inner and outer radii of the spherical 
shell are a and 6. Each region is assumed to have a dif- 
ferent complex propagation constant, complex dielectric 
factor, complex characteristic velocity, and perme- 
ability. These quantities are indicated by K, & V, and 
u, respectively, with subscripts 1, 2, or 3, depending on 
the region. The incident plane wave is propagated in 
the positive z direction, and its electric vector is linearly 


Ej | 
REGION 3 







Si 

B; PoINTs ouT 

FROM PAPER 
SPHERICAL SHELL 


Y-AXIS POINTS 
OUT FROM PAPER 


Fic. 1. Plane wave incident upon a sphere with a concentric 
spherical shell. 


* This solution was worked out independently by each of the 
present authors. 

Tt Part-time research assistant, MacDonald Physics Laboratory, 
McGill University, Montreal, Canada. 

!G. Mie, Ann. Physik 25, 377 (1908). 

2 J. A. Stratton, Electromagnetic Theory (McGraw-Hill Book 
Company, Inc., New York, 1941), pp. 563-573. 


polarized in the x direction. Time dependence of the 
form e’*‘ is assumed ; but it will not be written explicitly, 
since it does not enter directly into the analysis. 

With the conditions stated above, the expressions 
for -the incident plane wave have the same formt as 
that given by Stratton.” 

It is seen that 


2n+1 mM .@ 


E,;= = Ey > (- j” “uno (1) 


n=1 n(n+ 

Eo\ = 2n+1 wy () 
B;= -(~) (— j)"—_——— Mein— m n)y 2 
V3 ms ware u— jam), @) 


where Ep is the amplitude, V3=(£3u3)~} is the complex 
characteristic velocity of region 3, and 


aft 1 COs | 
mime ’ (K3R)P, '(cos8) = hla 
a sin@ sin 
(1) dP,,\(cos@) sin . 
—s, (K3R)———-__ ¢,_ (3) 
cos 
a n(n+1) wm . sin | 
No, = ———2n (K3R)P,'(cos8) oh 
=“ cos 
1 (1) dP,,\(cos6) sin . 
+——(K3Rz, (K;R) /—————_ oly 
K3R dé cos 


1 Os 
——- a Rs, (K 3R) |'Pr (cos) dis. (4) 
K3R siné sin 


The primes at the square brackets indicate differenti- 
ation with respect to the argument K;R. Here ii, h,, 
and i; represent unit vectors in the directions R, 0, and 
¢ of a spherical coordinate system with its origin at the 
center of the sphere ; P,'(cos@) is the associated Legendre 
polynomial of the first kind, first order, and mth degree; 


(1) 
Zn (KR) is the spherical bessel function of the first 


t Note that Stratton assumes e~*“* time dependence. Therefore, 
in comparing equations written here with those of Stratton, replace 
j by —i. 
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kind, which is related to the ordinary bessel function of 
the first kind and half-integer order as follows: 


on (KR) = (r/2KR)'Jn43(KR). (5) 


The induced secondary field must now be constructed 
in three parts, one applying in each of the three regions 
defined previously. These parts are written as expan- 
sions similar to those for the incident wave, but with 
unknown amplitude coefficients. The parts applying 
outside the shell and inside the sphere will be called the 
scattered and transmitted fields, in accordance with the 
terminology used for the single-sphere problem. They 
will be indicated by using subscripts s and ¢. The formal 
expansions for these fields are identical to those for the 
single-sphere problem, although the amplitude coef- 
ficients have different values here. Thus, one can write 


2n+1 
n(n+1) 





2 . (3) . (3) 
E,= Eo > (—j)” (@,°Moint+ jon*Mein), (6) 
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J 37 n=l n(n+ 1) 


valid when R>b, and 








2 _ 2n+1 Mo . @ 
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valid when R<a. The functions 


(3) (3) 


me, and Mo,, 


are obtained by replacing 


(1) (3) 
fn (K3R) by zn (K3R) 
in Eqs. (3) and (4): 


. 
ce (KR)=(x/2KR)'H43(KR) 


is the spherical hankel function of the second kind. 
These functions are necessary in the solution for the 
scattered field, since this solution must be regular at 
infinity and must satisfy the radiation condition. Bessel 
functions of the first kind only are used in Eqs. (8) and 
(9), since the solution must be finite at the origin. In 
addition, K3 is replaced by K, in Eqs. (8) and (9). 

In the restricted region of the spherical shell, it is 
necessary to use bessel functions of both the first and 
second kinds. Here the solution may be written 


2n+1 (1) (3) 
(@n"Moin+ An°Moin 
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with the understanding that KR is the argument of 
the vector functions. 


SOLUTION FOR THE SCATTERING AMPLITUDE COEFFICIENTS 


Equations (6) through (11) represent a formal solution for the induced secondary field. All that is needed to 


complete the formal solution is the evaluation of the eight amplitude coefficients. This is done by applying the 
boundary conditions at the two surfaces of dielectric discontinuity. 


The boundary conditions at R=a are 


i, XE.=iXE,,, (12) 
(1/1) <XB.= (1/2)i1 <B.., (13) 
and at R=} are 
i XE,.=1,X (E;+E,), (14) 
(1/u2)i:X B,.= (1/us)i:< (B+B,). (15) 


These lead to two sets of simultaneous equations involving four unknowns each, as follows: 
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These two sets of simultaneous equations may be solved for the eight amplitude coefficients. However, since the 
main interest of this paper is in the scattered field, only the solution for a,* and b,° will be given explicitly here. 

In writing down the solution for the scattering amplitude coefficients, it is convenient to introduce certain 
simplifications in some of the terms of the above equations. Thus, it will be recognized immediately that (K/u)=oY, 
where w is the angular frequency, and Y=(£/)! is the complex characteristic admittance of the region under 
consideration. With this substitution, w drops out of the equations. In addition, it is convenient to let y= K;b, 
a=K;a, Ni=K,/K;3, N2=K./K; and to introduce the new functions 
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These equations may be expanded in the form 
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The above solution may be checked by letting a=v. In this case, Eqs. (28) and (29) become 
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which are identical with the Stratton-Mie solution for the scattering from a single sphere. 

It should be noted that the y-functions used here are closely related to the logarithmic derivative functions used 
by Aden® to simplify the computations in the case of scattering from a single sphere. These latter functions may 
also be used in the formulation of the present problem. However, such a procedure is more artificial here and does 
not yield the computational advantages that are manifest in the single-sphere problem. 





THE FAR ZONE FIELD AND THE SCATTERING j 
PARAMETERS Ren Vidiam ( =) Eye7iks 

Since the present problem was formulated in such a KK 
way that the equations for the scattered field are «© In+1 P,.\(cos@) dP,,\(cos6) 
formally the same as the equations in the single-sphere x = - +5,# | (37) 
problem, it is possible to utilize directly the further n=l n(n+1)t sin@ dé 
equations from that problem. Thus, the far zone scat- 
tered field in spherical component form is 
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* A. L. Aden, “Electromagnetic scattering from metal and water Pi P.1 
- cael Technical Report No. 106, Cruft Laboratory, Harvard 0 dnt [ a fg n'(cos@) 1 “ =] (38) 
7° ron : FY a : ; 1 n 1 . ’ 
(ost (1950); condensed version, J. Appl. Phys. 22, 60 a= n(n-4 1)L d6 sind 
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the total scattering cross section is 
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and the backscattering cross section is 
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The radiation pattern and the gain were measured for four series of dielectric antennas of circular cross 
section. These antennas were excited by a TE; mode in a circular metallic wave guide at a frequency of 
9275 mcps. Only one parameter was changed in each series of.rods so that the effect of this parameter on the 
pattern and gain could be studied. Measurements of gain are presented for two further series of rods in 
which the length is varied. The measurements showed that a polystyrene rod of length 6X0 and constant 
cross section has maximum gain and best pattern for a diameter of approximately 29. On the other hand 
when the taper of a Lucite rod was varied while the length and the maximum diameter were held constant 
at 60 and 0.870, respectively, the gain and the pattern behavior improved as the minimum diameter was 
decreased to zero. It is shown that losses in the dielectric produce a significant decrease in the gain. 


I. INTRODUCTION 


UMEROUS workers'~’ have investigated the di- 

electric rod antenna both experimentally and 
theoretically. The present paper is a report of an ex- 
tensive experimental investigation of four series of di- 
electric radiators of circular cross section in which the 
range of physical dimensions was larger than usual. 
In addition the radiation characteristics are correlated 
with measurements of the fields on the surface of the 
radiators. The data are arranged so as to be of value 
both for design work and for guidance and verification 
of theoretical work. 


II. PARAMETERS STUDIED AND CHARACTERISTICS 
MEASURED 


All of the measurements were made at a frequency of 
9275 mcps or a free space wavelength, Ao, of 3.20 cm. 


* The work described in this paper was done at the Defense 
Research Laboratory, under the sponsorship of the Bureau of 
Ordnance, Navy Department, Contract NOrd-9195. 

t Now at Texas Technological College, Lubbock, Texas. 

' Peter Mallach, “Dielectric directional antennas for dm and 
cm waves,” Air Materiel Command Report, F-TS-2223-RE, 
February, 1948. Translated by P. L. Harbury of Harvard Uni- 
versity, Cambridge, Massachusetts. 

2 G. Mueller and W. Tyrrell, Bell System Tech. J. 26, 837 (1947). 

* D. F. Halliday and D. G. Kiely, J. Inst. Elec. Engrs. (London) 
64, 610 (1947). 

*G. Wilkes, Proc. Inst. Radio Engrs. 36, 206 (1948). 

5 R. B. Watson and C. W. Horton, J. Appl. Phys. 19, 661 (1948). 

® R. B. Watson and C. W. Horton, J. Appl. Phys. 19, 836 (1948). 

7 Horton, Karal, and McKinney, J. Appl. Phys. 21, 1279 
(1950). 


One end of the dielectric rod radiators, which were 
constructed of Lucite or polystyrene, fitted snugly into 
the open end of a hollow metal wave guide excited in 
the TE,, mode. This TE; wave excited a hybrid mode 
(neither TM nor TE) in the dielectric rods which will 
be called the HE,, mode. The rods extended into the 
metal guide for approximately four wavelengths and 
the end of each rod was tapered gradually to a point 
in order to achieve a suitable impedance match. 

Three sets of rods were made so that each set showed 
the effect of one parameter at a time. Series A was made 
of polystyrene and showed the effect of changing the 
diameter of an untapered rod. Each rod was connected 
to the metallic guide through a short metal-covered 
transition. Series B and C were made of Lucite and 
showed the effect of changing the taper and the length, 
respectively. A D.series consisted of a single polystyrene 
rod made in the shape recommended by Mueller and 
Tyrrell.’ 

The characteristics which were measured were the 
radiation patterns (E- and H-planes), the absolute 
gain, and the electric field distribution on the surface 
of the radiating rods. A klystron type of transmitter 
and a superheterodyne receiver were used to measure 
patterns and the results were recorded continuously. 
The gain was measured by the reciprocity method.’ 


8C. G. Montgomery, Technique of Microwave Measurements 
(McGraw-Hill Book Company, Inc., New York, 1947), Vol. 11, 
Radiation Laboratory Series. 
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The surface electric fields were measured by means of 
a small probe antenna moved over the surface of the 
rods. 

The theory of the dielectric wave guide indicates that 
the same HE,, mode should exist in the dielectric 
whether it is generated by a TE,, or a TM, mode in 
the metallic wave guide that serves as a source. Con- 
siderable measurements were made on dielectric an- 
tennas excited by the 7M, mode and the results were 
essentially the same as those obtained when the rods 
were excited by a TE; mode. 


Ill. MEASUREMENTS OF RADIATION PATTERNS 
AND SURFACE FIELDS 


In an earlier paper’ it was shown that the radiation 
field of a dielectric rod can be expressed as the product 
of two terms. One of these terms depends primarily on 
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Fic. 1. Patterns in the H-plane for the A series of rods. Length 
6Ao0; material polystyrene; diameter, d, constant. 


the diameter of the rod while the other depends pri- 
marily on the length. It is evident that to a first ap- 
proximation dielectric rods of the same length should 
have the same number of lobes in the pattern. However, 
the envelope of the lobes will be a function of the di- 
ameter of the rod. This can be seen by reference to 
Figs. 1 and 2. Figure 1 shows the radiation patterns in 
the H-plane for the A series of antennas. Figure 2 
shows the envelopes of the patterns in the E-plane for 
the same series of rods. It can be seen from Fig. 1 
that the number of lobes remains constant while Fig. 2 
shows that the directivity of the envelope increases 
with the diameter. 

Figure 3 shows the radiation patterns for the B series 
of antennas. The outstanding characteristic is that the 
amplitudes of the minor lobes are attenuated more 
sharply as the amount of taper increases. Figure 4 
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Fic. 2. Envelopes of the patterns in the E-plane for the A 
series of rods. Length 69; material polystyrene; diameter, d, 
constant. 


shows a graph of the normal component of the electric 
field on the surface of the rod which was tapered to a 
point as a function of the axial length. It is seen that 
the amplitude of the field has some resemblance to a 
sine wave. The field on the corresponding rod with no 
taper is uniform. Figure 5 shows the field on another 
rod of the B series with a ratio of minimum diameter/ 
maximum diameter equal to 3. Here it is seen that the 
field intensity increases monotonically as the diameter 
becomes less except at the very extremity. In the case 
of the field shown in Fig. 4, the field strength rapidly 
decreases after the midpoint along the rod because 
most of the energy has already been lost by radiation. 
Very little of the decrease is the result of internal losses 
in the rod. 

Figure 6 shows the radiation patterns for the C series 
of antennas. As one would expect, the number of lobes 
increase with length (but not uniformly) and the width 
of the major lobe decreases. However, the minor lobes 
decrease in intensity with increasing length. This is 
not the case with rods of uniform cross section. Figure 7 
shows the normal component of the electric field on the 
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Fic. 4. The normal component of the electric field at the sur- 


face of a rod in the B series. Length 6\0; maximum diameter 
0.870; minimum diameter 0; material Lucite. 
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Fic. 5. The normal component of the electric field at the sur- 
face of a rod in the B series. Length 69; maximum diameter 
0.879; minimum diameter 0.33A9; material Lucite. 


10Ao rod as a function of distance along the rod. Again 
the surface field approximates a half-period of the 
sine wave when the rod is tapered to a point. 

Figure 8 shows the radiation pattern of the antenna 
in the D series. The plot of the electric field along the 
antenna shown in Fig. 9 does not approximate a sine 
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Fic. 6. Patterns for the C series of rods. —— E-plane. 
---- H-plane. Maximum diameter 0.879; taper 1.04 inches per 
foot; material Lucite; length L. 
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Fic. 7. The normal component of the electric field at the sur- 
face of a rod in the C series. Length 10\9; maximum diameter 
0.870; minimum diameter 0; material Lucite. 


wave as well as the fully tapered rod shown in Fig. 4. 
However, the radiation patterns were very similar. 


IV. MEASUREMENTS OF GAIN 


Figure 10(a) is a graph of absolute gain (i.e., gain 
above an isotropic radiator) as a function of diameter 
for the rods in the A series. It is seen that the gain in- 
creases uniformly with diameter to an optimum value 
at 2.0\o. This indicates that rods considerably larger 
than those normally used would be desirable both from 
the standpoint of gain and pattern. This improvement 
in gain is undoubtedly due to an increase in efficiency 
in exciting the rods. 

Figure 10(b) shows how gain is affected by taper in 
the rods in the B series. Here it is evident that opti- 
mum gain is obtained by tapering the rods to a point. 

Figure 10(c) is a graph of the gain as a function of 
length for different tapers. The solid curve is for the C 
series of rods which are made of Lucite and have a 
uniform taper independent of length. The effect of 








Fic. 8. The radiation pattern in the E-plane of the rod in the 
D series. Length 69; maximum diameter 0.87\9; minimum di- 
ameter 0.44X0; material polystyrene. The first half of the rod is 
uniformly tapered while the second half is constant in diameter. 
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Fic. 9. The normal component of the electric field at the sur- 
face of the rod in the D series. Length 6\9; maximum diameter 
0.87\0; minimum diameter 0.44d9; material polystyrene. The 
first half of the rod is uniformly tapered while the second half is 
constant in diameter. 


taper can be observed by comparing this curve with 
the dashed curve which shows the gains of a series of 
Lucite rods whose cross section has a constant diameter 
of 0.870. (This paper does not contain any other data 
for this series of rods.) The dashed and the solid curves 
indicate that the presence of the taper significantly 
increases the gain of the antenna. 

The dotted curve in Fig. 10(c) shows the gains of a 
series of polystyrene rods whose cross section has a 
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Fic. 10. Measurements of absolute gain. (a) The A series; d is 
the diameter. (b) the B series; R= minimum diameter/maximum 
diameter. (c) Gain versus length, L, for three series of rods. —— 
C series; diameter constant at 0.87\9; material Lucite. — — — 
Lucite rods with constant diameter of 0.879. - - - — polystyrene 
rods with constant diameter of 0.87)o. 


constant diameter of 0.879. Thus the dotted and the 
dashed curves apply to two sets of antennas which 
are identical except for the material. The smaller gains 
of the Lucite antennas must be attributed to the losses 
in the dielectric. 
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Fic. 5. The normal component of the electric field at the sur- 
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0.87A9; minimum diameter 0.33A9; material Lucite. 


10Xo rod as a function of distance along the rod. Again 
the surface field approximates a half-period of the 
sine wave when the rod is tapered to a point. 

Figure 8 shows the radiation pattern of the antenna 
in the D series. The plot of the electric field along the 
antenna shown in Fig. 9 does not approximate a sine 
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wave as well as the fully tapered rod shown in Fig, 4, 
However, the radiation patterns were very similar. 


IV. MEASUREMENTS OF GAIN 
Figure 10(a) is a graph of absolute gain (i.e., gain 


above an isotropic radiator) as a function of diameter 
for the rods in the A series. It is seen that the gain in- 
creases uniformly with diameter to an optimum value 
at 2.0\0. This indicates that rods considerably larger 


than those normally used would be desirable both from 
the standpoint of gain and pattern. This improvement 
in gain is undoubtedly due to an increase in efficiency 
in exciting the rods. 

Figure 10(b) shows how gain is affected by taper in 


the rods in the B series. Here it is evident that opti- 
mum gain is obtained by tapering the rods to a point. 


Figure 10(c) is a graph of the gain as a function of 


length for different tapers. The solid curve is for the C 
series of rods which are made of Lucite and have a 
uniform taper independent of length. The effect of 








Fic. 8. The radiation pattern in the E-plane of the rod in the 
D series. Length 60; maximum diameter 0.879; minimum di- 
ameter 0.449; material polystyrene. The first half of the rod is 
uniformly tapered while the second half is constant in diameter. 
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Fic. 9. The normal component of the electric field at the sur- 
face of the rod in the D series. Length 6\9; maximum diameter 
0.8740; minimum diameter 0.449; material polystyrene. The 
first half of the rod is uniformly tapered while the second half is 
constant in diameter. 


taper can be observed by comparing this curve with 
the dashed curve which shows the gains of a series of 
Lucite rods whose cross section has a constant diameter 
of 0.87Ao. (This paper does not contain any other data 
for this series of rods.) The dashed and the solid curves 
indicate that the presence of the taper significantly 
increases the gain of the antenna. 

The dotted curve in Fig. 10(c) shows the gains of a 
series of polystyrene rods whose cross section has a 
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Fic. 10. Measurements of absolute gain. (a) The A series; d is 
the diameter. (b) the B series; R= minimum diameter/maximum 
diameter. (c) Gain versus length, L, for three series of rods. —— 
C series; diameter constant at 0.879; material Lucite. — — — 
Lucite rods with constant diameter of 0.87A0. - - - — polystyrene 
rods with constant diameter of 0.87o. 


constant diameter of 0.879. Thus the dotted and the 
dashed curves apply to two sets of antennas which 
are identical except for the material. The smaller gains 
of the Lucite antennas must be attributed to the losses 
in the dielectric. 
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Design and Analysis of a Conservative Dynamic Load Simulator 
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The design for a highly versatile dynamic load simulator is described, in which the correct resultant load is 
computed electrically and applied through a hydraulic servomechanism. One of the significant features of 
this machine is that virtually all of the power absorbed from the unit to which the load is applied is conserved 
in the form of electrical energy and returned to the power lines. A complete analysis of this electrohydraulic 
system is presented for the case in which coulomb friction is negligible. Two basic applications for this 
simulator are a servo test stand designed primarily for testing autopilots and: an aircraft control load simu- 


lator for both design and test of control systems. 





I. INTRODUCTION 


OWER-ABSORBING equipment for producing 

steady-state loads in connection with endurance 
life or performance tests is commonly used by many 
industries. The greater portion of routine and experi- 
mental testing of this nature at the present time is per- 
formed on dynamometers designed to operate under 
invariant conditions and equipped to dissipate the 
output energy in the form of heat. 

There are numerous practical applications in which 
machines are required to produce loads simulating 
dynamic forces or torques. The conventional approach 
for the description of a physical problem involving 
dynamic torques has been to visualize the torques as 
individual components which may be designated in the 
following manner: 


Ig+Co+KoFD=F (1) 


where F=total torque. K¢=spring torque. C(d¢/dt) 
=C¢=viscous damping torque. I[(d°¢/df?)=Ig=in- 
ertial torque. D=coulomb friction torque. K=spring 
coefficient. C= viscous damping coefficient. = moment 
of inertia of moving parts. ¢=angular displacement. 
The sign of the coulomb friction term of (1) is deter- 
mined by the direction of the velocity ¢. 

A linear displacement system in which forces replace 
torques and distances replace angles may be analyzed 
in a form identical to (1). These linear forces or rota- 
tional torques can be obtained through the use of equip- 
ment such as springs, weights, dampers, brakes, etc. 
There are certain obvious limitations which must be 
recognized when these mechanical components are used, 
such as (a) size of the mechanical parts may be cumber- 
some, if not prohibitive ; (b) the adjustment range of the 
parameters may be too limited ; (c) the time required to 
make these adjustments may be excessive; (d) limita- 
tion on the introduction of nonlinear characteristics, 
primarily because these features may involve a great 
amount of complication. 

The particular problem of interest to the writers has 
been that of providing a highly versatile machine which 

* E. H. Gamble, Electronics Division, Carlstadt, New Jersey. 


t B. W. Hatten, Westinghouse Company, Pittsburgh, Pennsyl- 
vania. 


may be used to determine in the laboratory the oper- 
ating and dynamic characteristics of aircraft servo- 
mechanisms and controls. The aircraft controls designer 
must demonstrate that a particular servo design is ade- 
quate for actual flight controls when subjected to all 
existing inertial, hinge moment, and damping effects 
which the servo will experience in the flight operation of 
a particular aircraft. Moreover, a satisfactory test 
equipment should be useable for testing the control 
servomechanisms for a variety of aircraft having a wide 
range of power and response requirements. 

The power requirements for the control servos used 
on large subsonic and transonic aircraft are quite large, 
A machine capable of testing and loading these servos 
must be capable of dissipating large amounts of energy. 
If this energy is converted to heat, the cooling of sucha 
machine is an expensive proposition. 

To satisfy these needs, the writers have devised a load 
simulator believed by them to be original and possessing 
noteworthy and novel features in which the correct 
resultant load is computed electrically and applied 
through a hydraulic servomechanism. One of the signifi- 
cant features of this machine is that virtually all of the 
power absorbed from the unit to which the load is 
applied is conserved in the form of electrical energy and 
returned to the power lines. The need for any heat 
dissipation apparatus is eliminated by the use of variable 
displacement hydraulic pump which efficiently adds or 
removes energy from the loading shaft of the equipment 
to be tested. 


II. DESCRIPTION OF LOAD SIMULATOR SYSTEM 


The particular item being tested may be an actuator 
or some similar piece of equipment, and will be desig- 
nated in this discussion as the test item. The test item 
may assume a variety of forms, which are operated 
electrically, hydraulically, pneumatically, or in some 
mechanical manner. For the purposes of this paper the 
assumption will be made that the ultimate energy sup- 
ply for the test item is electrical and that the output of 
the test item is a torque on a shaft which rotates through 
an angle ¢. This shaft is designated as the loading shaft. 

The basic operation of the simulator is indicated in 
Fig. 1. This figure describes the flow of mechanical 
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DESIGN 


energy to and from the test item and the conversion and 
conservation of this energy in the form of electrical cur- 
rent which is returned to the power lines. 

The simulated load which opposes the test item out- 
put torque is applied to the loading shaft by means of a 
hydraulic motor. A gear train having a speed reduction 
ratio R: 1 may be used between the motor and the shaft, 
in which case the motor displacement is R¢. Suitable 
pickoffs are connected to the loading shaft to provide 
electrical signals proportional to ¢ and its first two time 
derivatives. These pickoff signals are multiplied by the 
appropriate physical constants to obtain signals pro- 
portional to the spring, friction, and inertia torques. The 
latter signals are added in a summing type amplifier, the 
output of which represents the desired load torque. 

The physical load torque being applied by the 
hydraulic motor is detected by a pickoff which transmits 
an electrical signal proportional to this torque. This 
pickoff signal voltage is compared with the output of the 
Electronic Computer in an error measuring circuit. The 
error signal V, is transmitted to the servo controller 
which consists of an internal servo loop used to position 
the stroke mechanisms of the variable-displacement 
hydraulic pump. 

The hydraulic pump is mechanically driven at a 
nearly constant speed by an induction motor. The pump 
controls the flow of oil through the hydraulic motor 
which, in turn, produces the simulated load. The func- 
tion of the pump is unique in one particular, that it 
efficiently adds or removes energy from the loading shaft 
and eliminates the need for any heat dissipation appa- 
ratus. This function is accomplished by driving the 
pump with an induction-type motor which can also 
function as a generator. Hence, energy from the loading 
shaft will be transmitted back to the electric power lines. 
One can see from an inspection of the characteristic 
curve given in Fig. 2 that such a motor has a rather 
steep torque-speed curve. For this reason, the motor will 
function in the desired manner with a relatively small 
change in speed. 

For example, when the pump yoke is in its neutral 
position, the hydraulic load on the induction motor is 
zero and the motor rotates at nearly synchronous speed. 
When the load servo demands that a torque be applied 
to the loading shaft in a direction to aid the shaft’s 
motion (as in simulating a restoring torque from a 
spring), the pump functions in the conventional manner 
—power being transferred from the induction motor to 
the loading shaft. However, when a torque is required to 
oppose the motion of the loading shaft (as in simulating 
viscous friction), the pump stroke must be reduced 
sufficiently to cause the volumetric displacement of the 
pump at synchronous speed to be inadequate to handle 
the flow of oil from the hydraulic motor which, at this 
time, is being driven as a pump. As a consequence of this 
condition, the induction motor tends to overspeed, be- 
comes a generator and transfers energy from the loading 
shaft back into the electric power lines which were the 
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Fic. 1. Operational diagram of load simulator showing energy flow. 


original source of the energy. In addition, some energy 
will be stored mechanically by the flywheel effect of the 
motor rotor. However, there will be small energy losses 
in each hydraulic and electrical component. The greater 
portion of the energy, in any case, will be conserved 
electrically. 


Ill. ANALYSIS FOR THE CASE D=0 
NO COULOMB FRICTION 


A. Over-All Electrohydraulic Servo System 


The electrohydraulic system indicated in Fig. 1 is a 
servomechanism and must be investigated from the 
standpoint of response time and stability. In order that 
satisfactory design be adopted for the dynamic load 
simulator, the system must be suitably stable and 
possess the desired frequency response. Therefore, the 
proper components must be chosen and connected to- 
gether in such a way that the desired objectives will be 
accomplished. This section of the paper will describe the 
system in more detail and analyze in particular the case 
where coulomb friction is negligible. A later section will 
deal with the case where coulomb friction is important. 
Of particular interest to the writers is the design of a 
load simulator to operate over the frequency range from 
0 to 16 cps. 

The arrangement of the necessary components in the 
electrohydraulic system is shown in the block diagram 
of Fig. 3 which shows the flow of signal or information 
throughout the system. The test item introduces in re- 
sponse to an applied signal a displacement ¢ at*the 
loading shaft; with the proper adjustment of the load 
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Fic. 2. Induction motor torque as a function of slip. 
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Fic. 3. Block diagram of load simulator servo system for D=0 showing information flow. 


parameters the electronic computer calculates the 
torque which the torque-feedback controlled load servo 
should produce to simulate the desired operating loads 
on the test item. . 

The electronic computer contains elements corre- 
sponding to each function which the load servo must 
represent and is essentially a mechanical to electrical 
transducer. A constant or time-programmed signal 
voltage is applied across a potentiometer whose brush is 
positioned by the loading shaft. This brush voltage is a 
function of the shaft displacement ¢. For a linear po- 
tentiometer this voltage is proportional to ¢. The type 
of potentiometer recommended by the writers and 
employed in the experimental prototypes of this ma- 
chine is described by G. A. Korn.' A signal proportional 
to shaft velocity ¢ is produced by an electrical motor- 
generator whose rotor is turned by the loading shaft. 
The two-phase motors which are widely employed in 
analog computer servos may be used in this application. 
One winding of the stator is connected across a fixed ac 
voltage source. The voltage induced by rotor motion 
into the other winding will be proportional to ¢. A small 
dc generator or tachometer may be used. In this case, 
the rotor current or voltage is proportional to the rotor 
velocity. A detailed knowledge of the generator charac- 
teristics found from experimental data enables one to 
find the proportionability constants or parameters. The 
shaft acceleration ¢ may be obtained by using a strain- 
gauge accelerometer such as that used in the flight test 
instrumentation of aircraft or by electrically differ- 
entiating by means of a differentiating dc amplifier the 
output voltage of the tachometer proportional to ¢. 

These signals are caused to represent the proper 
amount of spring, viscous damping, and inertia torque, 
respectively, by varying the signals in direct ratio to the 


1G. A. Korn, Rev. Sci. Instr. 21, 77-81 (1950). 


physical constants K, C, and J, either by using electronic 
multipliers or by changing the voltage level of each 
element. These voltage levels may be programmed with 
time, thereby simulating the varying load conditions 
existing over an entire flight. Any of these pickoffs could 
be employed separately or in combination with the other 
pickoffs, according to the types of load to be simulated. 
For example, one may wish to operate a servo into a 
varying but purely spring load. In this case J and D are 
zero, and K only has value. The potentiometer pickoff 
for the angular displacement only will then be required. 

The actual torque F'4 on the loading shaft is measured 
by means of a torque pickoff which produces a pro- 
portional electrical signal V4. The error V, between 
these two voltages is determined in an appropriate cir- 
cuit and transmitted to the servo controller. It will be 
noted that an input signal V; can be superimposed on 
the error V, to produce a change in the neutral position 
of the effective spring. This is equivalent to adding a 
fixed voltage to K,K¢. 

The servo controller constitutes an internal loop in 
the load simulator servo system. This servo controller 
transforms the error V, into an angular displacement of 
the yoke of the variable pump. A key factor in the 
operation of the system will be the physical reproduction 
of a satisfactory transfer function for this internal closed 
loop. The error V, between the voltage V, and the 
feedback signal V, from the pump yoke is amplified and 
transmitted to the valve actuator which in turn pro 
duces a displacement X, of the transfer valve. Conse- 
quently, the valve supplies a flow of oil Q, to the control 
cylinder which displaces the yoke of the variable 
displacement pump through an angle @,. The pump 
controls the flow through the hydraulic motor, thus 
creating the pressure on the loading shaft required to 
produce the load torque F4. The transfer function o 
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DESIGN AND ANALYSIS OF 
each component will be determined separately before 


considering the complete loop. 


B. Component Transfer Functions 


The transfer function of a frequency-sensitive com- 
ponent is commonly indicated by the symbol KG(s), 
where K represents the static gain and G(s) contains all 
frequency-variant factors. The symbol s represents the 
laPlace complex frequency variable. At zero frequency, 
the quantity G(s) becomes unity by definition. Hence, in 
certain components having relatively high frequency 
response characteristics, the value of G(s) remains 
virtually equal to unity at all frequencies in the servo 
range. 

In this analysis, the respective transfer functions will 
be designated by the following subscripts: @ indicates 
power amplifier, b indicates pump stroke pickoff, c indi- 
cates pump control cylinder, d indicates transfer valve 
actuator, f indicates torque pickoff, m indicates hy- 
draulic motor, o indicates summing amplifier, p indicates 
variable pump, and 2 indicates transfer valve. Further- 
more, the factors in each transfer function will be ex- 
pressed in terms of the following customary symbols: 
w,=natural frequency (rad/sec), ¢=viscous damping 
ratio (dimensionless), and T= time constant (sec). 

The following symbols with the correct component 
subscript describe the important physical quantities: 
V=volume of oil under compression (in.*), B=bulk 
modulus of oil (psi), P=cylinder differential pressure 
(psi), A=piston area (in.*), L=leakage coefficient 
(in.*/sec-psi), J = moment of inertia about main shaft of 
rotation (in.-lb-sec?), r=lever arm (in.), f=coefficient 
of viscous friction about main shaft of rotation (in.-lb- 
sec), d=volumetric displacement of hydraulic unit 
(in.3/revolution) and R=speed reduction gear ratio. 


1. Computer 


The outputs of the three test item pickoffs are multi- 
plied by the respective coefficients (spring, damping, 
and inertia) and the total of these inputs to the summing 
amplifier is multiplied by its transfer function to produce 
the voltage Vp which represents the desired torque. 


Vp=K .G)(s)[KiK+KCs+K3ls*]o. (volts) (2) 
2. Error Detector 


The error voltage V, is the algebraic sum of the three 
inputs to the error detector. It is apparent that if the 
displacement @ is positive, a positive spring torque 
ordinarily should exist. Thus, the voltages Vp and V4 
would both be positive, and the latter voltage should be 
subtracted from Vp to reduce the error voltage V, to 
zero (i.e., zero if V;=0). A similar argument applies to 
the polarity of the derivatives of ¢. The voltage V;, 
being an external input, appears directly in the error. 
The equation of the circuit may be written as 

V.= Vit Vp- Va. 


(volts) (3) 


A 
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3. Amplifiers and Pickoffs 


The respective transfer functions of all amplifiers and 
pickoffs depend upon the gain and compensation re- 
quired by the closed loop and will be designated as 
KG(s), employing the correct subscript. 


4. Valve Actuator 


The previous experience of the writers? with several 
different designs of such actuators indicates that a 
typical transfer function will be of the following form: 








Ka 
KiGa(s)= 
[(s?/wna*) + (26a/wna)s+1](Tas+1) 
Ka 
= —. (in./volt) (4) 
ga(T as+1) ; 


5. Transfer Valve 


For practical purposes, the transfer function of the 
valve may be considered as producing a flow propor- 
tional to the displacement. Thus 


Q,=K,X,_ (in./sec), 


(5) 
where K, is the valve gain in in.*/sec per unit displace- 
ment of the valve. Actually, the value of K, decreases as 
the load pressure becomes larger, and there is usually a 
small dead spot in the middle of the valve travel. How- 
ever, the effect of both these characteristics can be made 
small, and they will be neglected in order to linearize the 
equations. This step permits the system to be analyzed 
in accordance with the procedure used by Brown and 


Campbell.* 
6. Control Cylinder 


The oil flow through the control cylinder which 
actuates the pump yoke differs from that through the 
valve only by the amount of loss resulting from leakage 
Q. and compressibility Qc. 

0.=Q.—-Q1—-Qe (in.3/sec). (6) 


The leakage flow Q7 of the cylinder is considered as 
going through an external path and is assumed to vary 
directly with the cylinder differential pressure P,. 


Q.=L.P,. (in.*/sec). (7) 


The oil compressibility ‘flow’ Qc is equivalent to the 
rate of volumetric change in the oil under pressure and 
is defined as follows: 


Qc=(V./B)sP. (8) 
The oil flow Q, through the cylinder equals the piston 


(in.*/sec). 


2 “Rotary hydraulic transfer valve and actuator,” Report No. 
R 49-13 to U.S.A.F., Curtiss-Wright Corporation, Airplane 
Division, Columbus, Ohio, May, 1949. 


3 Brown and Campbell, Principles of Servomechanisms (John 
Wiley and Sons, Inc., New York, 1948). 
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velocity times its area. Thus 


Qe=A-sy (in.*/sec), (9) 


where the piston displacement y is in inches. If the 
piston actuates the pump yoke through an arm length r, 
inches, one may let y=r.6,, where 0, is the angular dis- 
placement of the yoke in radians, and write 


0.=17-A SO 


provided that @, is small. 

Furthermore, the torque on the yoke is expressed as 
P.Ag#f.. The experience of the writers and other in- 
vestigators indicates that the inertia and friction forces 
on a pump yoke are actually small at frequencies as low 
as 16 cps. They have found that the principal forces 
opposing the displacement are hydraulic and result from 
various design conditions such as the porting arrange- 
ment. However, for the purposes of this analysis, the 
motion of the yoke is described by 


P.Age=(J-s°+f-s)0, (in.-lb), (11) 


where J, and f, are about the yoke pivot axis. Substi- 
tuting (7), (8), and (10) into (6), 


r.A $0,=Q0,—(V./B)sP.. (12) 


Solving (11) for P., substituting in (12), and solving for 
the transfer function, one finds that 


(in.*/sec), (10) 





Oy Ae 
KG(s) = . 
Q, (reAc)*s+(Le+V-s/B)(J-s?+ fs) 


K. K. 
i. SL (s?/wne?)+ (2¢./wne)S+ 1] s Sde 


radians 
( ), (14) 
in.*/sec 


} 
(r.2A +f) 








where 
reAe B 
— 
r7Ae+feLe JV. 
J L.B+]feVe 
LIV BreAet+fel.)} 
7. Variable Pump 


The delivery of the pump varies directly with the 
pump speed and the sine of the yoke displacement angle. 
If one assumes that the torque-speed characteristic of 
the induction motor driving the pump is sufficiently 
steep to maintain a substantially constant speed, and 
that the yoke angle is small, the pump output flow is 


Q,=K,0,. (in.3/sec). (15) 








8. Hydraulic Motor and Gear Train 


The determination of the transfer function of this 
combined component is similar to that for the control 


H. GAMBLE AND B. W. 





HATTEN 


cylinder. In a hydraulic type motor with a speed reduc- 
tion gear ratio equal to R:1, the relationship between 
output torque and pressure is expressed as follows: 


Torque= Pm(dm/27)R=TmXR. (16) 


This relationship exists since mechanical work must 
equal hydraulic work. For N revolutions of the motor 


Mechanical work = Force X 227rn,N = TmX2aN 
Hydraulic work= P,,X Volume= P»,Xd,N. 


Therefore, 
T a= Pubbn/ 22. 


In a comparison of the torque given in (16) with that 
used in (11), it is apparent that when a cylinder working 
at a radius 7, is used to apply torque to the loading 
shaft, the torque is expressed as P»Amfm. Since the 
torque applied by the geared hydraulic motor is given as 
PndmR/27, it follows that the transfer function of the 
motor will be identical in form to that for the control 
cylinder except that the factor Amfm is replaced by 
d,R/2x. The counterparts of (11) and (14) are 


P(dmR/2r) = (Ims?+ fms)O=F 4 (in.-lb), (17) 


where J,, and f, apply to all parts mechanically con- 
nected to the hydraulic motor, but referred to the 
loading shaft, and 


o Km 
K,G»(s)=—= 
Oy SL (S?/wnm?)+(25m/Wnm)s+1] 


radians 
eae — (: ) (19) 
gm \in.3/sec 

















where 
dmR 
~ Mnf (dmnR/20)?+ finLm] 
B 3 
Onn = V [(dnR/2n)*+ fan] 
I mLmB+ fmV m 
on= 


2{ImV mBL(dmR/2)?+ fnkm]}* 
Equations (17) and (18) apply only to the condition that 
there is no input torque from the test item. 

C. Servo Controller Loop Transfer Function 


The open-loop transfer function of this internal loop 
is 








Ve K.KaK .K-KiG.G» 
—=K,GaKGiK.K GK Go= . (19) 
V. s(Tast+1)qage 

The closed-loop transfer function is 
a) K,KaK.KGa 
—- : = KG;,(s). (20) 


V. s(Tast+1)qage+KaK aK K-KiG.G> 
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It will be evident later that this is a key function in the 
operation of the system. A high speed of response is re- 
quired of this loop. 


D. System Transfer Function 


For the purpose of determining stability, the servo 
will be visualized as receiving an input signal V; at a 
time when there is no torque produced by the test item. 
In this situation, the hydraulic motor and gear train 
must work against the inertia and friction inherent in 
the loading shaft and the test item. The pressure de- 
veloped to do this work constitutes one feedback; the 
displacement ¢ the other. Later, the closed-loop transfer 
function will be determined on the basis of the torque 
F, resulting from an input ¢ from the test item. In this 
case, it is the test item which must overcome the loading 
shaft inertia and friction (in addition to opposing the 
torque from the hydraulic motor and gear train). 

Hence, the degree of stability is determinable from 
the following open-loop transfer function: 


ae K KG] (—) -) 
HEME} & 


Utilizing (2), (17), (18), and (20), the following relation- 
ship is obtained : 
VotVa 

Ve 








K,K.KaK.KGo 
s(Tas+1)qage+KaKaKvK-KiGiGs 


where 





fron (22) 


Km 
v(s)\=—[_K.Go(K K+ K.Cs+ K3Is?) 
Sm 
+ KG (J ms?+ fms) ]. 


By replacing s with jw, the frequency response charac- 
teristics may be obtained. Then the degree of stability 
and the speed of response may be determined by 
plotting these data in the form of a Nyquist diagram. 
Note that it will be pointed out later in this section that 
the value of C used in (22) must be modified to secure 
ideal performance and taken into account in the 
stability analysis. 

The closed-loop transfer function with ¢ as the input 
involves a solution for F4 in terms of Q, and ¢. In this 
instance, the hydraulic motor does not work against the 
loading shaft inherent inertia and damping, and one 
may write the counterpart of (12) as 


d»AR Fu 
—so=Q,- color e (23) 


2r 
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Since P,=22F 4/d,R, one finds that 
d,»R\? dmR Vin 
(—) so=—_0O>- (1.+—s) Pa (24) 
2r 2x B 


Also, for V;=0, (3) becomes 


V.=Vo—Va. (25) 
Then, using (2), V, is given by 
V.=KoGo(KiK+K2Cs+Kzgls*)@—K G;F 4. (26) 


Since 0,= K,K.G;(s)V., one can now write (24) as 





d,»R\? d,R 
( ) so=——K ,K Gi_KoGo(Ki,K+ K»Cs 
2a T 


Vim 
+ K3Is*)¢—KG;F 4 |— (1n+—s) a (27) 


Finally, the transfer function is given by 


(dm R/ 2m) K pK: KoGGo(KiK+ K-Cs+ K3Is*) 
Fa — (d,,R/21)?s 


—o . (28 
? (dmR/2m)K pK iK sG.Gst+ Ln+(Vm/B)s ™ 





Since in the ideal situation (28) should be a pure 
quadratic expression, we can make the following obser- 
vations concerning desired conditions in the operating 
range: 


(a) The products obtained by multiplying the controller trans- 
fer function K;,G;(s) by both Go and Gy should represent sub- 
stantially constant amplitude responses with negligible phase shift. 

(b) The time constant 2rV,,/B(d»RK,K;iK sG;Gs+2rLm) 
should be quite small. 


(c) The leakage coefficient L,, should be small relative to 
K,KiKsG,G;d,,R/2r. 

(d) The coefficient K2C should be supplemented by a quantity 
KC’ equal to d»R/2xK,K;KoG;Gp in order to eliminate the nega- 
tive term. This requires that the setting of the damping multiplier 
be changed from C to C+C’. 


If conditions (b), (c), and (d) are fulfilled, the first con- 
dition becomes much less critical. When all the above 
conditions are met, (28) reduces to the ideal form: 


F4 K,Go(s) 
—= (K,K+ K.Cs+ KJ’), 
@ KjG;(s) 


which is identical to (26), for V. equal to zero. 
The actual response will, of course, differ from the 


ideal in that it will be frequency sensitive to some 
extent. 





(29) 


IV. ANALYSIS FOR CASE D#¥0, COULOMB 
FRICTION PRESENT 


The effect of coulomb friction on the performance of 
three basic types of speed-regulating and position- 
regulating servomechanism systems was studied by 
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Fic. 4. Control load simulator. 


McCann, Lindvall, and Wilts making use of an electric 
analog computer. Coulomb friction has been defined as 
a constant force or torque, independent of speed, and 
always opposing motion. The discontinuous nature of 
coulomb friction considerably complicates the mathe- 
matical analysis of the performance of a system in which 
it is present. 

Den Hartog’ has analyzed a simple mechanical sys- 
tem with one degree of freedom when coulomb friction is 
present. McCann ef al.‘ have shown that this system 
with high friction ratios can grab, chatter, or be held 
still for appreciable periods of time. Moreover, coulomb 
friction can easily be simulated by electric analogs 
and introduced into a servo system. An exact solution of 
(1) requires a step-by-step integration of the equation of 
motion for the particular case under study which is 
beyond the scope of the present paper. 


V. APPLICATION OF BASIC DESIGN 


After a detailed study of the analysis and basic design 
of such a system, one naturally questions the feasibility 
of building a practical design for such a machine. To 
meet rigid requirements of frequency range and mini- 
mum angular phase error, some detailed development is 
required. A study of the component requirements and 
available components has been made. For operation up 
to 16 cps and large power requirements, some develop- 
ment of current model variable displacement hydraulic 
pumps and the hydraulic motor will be needed. These 
will be available in the near future. The other summing 
amplifiers, transfer valves, and pickoffs have already 
been developed to a suitable degree of accuracy during 
the postwar era by Navy and Air Force contractors. 


* McCann, Lindvall, and Wilts, “The effect of coulomb friction 
on the performance of servomechanisms,” AIEE Winter Conven- 
tion Paper, 1948. Technical Paper 48-84. 

5 Den Hartog, Mechanical Vibrations (McGraw-Hill Book Com- 
pany, Inc., New York, 1940). 
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Prototype models of two simulators are under experi- 
mental development. One is a servo test stand for the 
testing of autopilot and remote pilot servomechanisms 
and the other is for use with an electronic flight com- 
puter to simulate the characteristics of the control 
systems of a particular aircraft for training purposes. A 
brief description of each of these is given below. These 
are meant to be examples only of the use of such a 
versatile machine whose basic design is given in this 
paper. The detail design and experimental results for 
these examples will form the subject matter of a later 
paper. 


A. Control Load Simulator Used as a 
Servo Test Stand 


One significant application of the system described in 
Sec. II of this paper is that of a servo test stand. This 
apparatus will test the dynamic characteristics of 
numerous servomechanisms and is particularly suited 
for testing autopilot servos. Aerodynamic flight charac- 
teristics for a particular type of aircraft may be intro- 
duced into the electronic computer of the system, such 
that the proper reaction loads upon each autopilot servo 
will be applied as a function of flight time. This feature 
permits the dynamic test of an autopilot for both typical 
and extreme flights of the particular aircraft in which 
the autopilot is to function. Figure 4 depicts a three 
channel load simulator for testing autopilots. 


B. Aircraft Control Load Simulator 


A second application of the basic load simulator sys- 
tem is that of a control load simulator for use in the 
design or dynamic test of the control system of a 
particular aircraft. When used in the conjunction with 
an electronic flight computer which utilizes the aero- 
dynamic characteristics of the airplane and calculates 
the dynamic flight characteristics, the voltage signals 
proportional to attitude, altitude, and airspeed from the 
flight computer are transmitted to the computer section 
of the load simulator. These signals and breadboard or 
prototype model of the aircraft control system enable a 
controls engineer to determine the proper spring, 
damping, inertia, and coulomb friction characteristics of 
the control system which will provide the optimum 
dynamic loading and response. The load simulator sys- 
tem described in Sec. II is particularly adapted for such 
use, since most aircraft control systems are hydraulic in 
nature and its electrohydraulic system possesses high 
sensitivity at the null or zero position of the controls as 
do the aircraft hydraulic systems. The equipment 
shown in Fig. 4 also may be utilized to provide these 
outputs for the control loading of the rudder, elevator, 
and aileron of an aircraft. 


VI. CONCLUSIONS 


The basic design of a dynamic load simulator has been 
presented which conserves the energy absorbed from the 
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item under test and returns this energy to the electric 

wer lines. The analysis of this electrohydraulic system 
without coulomb friction present is complete. The sys- 
tem described will provide a near perfect loading of the 
output shaft of the item under load. An analysis which 
includes coulomb friction must follow a step-by-step 
method of integration of the equation of motion for each 
particular form of forcing function. The application of 
this design to a control load simulator is of particular 
significance, since it enables the aircraft designer to 
flight test his aircraft’s control system on the ground 
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where adjustments can be made to secure an optimum 
design. The electrohydraulic system described should 
find numerous other applications. 
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Internal friction studies have been made of an alloy 88 percent Mn—12 percent Cu. This alloy when 
annealed at 925°C and quenched to room temperature has a tetragonal structure of axial ratio 0.97. Twins are 
readily formed along the (101) and (011) planes. It has been suggested that relaxation across these twin 
boundaries would account for certain observed internal friction. Verification of this is shown. 


INTRODUCTION 


HE alloys of copper and manganese have been 
found to exhibit many unusual and interesting 
structural and physical characteristics. One of these, the 
high internal friction at certain compositions, has been 
the subject of interest to the authors of this paper. 
From Fig. 1,! the equilibrium diagram for the copper- 
manganese system, it is seen that copper and manganese 
form a continuous series of solid solutions at the high 
temperatures. In this y-region the structure is face- 
centered cubic at the copper-rich end and face-centered 
tetragonal at the manganese-rich end. The transition in 
structure occurs, in an apparently continuous manner, 
at 82 percent manganese as shown in Fig. 2.” These 
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Fic. 1. Equilibrium diagram for the system copper-manganese 
(after Dean). 





‘Dean, Long, Graham, Potter, and Hayes, Trans. Am. Soc. 
Metals 34, 458 (1945). 


* Dean et al., Trans. Am. Soc. Metals 34, 455 (1945). 


diagrams were drawn from structural studies on speci- 
mens of the alloy which were quenched from the high 
temperature 7-solid solution region to room tempera- 
ture. The structures obtained were assumed to be 
representative of the high temperature alloy. Investiga- 
tion of a sample of composition 88 percent Mn and 12 
percent Cu which had been quenched to room tempera- 
ture after annealing at 925°C has disclosed unusually 
high values of internal friction under certain conditions.* 
Zener* has suggested that the tetragonal alloy was 
caused by twinning along (101) and (011) planes, and 
that the internal friction observed could be attributed to 
the stress relaxation across the twin boundaries. Struc- 
tural studies of this alloy® showed that these twins did 
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Fic. 2. Lattice parameter as a function of composition of the 
-phase (after Dean). 


3 F. T. Worrell, J. Appl. Phys. 19, 929-933 (1948). 

*C. Zener, Elasticity and Anelasticity of Metals (University of 
Chicago Press, Chicago, 1948), p. 160. 

5 F. T. Worrell, J. Appl. Phys. 19, 930-931 (1948). 
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Fic. 3. Block diagram of apparatus used in internal 
friction measurements. 


indeed exist, but no conclusive evidence of their con- 
nection with the internal friction behavior could be 
found at the time. It has been the purpose of this in- 
vestigation to see if the internal friction could be 
attributed to this stress relaxation across the twin 
boundaries. 


MEASUREMENT OF INTERNAL FRICTION 


Measurements were made on a bar of this alloy } inch 
in diameter, and about 7 inches long using thesetup illus- 
trated in Fig. 3. The bar was supported by Nylon threads 
at the nodes for transverse vibrations. The positions of 
these nodes were determined from calculations, and the 
results checked by experiment. This experiment con- 
sisted of allowing the point of a pencil to rest lightly on 
the bar above the support. If no effect on the amplitude 
of vibration was observed, the support was considered 
to be located at the node. 

The transverse vibrations of the bar were excited by 
the action of an electromagnetic drive on a small piece 
of soft iron cemented to the end of the bar. This 
electromagnet was driven by a Hewlett-Packard oscil- 
lator which had a resistance network inserted in its 
R-C circuit to give the small variations in frequency 
that were necessary for this experiment. The vibrations 
of the bar were detected by an identical electromagnetic 
pickup at the other end, and transmitted to the vertical 
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Fic. 4. Effect of annealing on internal friction in 88 percent 
Mn-—12 percent Cu. 
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plates of a cathode-ray oscilloscope for observation of 
the amplitude. 

The bar and equipment shown in Fig. 3 were placed 
in an insulated box where temperatures could be varied 
from about —40°C to room temperature. The specimen 
was cooled by a blower which circulated the air under a 
copper tray, placed in the top of the box, containing a 
mixture of dry ice and acetone. By adding small 
quantities of dry ice, the temperature of the box could be 
maintained at a specified temperature for long periods of 
time. 

The measurements taken were the resonant fre. 
quency, vo, of the bar (about 700 cps in this case) and 
the frequency difference, Av, between points of half. 
maximum amplitude, as was indicated on the oscillo- 
scope. The internal friction was calculated from the 
equation® 


Q- = Ap ‘V3 vo. 


EFFECTS OF ANNEALING ON INTERNAL FRICTION 


The internal friction measurements were made on a 
specimen first after quenching from 925°C and then 
after successive 72-hour annealing periods at 625°C. 
Dean and his associates’ have stated that an anneal of 
24 hours at 625°C (or 72 hours at 525°C) is sufficient to 
give the room temperature equilibrium structure. Since 
the manganese reacts rapidly with the atmosphere at 
the elevated temperatures, and since it sublimes readily, 
the annealing was performed in an atmosphere of argon. 

In Fig. 4 is plotted the internal friction as a function 
of temperature. The curves correspond to the different 
heat treatments as indicated and show how the internal 
friction decreased as the annealing time was increased. 
The largest change in the internal friction is observed 
between the specimen annealed for 72 hours and that 
for 144 hours. A further anneal of 72 hours accomplished 
a considerably smaller change in internal friction than 
was produced during the second 72-hour anneal. 


METALLOGRAPH AND X-RAY STUDIES 


Since the greatest decrease in internal friction for the 
specimen was observed between the 72-hour and 144+ 
hour anneals, these specimens were investigated by 
metallographic and x-ray diffraction methods. 

Studies were conducted with a Leitz Metallograph at 
200X magnification with green light illumination ona 
surface that had been previously prepared by me- 
chanical polishing down to 4/0 paper and by electro- 
lytical polishing.* In addition, x-ray studies using 
chromium K, radiation were made of powder samples 
of the various specimens mounted in a 71.62-mm Picker 
camera. 

Specimens of this alloy quenched to room temper 
ture from 925°C showed a fine twin structure which has 


6 Reference 4, pp. 60-62. 
7 See reference 1, pp. 474-475. 
8 See reference 3, p. 930. 





_—s 





ipera- 
*h has 





THE ROLE OF TETRAGONAL TWINS IN 











(b) 

















Fic. 5. Typical surface appearance of 88 percent Mn— 12 percent 
Cu. (a) Quenched from 925°C (166*X). (b) Annealed 72 hours 
a (c) Quenched from 925°C. Fine twins within twins. 
(300X). 
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been previously reported.® A typical specimen is shown 
in Fig. 5(a).{ Metallographic examination of the surface 
of the specimen annealed for 72 hours revealed highly 
polished regions with the same fine twin structure 
(Fig. 5(b)). The remainder of the surface was obscured 
by a scale. It was suspected that this scale arose from no 
defect in the polishing procedure, but was evidence of 
chemical action of the electrolyte on a second phase in 
the alloy. The existence of a second phase was confirmed 
by the x-ray diffraction photograph of the specimen, 
which showed that the lines associated with the 
tetragonal structure of the quenched specimen were 
reduced in intensity, and new lines had appeared. 
Metallographic and x-ray examination of the specimen 
annealed for 144 hours showed almost complete dis- 
appearance of the tetragonal phase. The structure of the 
second phase was not investigated in detail beyond 
ascertaining that it was not tetragonal. 


RESULTS AND CONCLUSIONS 


Since successive annealing of the quenched specimen 
of the Cu-Mn alloy decreased the internal friction of the 
specimen, the source must be in the structure of the 
quenched specimen. Investigation by metallograph has 
shown that visible regions of the twinned structure de- 
creased in area as the annealing progressed. X-ray dif- 
fraction studies have shown that the amount of the tet- 
ragonal structure of the y-phase alloy is also reduced. 
Since this twinning is associated with the tetragonal 
structure, it can be said that the high value of internal 
friction of the Cu-Mn alloy, 88 percent Mn and 12 per- 
cent Cu, is due to the stress relaxation associated with 
the twins in the y-phase. 

The work herein reported is an extension of work 
submitted by one of the authors (A.V.S.) as a thesis in 
partial fulfillment of the requirements for the degree of 
Master of Science in Physics at Rensselaer Polytechnic 
Institute. The authors wish to express their indebted- 
ness to Professor A. A. Burr of the Department of 
Metallurgy for his help on some aspects of this problem, 
and for permission to use the facilities of the de- 
partment. 

We wish to express appreciation for the support given 
by a grant-in-aid from the Research Committee of 
R.P.L., and for the gift of high voltage equipment by the 
Behr-Manning Corporation. 

Figures 1 and 2 are reproduced by permission of the 
American Society for Metals. 


9 See reference 3, p. 930. 

t Since the work reported in reference 9, Bowles, Barrett, and 
Guttman have reported in Trans. Am. Inst. Mining Met. Engrs. 
188, 1478 (1950) on the similar indium-thallium system. They 
have found not only the {110} twins reported in 9, but fine twins 
within these. It is of interest to note that these were also found 
in one case with the copper-manganese alloy, but were not re- 
ported. The picture showing these is reproduced in 5(c). 
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Diffusion in Grain Boundaries and Their Structure 
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In order to gain an insight into the structure and properties of grain boundaries, diffusion along grain 
boundaries has been studied as a function of relative orientation of the grains and of the boundary. Experi- 
ments were made by comparing the depth of penetration of silver along the boundaries of columnar copper 
in the columnar direction [100]. The results indicate that, for angles of misfit between the grains, greater 
than 20° and smaller than 70°, diffusion along grain boundaries is greater than volume diffusion (near 700°C) 
and reaches a maximum around 45°. Outside of that range of relative orientation the rate of grain boundary 
diffusion is not greater than volume diffusion. These results are correlated with the information available 
about the energy of grain boundaries, and it is concluded that at small angles the dislocation model of grain 
boundaries describe satisfactorily the situation, while at larger angles the grain boundaries are better de- 


scribed in terms of areas of fit and misfit. 





1. INTRODUCTION 


LTHOUGH the arrangement of atoms in a lattice 

of a single metallic crystal is rather well under- 
stood, our knowledge concerning the structure of the 
intergranular region in a polycrystalline metal is still 
rudimentary. Since the contiguous grains have different 
orientations, there must be in the boundary region 
some disturbance of the crystalline structure, the de- 
gree of which, it has been intuitively predicted, varies 
with the orientation difference of the two grains. In 
order to test the various ideas concerning the structure 
of the grain boundary, much effort has been directed 
towards a study of properties of individual grain bound- 
aries as a function of the relative orientation of the 
adjacent grains. Yannaquis and Lacombe! have demon- 
strated the effect of orientation difference in increasing 
the rate of intergranular etching attack. Beck ef al.” 
have concluded that the rate of intergranular self- 
diffusion is influenced by the relative orientation from 
the finding that, during recrystallization, grain bound- 
ary migration is fastest when the orientation difference 
is the largest. Important quantitative evidence con- 
cerning grain boundary structure, presented by Dunn 
et al.** and Aust and Chalmers,® is the dependence of 
grain boundary energy on orientation: at small misfit 
angles the energy rises sharply with increasing orienta- 


tion difference and ceases rising at a relatively small ' 


angle. 

These and many other investigations of grain bound- 
aries have indicated that the properties of the inter- 
granular region depend on the relative orientation of 
the neighboring grains, differing more from the bulk 
material as the misfit angle increases. In the present 


* Now at Westinghouse Atomic Power Division. 

1 N. Yannaquisand P. Lacombe, Compt. rend. 226, 498-9 (1948). 

2 Beck, Sperry, and Hu, J. Appl. Phys. 21, 420-25 (1950). See 
also R. Smoluchowski, Phys. Rev. 83, 69 (1951). 

3C. G. Dunn and F. Lionetti, Trans. Am. Inst. Mining Met. 
Engrs. 185, 125 (1949). 

‘Dunn, Daniels, and Bolton, Trans. Am. Inst. Mining Met. 
Engrs. 188, 1245 (1950). 

6K. T. Aust and B. Chalmers, Proc. Roy. Soc. (London) 201, 
210-215 (1950). 


research, the rate of intergranular diffusion of Ag in 
Cu was studied as a function of the relative orientation 
of grains. The rate of diffusion indicates the degree of 
mobility of atoms and allows a comparison of “the 
tightness of fit” along various grain boundaries. This, in 
turn, permits the evaluation of the proposed models of 
grain boundaries. 


2. GRAIN BOUNDARY MODELS AND DIFFUSION 
MECHANISM 


In the dislocation model first proposed by Burgers, 
the mosaic boundary is pictured as an array of edge 
dislocations, each one formed by the termination of a 
plane of atoms. Read and Shockley,’ in their quantita- 
tive development of the model, which is applicable for 
boundaries separating grains of small orientation dif- 
ference, considered the total strain energy to be the 
sum of the elastic energies of the dislocations, the 
density per unit area depending on the relative orienta- 
tion. The general form of the energy versus orientation 
curve just described was explained as follows: at small 
angles the dislocations are far apart and all the strain 
energy is concentrated in the grain boundary and it 
increases with @. At larger angles, commonly considered 
greater than approximately 15°, the dislocation density 
is large, the model is no longer valid and plastic de- 
formation is distributed throughout the volume of the 
grain. The strain energy then no longer increases with 0. 

In his theory of viscosity of grain boundaries, Mott' 
proposed that grain boundaries are made up of islands 
of good fit separated by areas of disorder. He was suc- 
cessful in explaining on that basis the magnitude of the 
observed viscosity. Both these theories will be dis- 
cussed later in relation to the results of the present 
investigation. 

There are several possible mechanisms of volume 
diffusion: vacancy, simultaneous interchange of two or 
more atoms, interstitial, etc. Huntington and Seitz’ 


6 J. M. Burgers, Proc. Phys. Soc. (London) 52, 23 (1940). 


7 W. T. Read and W. Shockley, Phys. Rev. 78, 275-289 (1950). 
8N. F. Mott, Proc. Phys. Soc. (London) 60, 391 (1948). 


*H. B. Huntington and F. Seitz, Phys. Rev. 61, 315-25 (1942). 
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DIFFUSION IN GRAIN 
have shown that, at least in the case of copper, the 
first-mentioned alternative is the most probable. If we 
consider diffusion along a grain boundary, then the dis- 
tinction between the various mechanisms becomes less 

ronounced. For instance, in a distorted lattice it may 
. difficult or even impossible to distinguish between 
an interstitial position and a vacancy. In further dis- 
cussion we shall assume that grain boundary diffusion 
occurs by means of a vacancy mechanism, and we shall 
interpret our results from that point of view. In general, 
it is to be expected that with increasing misfit between 
two grains there will be more imperfections along grain 
boundaries favoring diffusion. 


3. EXPERIMENTAL PROCEDURE 
a. Choice of Material 


Inasmuch as the general grain boundary has five 
degrees of freedom, three to specify the orientation 
difference and two to fix the grain boundary direction, 
some kind of preferred orientation is necessary to re- 
duce the number of experimental variables. If, for 
example, a columnar structure is employed wherein all 
the grains have a common cubic direction and the 
grain boundaries are parallel to this direction, the 
degrees of freedom are reduced to two. Thus as seen in 
Fig. 1, a stereographic projection of a two-grain con- 
figuration on a plane perpendicular to the columnar 
direction, the relative orientation is completely fixed by 
the values of a; and ae, which are the angles between the 
grain boundary and the cubic directions of the grains, 
so taken that the sum of ai+a2, which equals @, is 
less than 90°. The complete relative orientation can be 
described by @, the orientation difference, and one of 
the a’s to fix the boundary direction. 

Since columnar Cu is readily available and since 
Berman and Harrington’® have shown that preferential 
intergranular diffusion of Ag in Cu can be demon- 
strated metallographically, this system was chosen for 
study. Columnar deoxidized copper was used of the 
following purity : Fe-0.0003 percent, Sb-0.0006 percent, 
Pb-0.0001 percent, Sn-Nil, Ni-0.0007 percent, Bi- 
0.00001 percent, As-0.0005 percent, Te-0.001 percent. 
The diffusion couples consisting of the above Cu and 
Cu with 2.4 atomic percent of Ag were welded in 
hydrogen in such a way that the interface was per- 
pendicular to the columnar direction and diffusion in 


each grain occurred in crystallographically equivalent 
directions. 


b. Precipitation Method 


In the first experiments the couples, after diffusion, 
were cooled slowly from the diffusion temperature so 
that the precipitate of Ag in Cu would delineate the 
diffusion zone after etching with NH,OH+H,0,. In 
Fig. 2, a section perpendicular to the interface, the 


”F,. Berman and R. H. Harrington, Trans. Am. Soc. Metals 
36, 362 (1945). 
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Fic. 1. Definition of angles used in relative orientation meas- 
urements. ¢ {100} poles as plotted. ¢(001) poles rotated into the 
center of the net. The diametric line represents the trace of the 
grain boundary in a plane perpendicular to the columnar direction. 


fast diffusion along grain boundaries is readily ob- 
served. The apparent faster volume diffusion in one 
grain than in the others, however, is not real. The 
difference is caused by an orientation etch effect as 
demonstrated by the fact that the diffusion zone is 
darker and its extent greater when the {100} planes of 
the grain are parallel to the surface of polish, an ob- 
servation in accord with the finding of Barrett, Kaiser, 
and Mehl" that the Ag precipitate platelets form on the 
{100} planes of the Cu matrix. The etch is, in other 
words, more sensitive when the plates of precipitate 
are parallel to the surface of polish. This suggests also 
that reported instances of anisotropic diffusion in cubic 
lattices may in reality be caused by a similar orientation 
etch effect. 

To eliminate this effect of orientation, the method 
was changed: The diffusion couples were analyzed by 
observing sections cut parallel to the interface, thus 
perpendicular to the columnar direction. On such sec- 
tions the (100) planes in all grains are parallel to the 
surface of polish and the orientation etch effect is 
eliminated. In Fig. 3, a section 0.05 mm from the inter- 
face, the precipitate is evenly divided between both 
sides of the grain boundaries as expected. 

It is to be noted that the grain boundaries along 
which there is Ag precipitate do not appear as isolated 
segments, but are in general connected in series. This 
appearance is apparently caused by the mode of freez- 
ing of columnar structures by which the grains solidify 
in colonies or clusters of similar orientation. This sub- 
ject will be discussed more fully in a forthcoming 
publication. 


4 Barrett, Kaiser, and Mehl, Trans. Am. Inst. Mining Met. 
Engrs. 117, 39 (1935). 
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Fic. 2. Diffused 260 hours at 675°C. Slow cooled. 
NH,OH+H.O, etch. 100. 


c. Solid Solution Method 


The drawback of the above precipitation method was 
the complicated influence of the rate of cooling on the 
amount of precipitate as determined by the shape of 
the solvus line of silver in copper. A search was made 
for a method which would reveal differences in silver 
content in a solid solution of silver in copper. Since 
silver is more noble than copper and is etched more 
slowly, it was concluded that at a lower temperature 
this difference will be still more accentuated. Indeed, 
by lowering the temperature of the etchant, 33 percent 
nitric acid, to 2-4°C the etching velocity was greatly 
reduced and a clear outline of the silver rich areas ob- 
tained. Admittedly, the actual concentration of silver 
along this “etching limit” is not known, but it is very 
well reproducible and makes possible the necessary 
comparisons between the various grain boundaries. 

The above method being much superior to the pre- 
cipitation method, all diffusion couples were now rapidly 
quenched from the diffusion temperature so as to re- 
tain the state of solid solution. They were then subse- 
quently immersed in the cool etchant, as described 
above, and the diffusion zones measured. This is illus- 














Fic. 3. Same as Fig. 2. 0.050 mm from interface. 50X. 
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trated in Fig. 4. These photographs show also how the 
course of diffusion may be followed along individual] 
grain boundaries indicating the decreasing amount of 
silver at increasing distance from the interface. Since, 
as yet, no good method has been developed to measure 
quantitatively the diffusion coefficient along an jp- 
dividual boundary,f the distance from the interface at 
which silver could no longer be detected metallo- 
graphically was denoted as the “distance of penetra- 
tion,” and plotted as a function of orientation. 

The experimental procedure consisted of lapping 
consecutive sections parallel to the interface at intervals 
of roughly 0.025 mm and etching and photographing 
them repeatedly until the bright bands at the grain 
boundaries could no longer be detected. The distance 
of penetration was arbitrarily taken as the distance 
halfway between the last section where the Ag still 
could be detected and the next section. Thus the error 
is about +0.015 mm. 

The relative orientations were determined by Laue 
back reflection patterns for each grain, plotting the 
{100} poles of adjacent grains on a stereographic pro- 
jection, and measuring the orientation difference angle, 
6, as defined in Fig. 1. 


4. RESULTS 


Diffusion runs were carried out at three temperatures, 
673°, 700°, and 725°C, the durations of which were 
calculated by means of the formula 


TiT>2 ty x2 


Q=R re 
Ti- T2 te x; 








to get approximately the same amount of intergranular 
diffusion at each temperature. Here Q is the activation 
energy in calories per mole, x is the distance of intergran- 
ular penetration after the time ¢ at temperature 7. The 
value for Q used, 23,000 cal/mole, is that obtained by 
the approximate method outlined in a previous pub- 
lication.” 

The results of the measurements at 725°C are pre- 
sented in Fig. 5 in which the depth of penetration, 1, 
is plotted as a function of the angle 0. The somewhat 
irregular distribution of points disappears if only those 
pairs of grains are taken into account which deviate 
from the perfect columnar structure™ by less than 8 
degrees (i.e., their [100] directions are parallel within 
this limit). Curves drawn through the data show then 
an increase in diffusion rate with increasing orientation 
difference up to the maximum which occurs at 45°. 
According to these plots, the maximum extent of inter- 
granular diffusion 0.52 mm, is approximately 2} times 

t Work is now in progress at this laboratory to measure the 
diffusion coefficient along an individual boundary by means of 
radioactive tracers. 

2M. R. Achter and R. Smoluchowski, Physics of Powder Mdal- 
lurgy (A Symposium held at the Sylvania Electric Products 
Company) (McGraw-Hill Book Company, Inc., New York, 1951). 


( 13M. R. Achter and R. Smoluchowski, Phys. Rev. 83, 16 
1951). 
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DIFFUSION 


the extent of volume diffusion, 0.20 mm, but it is not 
larger than the latter until a minimum orientation 
difference angle in the neighborhood of 20° is reached. 

Although the effect of orientation difference of the 
grains is readily apparent, the same is not true concern- 
ing the influence of the direction of the grain boundary. 
There seems to be some indication that the diffusion 
rate is smallest when the grain boundary direction bi- 
sects the angle @ and increases as the direction deviates 
more from the bisector; there were, however, so many 
exceptions to this rule that more measurements, de- 
signed to investigate this effect, must be carried out 
before any definite conclusions can be drawn. 


5. DISCUSSION 


Before discussing the grain boundary structure based 
on the diffusion measurements described above, it 
might be well to assess the reliability of the data. The 
relatively large distance, 0.025 mm, between the suc- 
cessive cuts of the sample makes the data not precise 
enough to come to any definite conclusion regarding 
the effect of temperature or of boundary direction on 
the orientation dependence of the diffusion rate except 
in isolated instances. There are, however, two results 
which stand out: the first is the absence of any measur- 
able preferential intergranular diffusion, at the em- 
ployed times and temperatures, between the angles of 
misfit of 0° to approximately 20° and approximately 
70° to 90°; the second is the increase in diffusion rate 
with increasing angle of misfit up to a maximum at 45°. 

As schematically pictured in Fig. 6, the diffusion 
data are quite different from the dependence of energy 
on orientation as obtained by other authors, a differ- 
ence which must be explained in terms of the effect of 
orientation on lattice matching or tightness of fit across 
the boundary. The explanation of the form of the energy 
curve based on Read and Shockley’s dislocation calcu- 
lation has already been given. The main features of the 
diffusion curve can also be explained tentatively by 
the dislocation model, assuming as previously stated, 
the vacancy mechanism of diffusion. The size of the 
void space associated with an edge dislocation is con- 
sidered to be approximately one-half of an atomic 
diameter and thus too small to play the role of atomic 
vacancy accelerating the passage of the Ag atom which 
is even larger than a Cu atom. The grain boundary 
diffusion is thus not faster than usual volume diffusion. 
When, however, at large angles the dislocations get so 
close together that they interact and loose their identity, 
then it is likely that such a “cluster” of dislocations 
provides the necessary void space to accommodate Ag 
atoms and the diffusion is greatly accelerated. 

It is also possible to describe the intergranular struc- 
ture in terms of a simple physical picture which is 
analogous to the usual theory of age hardening. Ac- 
cording to this theory, coherency strains between the 
growing precipitate particle and the matrix harden 
the material; when the strains become too high, co- 
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(a) 


(b) 





(c) 





Fic. 4. Sections parallel to the interface. (a) Diffused 194 
hours at 700°C and quenched. Etched by immersion in 33 percent 
HNO; at 2°C. 0.21 mm from interface. 8X. (b) Same field as in (a). 
0.36 mm from interface. 8X. (c) Diffused 260 hours at 675°C and 
slow cooled. Etched more heavily than (a) and (b). 0.20 mm from 
interface. 16x. 


herency is lost and softening results. In a similar manner 
one can argue that at small misfit angles the lattices of 
two grains are “coherent” across the grain boundary 
and the coherency strains, which increase with increase 
of orientation difference, are manifested as the meas- 
ured grain boundary energy. The drop in this energy 
at the critical angle indicates that when the strains be- 
come intolerable, coherency is partially lost. Some co- 
herency is, however, retained, since even at large angles 
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Fic. 5. Extent of diffusion, X, vs orientation difference @. 141 
hours at 725°C. — - — Extent of volume diffusion. Extent 
of intergranular diffusion. X (100) directions parallel within 8°. 
@ (100) directions not parallel within 8°. 
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Fic. 6. Schematic representation of the relative energy, E/Eo, 
and distance of intergranular penetration, X, as a function of the 
orientation difference, 6. The dashed line in the lower curve repre- 
sents the extent of volume diffusion. The maximum in the energy 
curve is indicated in the data of Dunn ef al. (see reference 4) and 
the theoretical calculations of Read and Shockley (see reference 
7). The plateau is indicated in the data of Aust and Chalmers (see 
reference 5). 


an appreciable amount of strain energy persists in the 
boundary. The persistence of a portion of the strain 
energy can be interpreted as meaning that part of the 
boundary is still coherent and part is not, an interpreta- 
tion which agrees with Mott’s conception that a grain 
boundary is constituted of islands of good fit (we would 
say coherent areas) separated by areas of misfit. The 
same line of reasoning explains the main features of 
the diffusion curves. At small angles the boundary is 
coherent resulting in good fit or good lattice matching, 
which prevents faster intergranular diffusion; in large 
angle boundaries, where the island model applies, 
faster diffusion can take place around the coherent 
areas, through the areas of disorder which may be 
identified with areas of high dislocation density. With 
increasing angle of misfit the proportion of areas of 
disorder increases and facilitates diffusion in accord 
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with experiment. This would also account for a drop 
in boundary energy due to the decrease of areas of 
coherency as implied by some experimental results,‘ 

The cusp indicated in Fig. 5 is also consistent with 
the idea that grain boundary structure is orientation 
dependent. At certain discrete angles, especially when 
the grain boundary is a rational plane of symmetry of 
the two grains, lattice matching should be especially 
good. The presence of the cusp at @ equal to 53°, a con- 
figuration corresponding to twinning on a (210) plane, 
the most rational plane possible in the (100) columnar, 
structure, is evidence of the improved degree of order 
in such a grain boundary. Such cusps have been dis- 
cussed by Read and Shockley.’ A further development 
of the two models of a grain boundary incorporating 
additional diffusion data will be given in a later 
publication. 

Based on the observed dependence of intergranular 
diffusion rate on orientation difference, as well as on 
energy data reported elsewhere, it can be concluded 
that grain boundaries can be classified into two types: 
small angle coherent boundaries in which the diffusion 
rate is the same as in the bulk of the grain, and large 
angle incoherent boundaries in which the amount of 
incoherency increases with the degree of orientation 
difference and reaches a maximum at 45°. 


6. SUMMARY 


1. It is shown that the rate of intergranular diffusion 
is strongly influenced by the relative orientation of the 
grains. 

2. The variation of the diffusion rate indicates a 
greater degree of disorder in grain boundaries with 
greater orientation difference. 

3. An interpretation of the diffusion data is given in 
terms of a grain boundary model which correlates the 
dislocation model applicable at low angles of misfit 
with the “island” model applicable at large angles of 
misfit. This interpretation is also in agreement with the 
dependence of grain boundary energy on the angle of 
misfit. 
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The electrical properties of selenium have been studied in their dependence on history of preparation, 
ie., as functions of temperature of initial liquid selenium, quench procedure, nucleation procedure, tem- 


perature of crystallization, and time of crystallization. 


The electrical properties of melt single crystals and the microcrystalline matrix of these crystals have 


been compared. 


The influence of nonmetals S, O2, I2, in small quantities on the electrical properties has been determined. 

Carrier densities and effective mobilities as functions of temperature of measurement have been estimated 
from data on thermoelectric powers and resistivities. The order of the carrier densities have been checked 
by measurements of the voltage dependence of capacity of completely vacuum prepared rectifiers. 

The frequency dependence of resistivity has been determined to gain insight into the nature of the ef- 


fective mobilities. 





INTRODUCTION 


HE work at the University of Pennsylvania on 
the electrical properties of selenium has been re- 
ported, in part, in references 1 and 2 which are con- 
cerned with the conductivity of the liquid, and the 
properties of single crystals from the melt, respectively. 
The present paper is written to describe the electrical 
properties of the most common form of selenium, the 
microcrystalline or polycrystalline hexagonal selenium 
and correlate the results with the conclusions reached 
in reference 2. 


MICROCRYSTALLINE SAMPLES 


The selenium used was the same as that described in 
references 1 and 2. Non-volatile impurities were in the 
order of parts, Fe, Te, Pb, Cu, in fractional parts per 
million. Volatile impurities were not analyzed. The 
method of purification was the following: 


1. Distillation of selenium at atmospheric pressure; 
2. oxidation in vapor state to SeO:; 3. resublimation of 
the SeO2, followed by reduction in the vapor state with 
NH; to elemental selenium; 4. distillation of the se- 
lenium in quartz apparatus followed by condensation 
and shotting in distilled water. 


The selenium was kept in a closed but not air tight 
glass vessel. The nonmetal impurities present were 
probably S, X2, P, O2, and Nz. The work on the liquid 
has shown that the selenium was free of any appreciable 
amounts of phosphorus and halogen. Furthermore, the 
purification procedure above would not have introduced 
these impurities. Sulfur, oxygen, and nitrogen were 
present to some extent. As described later, sulfur in 


*The major part of this work was conducted under contract 
NObsr 42487 with the Bureau of Ships, Department of the Navy. 
The remaining portion and continuing studies are under contract 
DA-26-039-SC-5426 with the Signal Corp. U. S. Army. 

'H. W. Henkels, J. Appl. Phys. 21, 725-731 (1950). 

*H. W. Henkels, J. Appl. Phys. 22, 916 (1951). 

*H. W. Henkels, Phys. Rev. 77, 734-736 (1950). 

*H. W. Henkels, “Summary of recent work on the conductivity 
of selenium,” NObs 34144 Bu Ships Tech. Rept. No. 9. 


amounts to 0.1 percent does not appreciably affect the 
electrical properties of crystallized selenium. Therefore, 
the possible important nonmetal impurities were O2 
and N2. 

The selenium was further prepared by heating molten 
in a glazed porcelain container in vacuum to about 
2X10 mm Hg. 


EXPERIMENTAL ARRANGEMENTS AND 
PROCEDURE 


Initial studies of crystallization rates were made by 
inserting flat nickel plate electrodes into samples of 
amorphous selenium cooling in pyrex dishes prepared 
from square cross section glass tubing, and noting the 
resistances as functions of time of heating at 120°C 
and 210°C. Checks of the range of resistivities were 
made with potential probes. No contribution of end 
contact resistance was noted but in what follows all 
data taken without potential probes will be presented 
in terms of total resistance. It was not feasible to make 
potential probe measurements in these studies because 
of the desire to continuously note the resistivities during 
the initial crystallization and because of the changes in 
density occurring in the process. 





Fic. 1. Experimental arrangement for measurement of 
thermoelectric power and resistivity. 
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Fic. 2. Hysteresis loops in the resistivity of selenium. 


In all other studies of dc resistivity potential probes 
of nickel were employed. Good current contacts were 
obtained with carbon suspensions and lead foils on the 
ends of the specimens. 

The device used in measuring thermoelectric power 
and resistivity simultaneously is shown in Fig. 1. The 
temperatures of the ends of the samples were measured 
by chromel alumel thermocouples which were insulated 
from the Ni blocks except at small points at the sur- 
faces of the blocks in contact with the samples. Thermal 
gradients were maintained by heating one of the Ni 
blocks. Good thermal contact was insured by screw 
pressure. The entire mount was situated in the center 
of a thermostated oven. 

High frequency measurements were made with con- 
ventional ac bridges and Q meters. The specimens were 
in all cases mounted between circular parallel plates to 
maintain the field lines parallel to the axis of the sam- 
ples and avoid the so-called Howe effect of capacita- 
tively shorting them. Checks of the effectiveness of the 
parallel planes were made with the active element of a 
carbon resistor of dc value 100K (Fig. 5). A further 
check was made by varying the cross sectional area of 
a selenium sample by two and observing that the re- 
sistance values varied inversely by the same factor. 

The resistivities of all samples were checked with 
potential probes and the over-all resistance between the 
parallel planes at zero frequency calculated from ge- 
ometry. Measurements were made as a function of fre- 
quency to 200 mc in the absence of light at atmospheric 
pressure. 

The simultaneous measurements of thermoelectric 
power and resistivity of the same specimens were con- 
ducted afterwards under conditions which are better 
described later with the results. 





HENKELS 


EXPERIMENTAL RESULTS 


In that the electrical properties of selenium depend 
in a complex manner on purity and microcrystalline 
structure, it was necessary to study them during al] 
stages of sample preparation. The resistivity of liquid 
selenium has already been reported in reference 2. [t 
was found that pure selenium has a single reproducible 
resistivity and temperature coefficient from the melting 
to the boiling points. The resistivity is a function of 
temperature alone. Selenium doped with nonmetals 
Clo, Iz, and P, exhibited lower resistivities and tempera- 
ture coefficients. 


A. Resistivity of Solid Hexagonal Selenium as a 
Function of the Quench Rate of the Liquid 


In preparing selenium specimens for study and in 
some commercial processes, liquid selenium is cooled 
to an amorphous black form which must then be crystal- 
lized by heating for periods of time at any temperature 
in the range somewhat above room temperature to the 
melting point. 

Molten selenium is considered as consisting of long 
chains of selenium atoms analogous to the case of sulfur 
chains. The length of these varies with temperature de- 
creasing at higher temperatures. The possibility of 
rapidly quenching nonequilibrium chain lengths of 
selenium into the amorphous material and producing 
thereby a conceivable change in the electrical properties 
presented itself. However, it was found that the tem- 
perature of the liquid from which the samples were 
poured did not appreciably influence the subsequent 
resistivities, although a slight effect was noted above 
500°C. 

B. Resistivity as Function of Nucleation Procedure 

Selenium crystallizes at any temperature from about 
room temperature to the melting point with a speed 


that depends on: (a) nucleation rate; (b) grain growth 
rate. 
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Fic. 3. Resistivity of selenium samples crystallized with 
out quench and crystallized at temperature of maximu® 
nucleation. 
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As already discussed in reference 2 maxima exist in 
the temperature rate functions. For nucleation, the 
maximum occurs at about 110° or 120°C, for grain 
growth, somewhat below the melting point. Further- 
more, the electrical properties depend on the nucleation 
and grain growth temperatures. It is this last fact 
which has lead in the past to the wide variation in re- 
ported values of resistivities and the postulation of 
different types of hexagonal selenium. Studies were 
made of the changes in resistivities produced by dif- 
ferent procedures of nucleation. The following is given 
in brief review. When the liquid is cooled rapidly to the 
amorphous form at room temperature, the selenium 
nucleates so that the speed of complete crystallization 
depends on the temperature as determined by the grain 
growth rate. Slow cooling of liquid selenium to about 
210°C results in very slow crystallization because of the 
slow nucleation rate. The electrical properties of samples 
nucleated in different manners but crystallized at the 
same temperatures are different. The resistivity-tem- 
perature behavior of selenium samples nucleated during 
quenching to room temperature and crystallized at 
higher temperatures, is very complex and dependent 
on temperature rates of change. Maxima occur at 
different temperatures. Figure 2 shows the typical be- 
havior of the thermoelectric power and resistivity of 
such samples. Samples nucleated by cooling slowly from 
the melt exhibit negative temperature coefficients but 
retain very considerable hysteresis loops in the tem- 
perature dependence of resistivity. The lower portion 
of Fig. 3 shows the behavior of such samples. The upper 
portion of the same figure demonstrates the behavior 
of samples crystallized at lower temperatures. In these 
cases, the hysteresis loops are very small. Thus, it is 
evident that structural defects in selenium are playing 
an important part in determining the resistivity. More- 
over, such defects are not relieved at any temperature 
below the melting point. 

Certain impurities accelerate the nucleation, the non- 
metals O2, Xs, having the largest effects. Phosphorus 
prevents crystallization. 


C. Electrical Properties as Function of 
Crystallization Temperature and Time 


The previously mentioned factors which influence the 
resistivity of any selenium sample are of secondary 
importance when the effects of nonmetal impurities 
and crystallization temperature (as already indicated 
in Fig. 2) are considered. The period of crystallization 
is also of considerable importance. On crystallization 
the resistivity drops rapidly to a minimum then rises 
on subsequent heating at a rate depending on the purity, 
crystallization temperature, and nucleation procedure. 
Figure 4 shows the changes in resistance occurring dur- 
ing crystallization at temperatures of 120° and 213°C 
The small effect of liquid temperature mentioned above 
is also noted in the same sheet. 


More recent work has aimed at determining: 


1. The changes in the fundamental electrical prop- 
erties of selenium (mobilities, carrier densities, activa- 
tion energies) producing the effects noted (i.e., the de- 
pendence of resistivity on temperature and time of 
crystallization), and 

2. The variation in the former properties caused by 
the presence of small amounts of nonmetal impurities 
in the selenium. 


The resistivity in all cases depends upon the fre- 
quency. It was therefore necessary to examine this de- 
pendence concurrently. The thought was to correlate 
the fundamental properties with the frequency de- 
pendence of resistivity. The high frequency work also 
served to show any macroscopic defects in the samples. 





Fic. 4. Changes of resistance of selenium during crystallization. 


In preparing final samples for a study of item 1 pure 
selenium was heated molten five hours in vacuum to 
2X 10~* mm Hg. Samples were poured in air from liquid 
selenium at 350°C into clean glass boats at room tem- 
perature. The liquid cooled in about five minutes to the 
black amorphous condition. Glass cover plates were 
placed over the samples. The glass boats and covers 
were then inserted in a light press and the press in a 
thermostated oven. As the samples crystallized, the 
press formed the samples to flat bars of approximate 
dimensions 0.18X1.3X3 cm. Two series of samples 
were crystallized at temperatures of 120°C and 208}°C. 

In these series the period of crystallization was varied. 

One selenium boat without cover glass was held 
somewhat above melting point for 4 hours, then cooled 
in thirty minutes to 205°C and held at that tempera- 
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Fic. 5. Frequency dependence of resistivity of various selenium samples. 


ture 41 hours. This and another from the matrix of the 
selenium in which the single crystals were grown, were 
cut to roughly the same shapes as the other samples. 


Crystallization at 120°C 


Figure 5 shows the frequency dependence of re- 
sistance of samples (a) to (d). The minimum in the re- 
sistivity as a function of time of crystallization has been 
mentioned above. Samples (a), (6), (c) represent con- 
ditions before the minimum, (d) the condition beyond 
the minimum. 





Fic. 6. Temperature dependence of resistivity of selenium samples 


crystallized at low temperatures for different periods of time. 


The thermoelectric powers and resistivities of these 
samples were measured at room temperature. The speci- 
mens were then rapidly heated and observations made 
of the temperature dependence of resistivity from room 
temperature to about 120°C. The results are shown in 
Fig. 6. 

The thermoelectric power of sample (d) exhibited a 
maximum about 97°C. This maximum has been previ- 
ously found in the data for similarly crystallized samples 
in reference 3. 

The resistivity-temperature curves also exhibit in- 
flection points as previously noted in the above 
reference. 

The hole densities, calculated from the thermoelectric 
power were about the same for the various samples, 
6 or 7X10"*/cm’, at 40°C. With these values, the mo- 
bilities increase to a maximum about 0.02 cm?/volt sec, 
at the minimum in the resistivity crystallization time 
curve and slowly decrease on further heating. The 
mobilities increased roughly exponentially with tem- 
perature of measurement. 


Crystallization at 210°C 


The lower portion of Fig. 5 shows the frequency de- 
pendence of series (e) to (k). It was observed that, in 
the cases of large end contact resistance, the resistivities 
decreased radically from 100 kc to 1 mc to values repre- 
senting the bulk properties of the samples. The range of 
values of total bulk resistivity as indicated by the lower 
lines coincides with that obtained later in resistivity 
measurements in conjunction with the thermoelectric 
power measurements. Thus the samples were homo- 
geneous on a macroscopic scale. All samples decreased 
in resistivity continuously at higher frequencies. 


The resistivity temperature curves are shown I ' 


Fig. 7. Inflection points around 84°C were apparent in 
the incompletely crystallized samples. The continual 
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decrease in slope with time of crystallization is in- 
teresting. 

The thermoelectric power for sample (i) is shown in 
Fig. 8 together with values of resistivity, carrier den- 
sity, and mobility. The density of carriers in the series 
(e) to (h) at 40°C averaged 5X 10'*/cm*. This value is 
identical with the 6 or 7X 10'°/cm# obtained in samples 
(a), (b), and (d), within the error of the measurements. 

The mobilities increased steadily with longer periods 
of crystallization but no sample was selected which had 
been crystallized a sufficiently long time to pass the 
minimum. The range of mobilities was 0.12 cm?/volt sec 
to 0.38 cm?/volt sec from (e) to (k). However, even with 
sample (i) having a constant resistivity, the sharp in- 
crease in thermoelectric power with temperature re- 
sulted in a decrease in » and a consequent apparent 


At temperctures designated with arrows the carrier 
densities (Cm-3) were: 


(e) 46 x 10'> 
(f) 5.3 x10'5 
(9) 6.7 x 10'> 
(h) 5.2 x 10'> 





Fic. 7. Temperature dependence of resistivity of selenium samples 
crystallized at high temperatures for different periods of time. 


rapid increase in mobility with temperature. It should 
be noted, in this connection, that all the curves show 
inflection points about 125°C. 

Figures 6 and 7 together show several points. First, 
the samples of the same lot and melt of selenium all 
had the same density of carriers (~5 or 6X 10'°/cm*) 
at 40°C. Second, the mobilities at that temperature in- 
creased steadily on crystallization in series (a) to (d) 
and (e) to (h). 


Crystallization at 220°C—Single Crystal Matrix 


The properties of single crystals of selenium grown 
in a melt were discussed in (2). A sample was cut from 
the matrix of these crystals and designated sample (/). 
The frequency behavior of this sample is shown dashed 
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in Fig. 5. There was no great difference in the fre- 
quency behavior of the high resistivity samples (d), 
(a) and (j). 

The thermoelectric power was somewhat lower in- 
dicating a hole density of about 9X10'*/cm* at 40°C. 
The mobility was very low, approximately 0.015 cm?/ 
volt sec at the same temperature. The single crystals 
grown in this matrix, on the other hand, had hole 
densities about 2X10"/cm* and mobilities Lc axis 
about 0.06 cm?/volt sec, || ¢ axis 0.18 cm?/volt sec. 
Thus the density of carriers in the microcrystalline se- 
lenium was a factor of 30 or 40 higher than in the single 
crystal while, on the other hand, the mobility was lower 
by a factor of about 6. The “activation energy” was 
greater than for sample (/). 


D. Electrical Properties as Function of 
Nonmetal Impurities 


Selenium of melt A was reheated in vacuum to 
4X 10~ mm Hg for 4 hours in a flat dish. Two samples 
were poured in glass boats and cooled rapidly to room 
temperature. The remaining selenium was divided into 
three portions and cooled rapidly. The portions were 
doped individually with S, Se2I2(?), and SeO2 and stirred 
above the melting point for a short period of time. 
Samples were poured in glass boats at room temperature 
and subsequently crystallized at the temperature of 
120°C and 210°C for periods of time, 180 minutes and 
60 minutes, respectively (consult tables for desig- 
nations). 





Fic. 8. Thermoelectric powers, resistivities, mobilities, and carrier 
densities for a sample crystallized at high temperature. 
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Fic. 9. Frequency de. 
pendence of resistivities 
of selenium samples con- 
taining impurities of 
oxygen, sulfur, and jo- 
dine, crystallized at low 
temperatures. 
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Crystallization at 120°C 


Figure 9 shows the behavior of samples (/), (7), (p), 
(r), (uw), and (w) as function of frequency. It is noted 
that all the higher resistivity samples begin to de- 
crease in resistivity at a lower frequency, approxi- 
mately 10 ke. 

Figures 10 and 11 present typical data on the re- 
sistivities, thermoelectric powers, carrier densities, and 
mobilities of the same specimens. The following factors 
were noted from an examination of the data. 


1. The resistivities of all samples were approximately the same 
at room temperature. 

2. The “activation energies” were the same for all samples ex- 
cept («) which has a higher one. 

3. The thermoelectric power exhibited a maximum about 84°C. 

4. The mobilities varied as in Table I. 

5. The carrier densities varies as in Table IT. 


The density of carriers in the vacuum prepared 
sample was the same as in the series (a) to (i). In- 
creasing amounts of O2 and S lower the thermoelectric 
power somewhat and consequently an increase in hole 
density is calculated. A small amount of iodine impurity 
increased the hole density slightly but a larger amount 
decreased the density below that in the vacuum pre- 
pared sample (/). In all cases the mobilities increased 
sharply with increasing temperature of measurement. 











TABLE I. 
Impurity Mobility Mobility (max) 
Sample (atomic %) (40°C) (cm?/volt sec) 
(1) ee 0.010 0.060 
(n) 0.0337 O 0.009 0.035 
(p) 0.00566 I 0.011 0.11 
(u) 0.0695 I 0.024 0.20 
(r) 0.0338 S 0.007 0.04 
(w) 0.233 S 0.003 0.02 








Crystallization at 210°C 


Figure 12 shows the frequency dependence of samples 
(m), (0), (q), (s), (4), (v), and (x) crystallized at 210°C 
for one hour. Samples (a), (6), (c), (d), (2), (n), (), (u), 
(v), (w), of high “activation energy” are similar to (i), 
(7), (D), (m), (0), (9), (s), , (2), (~), of low “activation 
energy.” The behavior of samples containing sulfur was 
identical with that of vacuum prepared selenium. The 
lower resistivity samples were constant with frequencies 
to almost 50 mc. 

Figure 13 presents the data on temperature de- 
pendence of resistivity, for samples (m), (s), («). Similar 
data for specimens (m), (0), (¢) are given in Fig. 14 while 
Fig. 15 gives data on (m), (gq), and (2). 

The dependence of thermoelectric power on tempera- 
ture was just about the same in samples (m), (s), (x), 
(v), and (¢). The thermovoltage of (0) and (q) increased 
with increasing temperature but at a slower rate. This 
data is assembled in Fig. 16. 

Table III gives carrier densities at 40°C, 100°C and 
190°C for the samples. 

Table IV gives the corresponding mobilities. 

The “activation energies” were about the same for 
samples (m), (0), (q), (s), («). Larger amounts of im- 
purities of O. and I, produced samples which increased 
in resistivity with temperature. 

Figure 17 assembles data on carrier densities. 


DISCUSSION 


General 


Selenium is generally considered as a semiconductor 
because of its position in the table of the elements and 
because of its photoelectric and contact properties. 

In explaining the contact properties of rectifier 
selenium, Schottky and co-workers have employed 2 
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Fic. 10. Mobilities and hole densities for samples containing 
oxygen and iodine crystallized at low temperatures. 


model of the semiconductor in which the activation 
energy for hole conduction (arising from impurity levels) 
was very small. This construction was developed from 
and for the analysis of the behavior of rectifier contacts. 

He also developed, on the basis of the assumption 
stated, a method for obtaining the density of carriers 
in the contact boundary layer from measurements of 
capacity. Values of carrier densities between 10'°/cm 
and 10'*/cm were noted. 

A survey of the literature on the resistivity of selen- 
ium to 1947 is given in reference 4. Briefly it had been 
found that it was possible to obtain grey hexagonal 
selenium with resistivities varying by a factor of 10* 
and activation energies positive and negative up to high 
values. In addition, the observed “activation energies” 
are not constant but vary with many factors. The 
resistivity data alone has been insufficient to separate 
the temperature dependences of » and b and permit the 
desired evaluation of hole densities and mobilities. 

It is possible to obtain values of n from measurements 
of thermoelectric power and Hall effect. In the present 
report 7 has been evaluated from data on the former. 
The relation used was 


Q=k/e(log RT!—5.32), 


where R is the Hall constant. This applies to an ideal 
case of free electrons in a structureless isotropic medium. 
The same objection is raised to the interpretation of the 
Hall constant. On the other hand, reasonable results 
have been obtained in the past with other semicon- 
ductors. In addition, there appears to be no other gen- 
eral way of arriving at values of m, other than by meas- 
urement of contact capacity. As is evident from the 
statement above and the data presented under experi- 
mental results, the two methods agree as to order of n. 
Plessner® has recently stated better than order of mag- 





°K. W. Plessner, Nature 166, 1073 (1950). 
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Fic. 11. Thermoelectric powers and resistivities of samples 
containing iodine crystallized at low temperatures. 


nitude agreement between carrier densities of single 
crystals based on his Hall effect measurements and our 
own work on thermoelectric power. 

The possible sources of errors in the three procedures 
for determining m are not the same. 

We have prepared selenium rectifiers completely in 
vacuum (including heat treatment and application of 
counterelectrode) from the pure selenium discussed 
previously and obtained results for acceptor densities 
from measurement of the voltage variation of capacity 
between 10 and 10'*/cm*. It is thought that the 
effects of oxide surface layers were decreased in our 
procedure. The agreement of the values is satisfying 
but the meaning may perhaps be still somewhat ques- 
tionable in view of the possible presence of surface 
conditions different from those in the bulk material. 

Report 2 has presented the properties of single 
crystals. It was noted that the density of carriers as 
determined from thermoelectric power measurements 
decreased with increasing temperature and that the mo- 
bilities increased exponentially with temperature. Fur- 
thermore, since the behavior of the resistivities with 
field may be equally well explained by potential bar- 
riers influencing the mobilities and the carrier densities, 
it was not possible to choose between the two. In work- 
ing with microcrystalline selenium, there is more lati- 
tude in the variation of the crystallization procedures 


TABLE II. 








Hole density 40°C Hole density (min) 





Sample Impurity (atomic %) (cm~) 
(1) han 6X 1015 3X 105 
(n) 0.0337 O 8x 10" 7.5X 10" 
(p) 0.00566 I 7.5X 10% 4x 10 
(u) 0.0695 I 2.3X 1015 1.5 105 
(r) 0.0338 S 7X10" 4x 10" 
(w) 0.233 S 15x 10'5 7X 10% 
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and impurity contents. In the case of microcrystalline 
selenium, in view of the past successful application of 
the idealized classical relations and the good agreement 
of the only available means of determining n, there is 
reason to consider the results valid, at least in indicating 
general behavior. It will be the purpose in the following 
discussion to present first a reasonable explanation of the 
properties of the microcrystalline aggregate and second 
to correlate the results with those in reference 2. 


1. General Behavior of Electrical Properties As 
Function of Period and Temperature 
of Crystallization 


Series (a) to (i) has shown the general behavior of the 
samples in the process of crystallization at two tem- 
peratures. It has been determined by past investigators 
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Fic, 13. Resistivities of samples with sulfur crystallized 
at high temperatures. 


that the high resistance selenium transforms to the 
lower resistance forms (crystallized at higher tempera- 
tures) on heating at the higher temperature. Figure 7 
shows well the gradual decrease in the activation energy 
with time of crystallization at any temperature (to the 
minimum in the crystallization time curve). 

In the entire series the carrier densities remained 
approximately constant with perhaps a slight decrease 
with the crystallization at higher temperatures. 

The mobilities, on the other hand, increased by about 
two orders in the series (a) to (#) at room temperature, 

In sample (), with which a larger range of tempera- 
tures was possible, the acceptor density n decreased 
with temperature. Thus the combination of increasing 
b and decreasing n resulted in a resistance independent 
of temperature. In case of samples crystalilized at low 
temperature, n also appeared to decrease with tempera- 
ture but the temperature range was too restricted to 
show a large change. 

The mobilities in both cases increased with tempera- 
ture but the rate of increase was higher for the high re- 
sistivity sample so that at 125°C the mobilities were 
approximately the same. 

The independence of resistivity of sample (/) should 
be reexamined. A skeptical view of this data alone would 
assume that the application of the thermoelectric power 











TABLE III. 
Impurity Carrier density (1/cm‘) 

Sample (atomic %) 40°C 100°C 190°C 
(m) te 2x 10'5 1.2 10" 1X 10" 
(s) 0.0338 S 2.1 105 1.2 105 1.1 10" 
(x) 0.233 S 3X 10'5 1.8X 105 2x 10" 
(0) 0.0337 O 4.3X 10 1.1 1015 3x 10" 
(t) 0.293 O 4.1 1016 1.3X 105 7X10" 
(q) 0.00566 I 7X10" 3X 10" 1.5X 10% 
(v) 0.0695 I 1.1 10'6 4.2X 10'5 4x 10" 
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Fic. 14. Resistivities of samples with oxygen crystallized 
at high temperatures. 


is only valid for order of magnitude and n is actually 
constant. This assumption together with one of b con- 
stant would give the desired result. However, in view 
of the other curves and our data in past reports, this 
view of b is untenable. Furthermore, Fig. 5 shows that 
the resistivity decreases above 10 mc. It would be 
expected that with a further increase in room tempera- 
ture mobility and a corresponding decrease in rate of 
change of mobility with temperature, selenium samples 
would attain positive temperature coefficients of re- 
sistivity if n actually decreased with temperature. 
That this actually occurs was shown in the series of 
measurements with samples having added impurities. 
In addition, past work has shown the possibility of 
positive temperature coefficients. 

It thus appears that the complex resistivity of selen- 
ium results from the counteracting influences of ac- 
ceptor densities which decrease with temperature above 
room temperature, and mobilities which increase with 
temperature. Furthermore, after crystallization at any 
temperature and for any period of time, the hole densi- 
ties at room temperature are constant. The decreases 
in resistivity on crystallization and the concurrent de- 
creases in “activation energy” are caused almost en- 
tirely by large increases in mobility and decreases in 
temperature rate of increase of mobility. 

To our knowledge, there is no theory of the behavior 











TABLE IV. 
Mobility (cm?/volt sec) 
Sample Impurity (atomic %) 40°C 100°C 190°C 
(m) ee 0.033 0.10 1.8 
(s) 0.0338 S 0.044 0.14 2.4 
(x) 0.233 S 0.022 0.07 1.0 
(0) 0.0337 O 0.058 0.35 2.2 
(t) 0.293 O 0.11 0.3 2.2 
(q) 0.00566 I 0.052 0.24 0.86 
(v) 0.0695 I 0.7 1.7 8.4 











Fic. 15. Resistivities of samples with iodine crystallized 
at high temperatures. 


of the ideal lattice which predicts an exponential in- 
crease of mobility with temperature. The theory ap- 
plicable to the other elemental semiconductors predicts 
a decrease of impurity scattering with increased tem- 
perature at low temperatures, and an increase of lattice 
scattering at higher temperatures. Thus the mobility 
would be expected to decrease at temperatures near the 
melting point. On the other hand, if the acceptor levels 
are occasioned by lattice defects as indicated in report 
2 and the remarkable constancy of hole densities shown 
in this report, there may be reason to expect the per- 
fecting of the lattice at higher temperatures resulting 





Fic. 16. Thermoelectric powers of various samples, 
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Fic. 17. Hole densities of various samples. 


in a decrease in n. The work of Brill and Krebs® gives 
considerable support to this view. Since the exact 
nature of the impurity levels in selenium is unknown, 
the effect of such a perfection of lattice on the mobility 
cannot be determined. However, it may be reasonable 
to expect that in this temperature range lattice scatter- 
ing predominates so that the scattering from defects 
would be unimportant. In addition, the data just 
discussed shows that with n constant, b varies by orders 
so that defect scattering to account for the mobility 
and its temperature dependence seems out of the ques- 
tion. Since lattice scattering causes a decrease in mo- 
bility with temperature, the scattering, to give rise to 
the noted temperature dependence (exponential in- 
crease in mobility with temperature), must arise from 
inhomogeneities in the specimens. 

The crystallization process should be reviewed at 
this point. In the initial stages of crystallization, small 
crystallites of selenium are enclosed in a matrix of 
amorphous selenium. On further crystallization, the 
grains grow at the expense of the selenium at the grain 
boundaries. The rate of growth of the grains is a func- 
tion of temperature with a maximum at about 210°C. 
The degree of crystallization at any temperature cannot 
be determined with sufficient accuracy to enable calcu- 
lation of the extent of the disordered regions between 
the crystallites. It is evident that the extent of the dis- 
ordered regions is a function of the time of crystalliza- 
tion. From the similarity in appearance of the resis- 
tivity temperature curves for samples incompletely 
crystallized at high temperatures to those of samples 


* R. Brill and H. Krebs, “First interim report of results on the 
structure analysis of selenium rectifiers,” PB 21981. 
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crystallized for long periods at lower temperatures, it 
may be hypothesized that the extent of the disordered 
regions is also a function of the maximum temperature 
(below 210°C) to which the samples are raised. The 
fact that samples crystallized at lower temperatures 
become similar on being heated to a temperature T to 
those initially prepared at that temperature seems to 
indicate an “increased crystallization.” 

The constancy of the carrier densities during the 
crystallization process should be kept in mind. Our 
work has been conducted at two temperatures of crystal- 
lization 110°C and 210°C but the forms of resistivity 
curves for crystallization at any intermediate tempera- 
ture (work of Borelius ef a/.’) are regular. The “activa- 
tion energies” continuously decrease with increased 
temperature of crystallization. It is not considered that 
these disordered regions are of large enough extent to 
affect any density measurements. 


2. General Behavior of Hole Densities with 
Crystallization Conditions and Nonmetal 
Impurities at Fixed Temperature 


Consideration of Microcrystalline Nature 


In reports (2) and (3) we have noted the hole densi- 
ties in single crystals prepared from a melt of selenium. 
These were remarkably low (~10"/cm’ at 40°C) in 
view of the expected purity of our selenium. The con- 
stancy of the resistivity of selenium single crystals 
prepared from different sources and in different man- 
ners, with other evidence in (2) lead us to consider the 
acceptor levels as probably arising from lattice defects. 
An alternate explanation was given that large crystals 
of selenium grow with a very small number of impurity 
acceptors of oxygen or nitrogen, which number is inde- 
pendent of the manner of crystallization or purity of the 
selenium. As it is not meaningful to compare the re- 
sistivity of vapor crystals with the melt from which they 
were prepared, the growth of crystals in the melt pre- 
sented the opportunity of comparing the matrix of the 
material from which the crystals grew with the crystals 
themselves. 

The data have been given for sample (/). It is to be 
noted that the hole densities were a factor of about 30 
higher. The cause of the difference may be sought in 
the following. 


(1) The crystals in growing tend to purify themselves so that 
those last formed or at least the exterior surfaces of the crystals 
have higher acceptor densities. 

(2) The lattice defects are increased by reason of the influence 
of the discontinuities at the grain surfaces. The grain size is very 
small so that such an effect might be appreciable in microcrystal- 
line selenium but unnoticed in single crystals. 

(3) The thermovoltage measurements are influenced by the 
granular structure of the material. 


The last possibility was considered first. The previ- 
ously mentioned agreement of capacity and thermo- 


7 Borelius, Pihlstrand, Anderson, and Gullberg, Arkiv. Mat. 
Astron. Fysik 30A, No. 14 (1944). 
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ELECTRICAL PROPERTIES OF SELENIUM. II. 


electric power measurements was some support for dis- 
carding this view point. In addition, most of the 
acceptor densities noted in this report fall in the range 
10" to 10'*. Thus unless we make the additional assump- 
tion that only very small amounts of impurities (less 
than we added, i.e., 0.0337 atomic percent oxygen) are 
effective in altering the influence of the irregular grain 
boundary regions on the thermoelectric power (in- 
creasing the calculated hole densities a factor of thirty), 
while increased amounts of impurities, in the range of 
our preparations, produce very little changes in the 
thermoelectric power while having an extremely large 
effect on resistivities, possibility, (3) must be discarded. 

The other two considerations lead to the same models 
of the microcrystalline aggregate. Thus considering 
only the series with no intentional impurities the in- 
creased hole density arose from increases in the accep- 
tors as occasioned by the influence of discontinuities at 
grain surfaces or by oxygen or nitrogen impurities in the 
surface layers. 


Consideration of Nonmetal Impurities 


Under the last heading the constancy of hole density 


‘with extent of crystallization was emphasized. 


From Table II for crystallization at 120°C and refer- 
ence 3 it is noted that the hole densities increased 
slightly with added impurities of oxygen and sulfur. 
A very small increase with added iodine in amount 
0.00566 atomic percent was noticed but an increased 
amount of iodine produced an apparent decrease in 
hole density. Crystallization at 210°C produced speci- 
mens with a considerably decreased hole density. In 
these samples the hole densities increased very little 
with added sulfur, but were increased with increasing 
amounts of both oxygen and iodine. 

The data in Tables II and III are not directly com- 
parable nor are individual items in the table. As noted 
previously, the crystallization rates are dependent on 
purity so that comparing samples with the same periods 
of heat treatment omits the factors considered in the 
series of samples (a) to (h). 

All the data in Tables II and III for samples crystal- 
lized at 210°C showed the regular increase in hole den- 
sity with added impurity. As previously referred to in 
reference 2, the addition of sulfur in amounts to 0.233 
atomic percent produced no appreciable changes in 
hole density. The addition of oxygen in amount 0.0337 
atomic percent increased the hole density a factor of 2 
at 40°C which would indicate that the impurity of 
oxygen present in our vacuum samples was at most 
less than this amount. If we assume for the moment 
impurity acceptor levels, we can calculate the efficiency 
of the impurities in producing acceptor levels (Table V). 
The efficiencies are very low and in view of the con- 
stancy of the resistivity of crystals give some support 
to the ideas of surface layers of impurities discussed 
previously. 

Thus in the cases of oxygen and iodine impurities the 
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hole densities seem to be definitely shifted upward but 
the amount of the shift is quite small. It is interesting 
to note that oxygen has about the same influence as 
iodine. The effect of oxygen on the electrical properties 
of selenium has not been properly recognized in 
literature. 


3. General Behavior of Hole Densities 
with Temperature 


Examination of the previous data on samples crystal- 
lized at low temperatures shows that the hole densities 
decreased somewhat to 90°C or so, then apparently 
increased at about 110°C. This maximum has been 
noted in practically all of our samples and has been 
previously noted in reference 3. The small range of 
temperatures for thermovoltage measurements possible 
with samples crystallized at 120°C makes a definite 
statement difficult. The data for samples crystallized 
at 210°C shows a continual decrease in hole densities in 
all samples as shown in Fig. 16. From the figure, it is 
seen that the hole densities increase slowly in the region 
to about 85°C then increase rapidly. The data of samples 
(qg) to (0) do not fit the general behavior of the other 
samples. Our previous data on oxygen containing se- 
lenium also indicate a somewhat slower rise in hole 
density than is indicated by the general behavior in 
Fig. 16. Therefore, it is not possible without more 
data to carry the consideration further. The behavior 
can be summarized by saying that hole densities de- 
crease slowly to about 85°C and then decrease very 
rapidly. The spread of hole densities with different 
purity samples is considerably less at high tempera- 
tures than at room temperature. 


4. General Behavior of Hole Mobility-Magnitude 
and Temperature Dependence as Functions 
of Extent of Crystallization and Impurity 


It was shown in the above discussion that the acceptor 
densities remained constant during the process of 
crystallization. In the range of temperatures possible 
with the 120°C crystallized samples, the value of n 
increased but slightly. Thus with these samples, the 
resistivity and temperature dependence of resistivity 
reflect directly the change of mobility and its tempera- 
ture dependence with extent of crystallization. The 
curves for samples (m) and (/) are comparable. We 
notice that b increases from 0.01 cm?/volt sec at 40°C 
to about 0.06 cm?/volt sec at about 100°C for sample 


TABLE V. Efficiency of acceptor level production. 








(3) 
Acceptor levels 





(1) (2) assuming (3) 
Impurity Impurity added exhaustion (2) 
O (0) 1.24 1019 4.3 105 3.5X 1074 
2 (t) 1.06 107° 4.1 10'* 3.9K 1074 
Impurity 
I (q) 2.08 X 1018 7X10" 3.4X 10-3 
2 (v) 2.54X 1019 1.1 10'* 4.3 10-4 
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(1), 0.034 cm?/volt sec to 0.09 cm?/volt sec for (m) in the 
same range. About the same temperature dependence 
is noted for all the low temperature crystallized samples. 
The magnitude of the mobilities varies from 0.003 cm?/ 
volt sec to 0.024 cm?/volt sec at 40°C with impurities of 
oxygen, sulfur, and iodine. Reference (3) has reported 
mobilities ranging from 0.02 cm?/volt sec to 0.001 cm?/ 
volt sec with oxygen impurity. It should be noted that 
oxygen and sulfur decreased the mobilities while iodine 
increased them slightly. 

The mobilities in samples crystallized at 210°C varied 
again with purity. The values for vacuum prepared 
samples and those containing sulfur averaged about 
0.03 cm?/volt sec at 40°C, a factor of three larger than 
those of the corresponding 110°C crystallized speci- 
mens. The addition of oxygen increased the mobilities 
to 0.1 cm?/volt sec while iodine additions further raised 
the value as high as 0.7 cm?/volt sec at 40°C. 

In all cases the mobilities increased very rapidly with 
temperature but in the range below 100°C more slowly 
than those in samples crystallized at lower tempera- 
tures. The temperature dependences of mobilities in 
selenium with sulfur impurities were practically identical 
with those of pure specimens. Impurities of oxygen and 
iodine produced samples with higher mobilities and cor- 
responding lower rates of mobility change with tem- 
perature. 

It is to be noted that with higher impurities of I, the 
mobility approaches 1. However, in all cases there re- 
mained a rather rapid increase of mobility with tem- 
perature. This observation, together with our data on 
high frequency dependence of resistivity gives support 
to our previous conclusion that the observed dependence 
of resistivities on temperature arise from the presence 
of inhomogeneities in the samples. 


5. Nature of the Inhomogeneities in Selenium and 
Their Effective Resistance 


As stated before, it has been determined that these 
inhomogeneities are not of macroscopic order. The 
microcrystalline nature of the material gives rise to 
disordered regions at the grain boundaries. In addition, 
it seems reasonable that the crystallite surfaces have 
higher acceptor densities than the bulk of the crystallite 
previously discussed. 

However, reference 2 has demonstrated the fact 
that the inhomogeneities we are concerned with may 
not be associated with any obvious geometrical dis- 
continuities. The decrease of m and the increase of b 
were observed in single crystals. In the crystals we 
were not able to detect any essential irregularities. 

There are two points of interest about the frequency 
dependence of the resistivities. 


1. All samples decreased in resistivity with increasing fre- 
quency at a rate depending in a regular manner on the magnitude 
of the resistivity. 

2. The frequency dependence of resistivity was independent of 
the “activation energy” associated with the samples. Thus, speci- 
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mens crystallized at different temperatures but having the same 
de resistivity had the same frequency dependence to 200 mc. 


The microcrystalline aggregate may be represented 
by an equivalent circuit (Fig. 18) where R2 represents a 
linear resistance of a disordered region; R, a resistance 
arising from difference in work function; Rs; a dielectric 
loss resistance ; C,, and C2, the shunting capacities 
the resistivity of a perfect lattice of selenium. 

The first attempt at simplification would consist jn 
assuming R,=0, R;=«. With such a model and a 
fixed value of C2 the total dc resistance is given by 
pr=pSe+cR2, where c is the number of irregularities 
per unit length, the ac resistance by 


> Pg. 


cR2 cR. 


pT =pset+ ~ 
wC2R2+ 1 wC?Ri+ 1 





The observed frequency dependence does not vary 
in the prescribed way with frequency, however, unless 
a distribution of values of C2 and R» is taken, which 
procedure may be somewhat arbitrary. 

Assuming then that R2<.R; at low w and R3= Rw" 
(n positive), at higher frequencies: 


pT = pSe+cRow"/(1+ Ro’w?" #C,”), 


By appropriate choice of » better agreement with the 
data can be obtained. Since values of m are generally 
between 1 and 2, the degree of arbitrariness is somewhat 
reduced. 

With difference in work functions between. the two 
regions potential barriers can exist. These are repre- 
sented by R;, C;. There is little point in expanding the 
treatment because of the already large number of 
undetermined parameters. 

The general character of the disordered regions should 
be reviewed. 


1. They are of such a nature as to give a rather slow decrease 
in effective resistance with frequency. 

2. Their resistance and temperature dependence of resistance 
decrease continuously with extent of crystallization. 

3. Non metal impurities in incompletely crystallized samples 
have little effect on their resistance and temperature dependence 
of resistance. 

4. Their resistances and temperature dependences are greatly 
reduced in fully crystallized samples with certain nonmetal 
impurities. 

5. In crystals the dependence of resistivity on field follows the 
behavior of a field reduced potential barrier at sufficiently high 
fields. 


The above behavior is characteristic of that of po- 
tential barriers arising between materials of different 
work functions. A distinct possibility is the adjunction 
of amorphous selenium and hexagonal crystalline se- 
lenium. This view is supported by the consideration of 
the crystallization process and by the fact that X-ray 
analysis has indicated the presence of considerable 
amorphous selenium in microcrystalline samples. 
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The temporal changes of resistance of microcrystalline 
selenium during testing have been noted by a large 
number of investigators. The concurrent changes of 
thermoelectric power and resistivity have been noted 
in Fig. 2. The data showed that the carrier density de- 
creased with temperature and changed somewhat with 
period of time at high temperature but the hysteresis 
effects at lower temperatures were occasioned by 
changes in mobility. The temporal changes of resistivity 
have been noted by other investigators, at low tempera- 
tures. Henisch* has recently considered the presence of 
amorphous selenium as a possible explanation of the 
apparent high mobility of ions at very low tempera- 
tures. The temporal changes in single crystals have 
been investigated in some detail in Fig. 2. 

In short, the assumption of barriers arising between 
disordered regions of amorphous selenium and crystal- 
line selenium seems natural and in conformity with 
existing data. However, it should be noted again that 
the presence of distinct geometric discontinuities in 
single crystals have not been noted so that the addi- 
tional possibility of barriers of the same nature within 
the crystallites has not been eliminated. 


SUMMARY | 


(1) The resistivity of microcrystalline selenium pre- 
pared by crystallizing the black amorphous form is 
only slightly influenced by the initial temperature of 
the liquid from which the amorphous material is formed 
on rapid cooling. This influence was evident only when 
the initial liquid temperature was above 500°C. 

(2) Selenium nucleates below 210°C when the liquid 
selenium is rapidly cooled to room temperature with a 
speed depending on the presence of nonmetal impurities. 
Oxygen and the halogens accelerate nucleation, phos- 
phorus hinders it. 

(3) The resistivity of microcrystalline selenium de- 
pends to some extent on the temperature of nucleation. 
The resistivity also depends on the period and tempera- 
ture of crystallization as summarized in 4. Samples 
nucleated at the temperature of subsequent crystalliza- 
tion decrease in resistivity on crystallization to a mini- 
mum value which is only slightly exceeded on subse- 
quent periods of crystallization. With nucleation at one 
temperature and crystallization at another the re- 
sistivity rises quite rapidly beyond the minimum, the 
rate depending on the amount of nonmetal impurities. 

(4) The extent of final crystallization depends on 
the temperature of crystallization, increasing with tem- 
perature to 210°C. The hole densities remain constant 
during crystallization to the minimum mentioned in 3 
but the effective mobilities decrease by large factors. 

(5) The thermoelectric powers increase with tempera- 
ture of measurement, more sharply at higher tempera- 
tures. The values are very high for pure samples. 





*H. K. Henisch and J. Ewels, Proc. Roy. Soc. (London) 63, 
861 (1950) 
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(6) The hole densities calculated from the thermo- 
electric power decrease with increasing temperature of 
measurement above room temperature. At room tem- 
perature, the hole densities are almost constant with 
temperature (between 10' and 10'*) but the densities 
increased quite sharply as the melting point is ap- 
proached. This behavior may indicate that the number 
of acceptor levels decreases with a perfection of the 
lattice at higher temperatures. 

(7) The ‘‘activation energies” decrease with the ex- 
tent of crystallization to the minimum mentioned in 3 
above. Together with 6 this indicates an exponential 
increase of effective mobility with temperature and a 
decrease in the exponent with the extent of crys- 
tallization. 

(8) A large difference between the values of hole 
densities observed in single crystals and in the micro- 
crystalline matrix of these crystals indicates: 


(a) selenium crystals purify themselves on crystallization so 
that the impurities are probably present in surface layers; and/or 

(b) the geometric discontinuities in microcrystalline selenium 
lead to increased acceptor densities in surface layers. 


(9) In completely crystallized selenium, sulfur in 
amounts to 0.23 atomic percent does not affect the 
electrical properties of selenium. Oxygen and iodine 
increase the hole densities somewhat but at a very low 
efficiency of acceptor production. 

(10) The mobilities are increased and temperature 
rates of increase of mobilities are decreased with in- 
creasing amounts of nonmetal impurities of oxygen and 
iodine so that with the factors mentioned in 6, positive 
temperature coefficients of resistivity are obtained with 
impurities in the order of 0.1 atomic percent. 

(11) The range of effective mobilities varied with the 
various factors mentioned above over a range from 
about 0.003 cm?/volt sec to about 0.7 cm*/volt sec at 
40°C, and up to 8 cm?/volt sec at 190°C (sample with 
0.07 atomic percent iodine). 

(12) The adjunction of amorphous selenium and 
crystalline selenium has been suggested as producing 
potential barriers similar to those arising at the junc- 
tion of semiconductors with different work functions. 
The grain boundaries appear to be the natural location 
for such barriers but in view of the work on single 
crystals the possibility that similar barriers occur 
within the crystallites has not been eliminated. The 
frequency dependence of the resistivity of selenium 
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shows that the assumptions of linear resistances or di- 
electric loss between the crystallites do not explain the 
behavior unless rather arbitrary distributions of these 
parameters are taken. 
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Approximate closed form analytical expressions are derived for the total and differential scattering cross 
sections of cylindrical, prolate spheroidal, and disk-shaped scatterers which subject scalar plane waves to 


only a small phase shift. 


When a plane wave excites a sphere or infinite cylinder, the internal field of these objects can be expressed 
as an infinite series of bessel functions. By making an approximation in coefficient of the mth bessel function 
which is valid for either 27a/\o>>n (a=sphere or cylinder radius, \o>= wavelength of incident plane wave) 
or |m—1|<K1 (m=o/d1, \1= wavelength of plane wave in medium composed of material of scatterer), it 
is possible to obtain a closed form expression for the series. The scattered field is obtained by integrating 


over these internal fields. 


The scattered fields for long but finite cylinders and prolate spheroids are calculated by approximating 
their internal fields by those of infinite cylinders of the same radius. The scattered fields of thin disks are 
obtained by assuming the internal fields to be that of an infinite flat plate of the same thickness. 


1. INTRODUCTION 


HERE is considerable interest in the scattering 

of electromagnetic and sound waves by obstacles 
which subject an incident plane wave to only a small 
phase shift. In many important cases (for example in 
the scattering of light by virus proteins) the ratio of 
the propagation constant k,;=22/, in the scatterer 
(A1= wavelength of plane wave in scatterer) to that of 
its surrounding medium, kyp=22/Xo, is between 1.1 
and 1.4. Scattering measurements under these condi- 
tions are usually interpreted on the basis of the Ray- 
leigh'-Gans*-Born’ theory. In this theory one computes 
the wave function of the scattered wave by assuming 
that the obstacle is an assembly of infinitesimal scat- 
terers, each of which is excited by the field of the in- 
cident wave. The field associated with the scattered 
wave is then taken to be the linear combination of the 
scattered fields of each of the scattering elements of 
the obstacle. Of course, each infinitesimal element is 
excited by the internal field associated with its point 


* Most of this work was done while one of the authors (EWM) 
was a visiting professor at the Institute for Mathematics and 
Mechanics at New York University. It was partially supported 
by a contract between the Institute and the ONR and partially 
by the Bureau of Ordinance contract NORD 7386 with the 
Applied Physics Laboratory. 

! Lord Rayleigh, Proc. Roy. Soc. (London) 84, A25 (1911); 90, 
219 (1914). 

2 R. Gans, Ann. Physik 76, 29 (1925). 

3M. Born, Z. Physik 37, 863 (1926) and 38, 803 (1926). 


of location rather than by the incident field which 
would have existed there in the absence of the scatterer. 
The limitations of the Rayleigh-Gans-Born approxi- 
mation are well known and will not be discussed here. 

In the first paper* of this series an improvement over 
the R.G.B. theory was made for scattering by spheres 
for which k,/ko<1.5. The exact solution of this problem 
can be expressed as a series*:* of products of bessel 
functions and Legendre polynomials. 

The coefficients of these products are complicated 
combinations of bessel functions of the parameters 
ak, and ako (a=sphere radius). By making an approxi- 
mation in the denominators of these coefficients which 
is valid when ak,, ako>n (n is the ordinal number of 
the term coefficient being approximated) or when 
a(k,—ko)<1, it is possible to sum the series for both 
the internal and scattered wave functions by using the 
addition theorems for bessel functions. Although the 
resulting expressions are valid only when k,/ko<1.5, 
they have the advantage of being in a simple analytical 
form. They also furnish usable expressions in a range 
where the R.G.B. theory is no longer valid but where 
many experiments are conducted. 

The purpose of this paper is to generalize the results 


*R. Hart and E. Montroll, J. Appl. Phys. 22, 376 (1951). 
This paper will be referred to as S-I. 

5Lord Rayleigh, Theory of Sound (Dover Publications, New 
York, 1945), reprint edition, Vol. 2, p. 272. 

6G. Mie, Ann. Physik 25, 37 (1908). 
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SCATTERING OF PLANE WAVES BY SOFT OBSTACLES II. 


of S-I to finite and infinite circular cylinders, spheroids, 
and flat plates. The approximate internal and scat- 
tered wave functions for infinite circular cylinders are 
obtained by making the same approximations in the 
coefficients of the exact bessel function series that were 
made in the sphere problem. As in the sphere case, 
these approximations allow one to sum the required 
series to obtain simple analytical expressions for the 
wave functions. 

The scattered field and total scattering cross section 
of prolate spheroids and long finite circular cylinders 
will be calculated by assuming that the internal fields 
in these scatterers is the same as that which would 
exist in an infinite cylinder-with the same propagation 
constant and diameter (minor axis in the case of 
spheroid). Then the technique of R.G.B. will be im- 
proved by assuming that each infinitesimal element of 
the scatterer is excited by a field that is approximated 
by the internal field of the infinite cylinder. The scat- 
tered field of the obstacle will be expressed as the sum 
of the scattered fields generated by all the infinitesimal 
elements. 

The scattered fields of finite plates will be calculated 
by using the internal field of an infinite flat plate of 
finite thickness. 

The scalar wave function in a medium of variable 
propagation function is 


Vy+ k(r)y=0. (1) 


We shall be interested in situations in which a homo- 
geneous isotropic scatterer with well-defined boundaries 
is immersed in an otherwise uniform medium. We shall 
restrict ourselves to y’s which satisfy the boundary 
conditions (a) y and its first derivatives continuous at 
the boundary of the obstacle and (b) y and its first 
derivatives bounded at infinity. These correspond to 
(a) the quantum-mechanical case, (b) the scalar analog 
of nonconducting dielectrics in the electromagnetic 
case, and (c) scatterers whose density is the same as 
that of the surrounding medium in the acoustical case. 
The methods used here can be extended to the general 
acoustical situation provided that of the scatterer does 
not differ much from that of the surrounding medium. 
In this paper we shall be concerned with scalar waves. 

Let V represent the volume occupied by the scat- 
terer. Then we postulate 


k(r) = ky 
k(r) =k, 
Also, let so be a unit vector in the direction of propaga- 
tion of the incident plane wave and ~;=expiko(r- So) 


its wave function in the absence of a scatterer. If y,(r) 
is the wave function of the scattered field, 


¥(r)=y.(r)+expiko(r- So). (2) 
Then (1) can be rewritten as 


Vet ko s= Lhe? — k(n) Wy. 


if rv, 
if rcv. 
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It is known that the differential equation 
VY +hoy.= F(x, y; z) 


with our boundary conditions has the solution’ 





1 etko|R—-t| 
y.(R) - -—f F(x, y; z)dr, 
4nJ |R—r! 


where dr is the volume element of our three-dimen- 
sional space of interest and the integration extends over 
all space. This equation implies that (1) is equivalent 
to the integral equation 





1 etkolR—r| 
8 R)=-— ke —k dr. 3 
VAR)=——— J The OMT. Ge 


Now let RCV. Then ko?—?(R)=0. Hence, if we let 
¥(r)=¥(r) when rCV, we have 





ke ke e*ko|R—t| 
y.(R)= f p.(r)dr. (3b) 
4dr V | R— r| 


Frequently one requires ¥,(R) only when R is much 
larger than any dimension of the scatterer. Then 


r-s 
|R—r| =[R'+-r°—2Re-s,}=R| 1 —+0(R) | 


where s; is a unit vector in the direction of the line 
connecting the center of the scatterer and the point R. 
Hence for large R 


(ki? — ko”) 
¥(R)~— ein fF exp —(r-si)iko Wi(r)dr. (4) 
4rR Vv 


This equation is completely rigorous. If y,(r) is re- 
placed by the incident wave function, one obtains the 
Rayleigh-Gans-Born approximation. We shall use the 
approximate internal wave functions for spheres, in- 
finite circular cylinders, and infinite flat plates to derive 
the scattered wave functions ¥,(R) for various shaped 
obstacles. 

As an example of the method of application of (4) 
let us find the approximate scattered field of a sphere 
of radius a for which 1<,/ko<1.5. Under these condi- 
tions and the boundary conditions ¥;(a)+y,.(a)=y,(a) 
and {d/drLy.(r)+nrp.(r) }} raa=((0/0r)(rp(r)) |,—a, it 


was shown in S-I that 
¥(r) 

2k e'**1-*0)[ exp(ikir- So) — x exp(2iak)exp(— ikir-So) | 
7 (kitko)[1—« exp4iak, | 





’ 


(5) 
where x= (ki—ko)/(ki +k). When y¥;(r) is substituted 


7 Mott and Massey, The Theory of Atomic Collisions (Oxford 
University Press, London, 1933), p. 85. 
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into (4), the resulting expressions involve the integrals 


i= f expfir- [his bess dr (6a) 
V 


and 


h= f exp{ —ir-[ki80+ kos; |} dr. (6b) 
Vv 


To integrate the first we let the polar angle 7 be the 
angle between the fixed vector [k18:—os;] and the 
variable vector r which ranges through the entire 
sphere. Let | ki8o— kos:| = [2+ o?— 2kiko cosd |! =a, 
where @ is the angle between 8» and s;. Then the first 
integral is 


h=2r f f exp(irw; cosy)r? sinydydr 
0 0 


4r p™ 
= —f u sinudu 
w,* 0 








— (*). (a) (7a) 
= ~ aw), a 
(aw, 2) 0" 
where J;(x) is the bessel function of order $. Similarly, 
we have 
, 4a’ (*). (av,) 
2= = av), 
(av,)! 2 ~~ 
where 


Y= [k+ ket+ 2kiko cosé }}. (7b) 


By substituting (5) into (4) and employing (7a) and 
(7b), we obtain the following expression for the differen- 
tial scattering cross section, o,(@)dQ (the fraction of 
incident energy scattered into the solid angle dQ at 6): 


o,(6)dQ= R*|y,|2d2 


JP (dw. 
2 cos(2k,a) 


A,|? 





a) 
Jy(ae1)Fy(a01) J¥(an) 


. (av,)° 





» (8) 


(a’w 0)! 
A i2 2rk;?(ki— ko)?a® 
wip [1+«'— 2x? cos4ak,] 





This is exactly the result obtained in S-I by a somewhat 
different method. 


2. INTERNAL FIELD FOR INFINITE CYLINDER 


Let a plane scalar wave of wavelength Xo be scattered 
by an isotropic, infinite circular cylinder of radius a 
whose axis is coincident with the z axis of our coordinate 
system. Let the incident wave be propagated in the 
direction of a vector 89 which lies in the x—z plane and 
which makes an angle a» with the z axis of our coordi- 
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nate system. The wave functions of the incident, scat. 
tered, and transmitted waves satisfy the wave equations 


Vy+kry=0, r<a, ky=2n/r (inside), (9a) 
V¥+ky=0, r>a, ko=2x/do (outside). (9b) 


In the acoustical case y is the time-independent factor 
of the pressure wave of frequency v, p=y exp(— 2ziv!), 
At the boundary of the cylinder the pressure and radial 
component of the particle velocity must be continuous; 
that is, the quantities y and (—i/pc)[dy/d(kr) ] must 
be continuous. We represent the densities of the scat- 
terer and medium by p; and po, and the corresponding 
sound velocities by c; and co. Furthermore, we define 
the ratios m and g by 


m= k,/ko=co/c1 (10) 


The case g=1 corresponds to the quantum-mechanical 
boundary condition at r=a and leads to considerably 
simpler results. 

The general solutions of the wave equation in 
cylindrical coordinates are® 


Yn= Zar — Hi? Weintsioe, (11) 


where Z, represents either the mth bessel function, J, 
or the mth hankel function H,“. The incident, trans- 
mitted, and scattered wave functions can all be repre- 
sented as linear combinations of these functions. 

The incident wave is given by 


and g=pi/po. 


y= exp(ikor : So) = erkoz cosadpirkg sinag cos@ 
ee) 
= eikoz cosag ym ie” 9J (ro sina). 
n=—0 


(12a) 


We choose the following expansions for the transmitted 
and scattered waves: 


Vi= CX Aa ni", (rho sinag)eir*+itoz cosa, (12h) 


n=—oO 


es) 
w= ye A: ni" J (rk, sina, ei oti 2 cosa) 
n=—o 


(12c) 


where {A,.n}, {Aen}, and a are to be chosen so that 
the required boundary conditions are satisfied. 

In order for ¥;+y¥,.=y; at r=a for all z, a; must be 
related to ao, ko, and k; by 


Ro COSap= ky Cosa, (13) 
and A, , must be related to A,, through 

At, nJ n(aky Sina) — As, nH x (ako sinay) = J n(ako sinay). 
The continuity of (—i/pc)dp/d(kr) at r=a implies 


Copo sina 1 


At nJ n (ak; sina) 
C1p1 SiNa ; 
—A,, nH n' (ako sina) =J n' (ako sina). 


8 J. A. Stratton, Electromagnetic Theory (McGraw-Hill Book 
Company, Inc., New York, 1941). 
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SCATTERING OF PLANE 


We use the abbreviations 





k;*=k;sina; and m*=k,*/ko* (14) 
and the identity 
J n(x) © (x) —JFa(x)H n(x) =2/rix (15) 
to obtain 
2g/(miako*) 
Aun 


~ m*H, (ako*)J »’(ak:*)—gJ a(ak:*)[H a (ako*) 





—m*J ,'(aky*)J ,(ako*)+ gJ n(aki*)J » (ako*) 
Agn= , 
m*H,, (ako*)J n' (aki*) — gJ n(aki*) [Hn (ako*) |’ 


In order to sum the series (12b) and (12c) in closed 
form, we examine several approximations of the de- 
nominators of A,,, and A;,». A rather simple approxi- 
mation can be obtained when ak,* and ako* are much 
larger than x. When x>>n, 


H,,(x)~(2/rx)' expi[x—3(2+1)r], (15a) 


J n(x)~(2/xx)! coslx—}4(24+1)r]. (15b) 


Since hg = 1(Z n—1(X) — pa 1(x) 1, the 


nominator is approximately 


i(m*+ g)exp[ia(ko*— k,*) ] 
ma(k;*ko*)! 


required de- 





~~ 
r— 


X[1+ix«(—1)" exp(2iak,*)], (16) 
where ako*, ak,*>>n, and 


m*— g 


= wae 





K 


It is easy to show that this approximate expression 
approaches the correct one asymptotically for all » as 
g—1 and m*—1 (hence as k,*—>k,* and x—>0), for then 
the denominator D, approaches 


H, (ako*)J n'(ako*)—LHn™ (ako*) J n(ako*) 
= —2i/(nako*), 


which is exactly the limiting form of (16) under the 
same conditions. 


The approximate expression for the coefficient A: n 
is then 


2g(m*)} 
(g+m*) 
nen? — ko*) |[ 1—ix(—1)* exp(2iak,*) ] 
[1+ exp(4iak,*) ] 


Ata 
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If this is substituted into (12b), we obtain the following 
approximation to the internal wave function: 


2g(m*)} exp[ia(ki* — ko*)+-izko cosag | 
(g+m*)[1+  exp(4iak,*) ] 


~~ 





x > ine J .(rky*)[1 _ ix(— 1) "eon", 


n=—D 


But (12a) and 


e~irki*® ccs? as >” in(— 1)"J ,(rky*)e"® 


—oO 


imply that 


2g(m*)! exp[ia(ki*— ko*)+izko cosay | 


““(gtm*) [1-442 exp(4iaks*)] 





(18) 


X {eitt*r 089 ix exp(2iak,*)exp(—irk,* cos0)}. 


The scattered wave function y, is derived in Appendix 1 
by summing (12b). 


3. SCATTERED FIELD—CIRCULAR CYLINDERS, 
INFINITE AND FINITE 


When g=1, the scattered field for an infinite circular 
cylinder can be obtained by substituting (18) into (3b) 
and integrating over the entire cylinder. Then 


2(ky2— ko?) (m*)} exp[ia(ky*— ko*) | 
¥(R)= a 
4n(1-+ m*)[ 1+ R? exp(4iak,*) ] 





XK {Ts— inl e?2™"}, (19) 


where 


n= f f f exp(izko cosao)exp(ik,*r cos@) 
0 0 _ 


Xexp(iko| R—r|)| R—r|—rdrdédz, 


=f f f exp(izko cosao)exp(—ik,*r cosé) 
0 0 —n 


_ Kexp(iko| R—r|)| R—r|—'rdrd6dz, 


c= (m*—1)/(m*+1), |R-#|=[(Z—2)+ 8}, 


s?= R°+7r?—2Rr cos(@— ©). 


The coordinates of the interior points are (in cylindrical 
coordinates) (r, 8, Z) and those of the external ones are 
(R, 0, Z). The integration over z in J; and J, leads to 
a hankel function. It is performed by letting z—Z 
= § CSCao(sinht— cosap cosh/). Then 


[(z—Z)*+s? }!=s cscao(cosht— cosa sinht) = dz/dt, 
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so that® the integral over z becomes If the direction of incident wave is normal to the se 

2 il wee bol (g—-Z)2-+s? Bids cylinder axis (ao=}7, ki*=ki, ko*=ko, and m*=m), 

J pa 0 cosa)exp| ikoL ( a EP+ #3 ne then one has ae 

- [(s—Z)?+s*]}} 

- e 2ratk,(m— 1)*ko? J (awe) 
=exp(ikoz cosao) f exp(isko sina cost)dt ool ~2 x costakit«*]\ (awe)? 
= mi exp(ikoz cosa) (sko*). , J (aw2)J (ave)  J,?(ave) 
: : +2« sin2ak,; K (22) 

If the scattered wave is observed at a distance of awW2Vo (av2)? 1 
many cylinder diameters from the cylinder axis so that ™ 
sko*>>1 and R>r, then Notice that this is very similar to the corresponding 
s= RI1—(2r/R)cos(@— 0 / py result for spheres (Eq. (8)) except that the bessel func. 
pun BLA — (27 /Rconh al v RY! ‘ tions of order $ have been replaced by those of order one. 

~R[1—(r/R)cos(@— ©O)+O0(r/R)* }. The scattered field at large distances from a long but Ag 
Hence (15a) implies that the integral over z is, as finite circular cylinder of length L can be obtained ina pla 
Dine similar manner by substituting (18) (with g=1) into 
; (4) and integrating over the entire cylinder. Since P 
wil 9 (sko*)~i(2x/Rko*)exp(—im/4) r-8,=z cosh+r sin® cos(@—@) in cylindrical coordi | %4 
Xexp(iRko*)exp[ —irko* cos(6— @)]. nates (we still choose (z, 7, 6) as internal coordinates, 
Then but we use spherical coordinates to define the point of 
Is~i(2e/Rko*)! exp(— }wit+-iRko*+ikoZ cosas) ey eee 
or 4G a : onwee tt 2(ki?— ko?) (m*)* exp[ia(ki* — ko*) Jexp(ikoR) 
x exp! ir(k,* cos@— ky* cos(@— ©)) |rdrdé AK) . 
J J on tetalaaitt aie )”) 4 (1-+-m*)R[ 1+ exp4iak,* | : 
I y~i(2x/Rko*)* exp(— jrit+-iRko*+ikoZ cosay) X {Is—ix exp(2iak:*)Ie}, (23) by | 
2r na where 
x f f exp[ —ir(k:* cosé+ ko* cos(@— O)) |rdrdé. 
0 0 @ pnltr 4L 
, by § 
[= k * 6 + a 
But Eq. (A-II-3) in appendix II implies that J i in mee ont; = 
sult 
~ ; / * } / 9 9 . . 
Ty~i(2x/Rko*)"(2ra, w2) J s(awv2) Xexp[—ikor sin? cos(@— @) ] sphe 
Xexpi(Rko*—4a+koZ cosao), (20a) infir 
Iy~i(2e/Rko*)*(2xa/v2)J (avs) Xexp[izko(cosa»—cos®) |rdrdédz, (23a) | imp 
Xexpi(Rko*—4a+koZ cosao), (20b) @ ptt all T 
where I,= f J J exp(—irk,* cos@) 
w= ky*?+ ky*?— 2ko*k;* cosO, (20c) 0o%o Y-4L v.(R 
vor = ky*?+ ko*?+- 2ko*k,* cosO. (20d) Xexp[—ikor sin® cos(@— @) ] 

By substituting (20a) and (20b) into (19), we obtain Xexp[izko(cosa>—cos®) rdrdddz. (23b) 
the approximate differential scattering cross section che 
o,(Q)dO=R|y,(R)|2dO Equation (A-II-3) of Appendix II implies that are ¢ 

and 
J ?(aw2) P 2ra*L sin{4Lko(cosay—cos®) I 
= |B, |?) ——+ 2« sin2ak,* Iya ad. il ia 
we aw3{4LRo(cosay—cos®) } ies 
4 Ne 2 Jo 
ae) own lors) a (av2) dO, (21) 2ma?L sin{4Lko(cosay—cos®)} 
Wee Ve" I,= J \(a0s), (24b) 
where av3{ }Lko(cosas—cos®)} 
2ar(ki?— ho?)2a2ki*/(Ro*)? where ; 
| B,|?= ' we=k,+ ke? sin’>—2kok,* sin® cosO, (254) | wher 
(1-+-m*)*[ 1+ «4+ 2x? cos4a,k,* | | 
EE v3°= ky+ ko? sin?@+ 2koki* sin® cos. (25b) 

*G. N. Watson, A Treatise on the Theory of Bessel Functions 
Co Pe ee om iene: If we substitute (24) into (23), we obtain the fll} Tie 
we shall make references to it in the form Watson, p. 180, Eq. (10). ing expression for the differential scattering cros 
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section : 
o d2=R*|y,|*d2 





,sin*| 3 Lko(cosay—cos®) } {- 1°(aw3) 
=|c 
{3 Lko(cosap—cos®)}? | (aws)? 
J \(aw3)J,(av3)sin2ak;* J,?(av;) 
+2k +< , (26) 


a7w303 (av3)? 





where 
Lat(ke— ko?)?m* 


|c,|*= > 
(1+ m*)?[1+- 2x? cos4aki*+ x*] 





Again, in the special case ap=2/2 (normal incidence of 
plane wave) 





od L?atkyko?(1—m)? [= cos®) } {- 1°(aw3) 


6 dQ= 
(1+ 2x? cos4ak,+«*)| 3Lko cos® (aw3)? 





2x sin2ak,J;(aw3)J;(av3)  J,?(avs) 
+ +6}. 
a7w 303 (av3)" 
4. SCATTERING BY PROLATE SPHEROIDS 


Let us now compute the scattering of a plane wave 
by the prolate spheroid 


[(at+y%)/a}+-2/0=1 


by substituting the internal field of an infinite circular 
cylinder of radius a into (4) and integrating the re- 
sulting expression over the space occupied by the 
spheroid. As b/a increases, the spheroid approaches an 
infinite cylinder; hence the internal field approximation 
improves. 

The scattered wave function is 


. 2(k:?— ko?) (m*)} exp[ia(ki*— ko*) lexp(iRko) 
_ 4nr(1+m*)R(1+ % exp4iak,*) 
X(T1—ixls exp(2iak,*) |, (28) 





where J; and J, are integrals with integrands which 
are exactly the same as those in J; and J¢, Eqs. (23a) 
and (23b). The integration is taken over the spheroid. 

It is convenient to integrate over the coordinates 
{,, 6 defined by 


cosf=cosg/[1+¢ sin’¢ ]}, 
sinf= (1+)! sing/[1+¢ sin’¢ }}, 


n= (r/b)[1+2 sin’ }}, (29a) 
where (r, 8, ¢) are the usual spherical coordinates and 
c= (b’—a?*)/a?. (29b) 


The surfaces »=constant form a family of spheroids 
of which »= 1 corresponds to our scatterer. The volume 


‘element in this coordinate system is a*bn* sin¢dnd¢dé. 


The limits of integration are OC n<1, O< 0S 2x and 
O<¢<-7. In the new coordinates 


2a T 1 
I;= bat f f f exp(ik,*an sing cos@) 
0 0 0 


Xexp[ —ikoan sing sin? cos(6—®) | 


Xexp[ibnko cosf(cosay— cos) }y? sintdndtdé, 


2r 7 1 
Is=a’b f f f exp(—ik,*an sing cos6) 
0 0 0 


Xexp[—ikoan sing sin® cos(@—®) | 
XexpLibnko cosf(cosae—cos®) |y? singdndgdé. 


Integration over @ and ¢ has been carried out in Ap- 
pendix II. Equation (AII-6) implies that 


(22)%a’b 





f (nos) (moos) 


7 = 
we 
= a?b(24/w4)*J; (ws), 


I3= a*b(2m/v4)?J3(v4), 
where 


w4= a?k,**+- a2hk,? sin’-+ ko?b?(cosay— cos®)” 
—2a*k,;*ky sin® cos, (30a) 
¥s= 07h, + ah? sin’-+ k,?b?(cosap— cos)? 
+2a*k,*ky sin® cos@. (30b) 
Hence the scattered wave function at R, 0, ® is 
2(ki?— ko?) (m*)4a?b(27)! 
4a (1+m*)R[1+ 4 exp4iak,* ] 





v.(R, 0, &)~ 


X [ws Jy (ws) — ik exp(2iak;*)vg74J; (v4) |, (31) 


and the differential scattering cross section is a, 
= R?|y,|? or 


22(ky?— ko?) m*a‘b? 
(1+ m*)*[ 1+ x*+ 2x? cos4ak,* | 





X | cory (cos) —ixeiae*y AT; (04) |2. (32) 


5. SCATTERING BY CIRCULAR DISKS 


To obtain the scattered wave function of thin disks 
we shall use (4) and approximate y; by. the internal 
wave function of an infinite flat plate. Let us first find 
the internal wave function of the flat plate. We assume 
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the direction of the incident plane wave to be in the: 
x—z plane and at an angle a» with the z axis. We also 
assume the middle of the plate to lie on the x— y plane. 
In the absence of the plate the wave function of the 
incident wave is ~y=expiko(z cosao+ x sinao). Hence 
after insertion of the plate of thickness 26, the wave 
functions in the regions below, inside, and above the 
plate are, respectively, 


y =exp(ixko sinay)[exp(izko cosao) 


+ B exp(—ikoz cosao) ], 





2ko* exp i6(ki*— ko*) |[exp(izk:*)+ x exp(2i5k1*)exp(—ik*z) Jexp(ixk, sina) 
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y® =exp(ixk; sinay)[C exp(ikiz cosa) 
+D exp(—ikz cosa), 


Y= E exp(ikoz cosay)exp(ixko sina), 


where the constants B, C, D, E, and the angle a, are 
to be determined from the boundary conditions at the 
planes which separate regions of different propagation 
constants. We use the boundary conditions ¥ continuous 
and dy/dz continuous for every x. Then it is easily 
shown that 





y® 


’ (33) 


(ki*+ ko*)[1—«? exp(475k,*) | 


where k;*=k; cosa;, c= (ki*— ko*)/(ki*+Ro*), and ay is 
defined through 


ky sina,= ko sina. 


Now we shall derive an approximate expression for 
the scattered wave function of a thin circular disk of 
radius a and thickness 25. To this end we substitute 
(33) into (4) and integrate over the space occupied by 
the disk. Then as R, the distance of the observer from 
the center of the disk, approaches infinity, 








(k?— het) of Dkyteis(hita" 
4nR ——L(ky*+-ho*)(1—«? exp4iok,*) 


— X (9+ lio), 
b= f f f e'**1* exp(irko cosé sinag) 
—5”0 0 


Xexp{ — kos cos@+r sin® cos(6— O) ]}rdOdrdz, 


C) a 24 
tw=f f f exp(—izk:*) exp(irko cos@ sinao) 
-s%o9 Yo 


Xexp{ —ikolz cosh+r sin? cos(@— O) ]}rdédgdz. 


Here we have used the expression for r-s, given above 
Eq. (23). We use cylindrical coordinates for the internal 
variable and spherical polar coordinates for the position 
of the observer. 

The integration over z in Jy and Jy is elementary. 
The remaining double integrals over r and 6 are of the 
form (AII-1) of Appendix II. We find 





2 sin{ dL ki* — ko cos® }} (= 


“ == 
[k i* —_ Ro cos® | Y 








2 sin{ dL ki*+ ko cos® ]} / 27a 

= ( )r ay), 
[ki*+ ky cos® | v 

where 


v= ke*[sin’ao+sin*— 2 sinay sin® cosO]. (34) 








Hence the differential scattering cross section at an 
angle (®, ©) is given by 


o.(&, O)dQ= |y,|*R’dQ 
4a? (ke — ho?)?(Ro*)*dQ 
=J?(ay)— 
7 (Ri*+ Ro*)? 
sind(ky*— ky cos®) i xe?*®ki* sind(ky*+ ko cos®) 
(ki*— ko cos®) iy ky*+ ky cos® 


1+ x*— 2x? cos(46k*) 





2 














(35) 


In the special case of the incident wave being normal 
to the plane of the disk (a>=0), we have 


4a? J\?(ako sin®)(ki— ko)” 
sin? [1+ «4—2« cos4éhk; | 
xe?*®k1 sind(ky+ ko cos®) 
(ki + ko cos®) 


sind(ky— ko cosé) 
(ki— ko cos) 


2 


a,(&, O)= 














The Rayleigh-Gans-Born approximation to a,(®, 9) 
is obtained by letting &:—p in all terms of (35) except 
the (ki— ko)? term. This leads to 


a,(®, O) 
4a*(k,— ko)?J 2(ay)sin? { 5ko(cosao—cos®)} rn 
~ . f 


y?(cosay— cos®)? 





6. TOTAL SCATTERING CROSS SECTION 


The approximate total scattering cross sections cat 
be obtained by integrating our expressions for «,(6) 
over all solid angles dQ. The formulas for o,(O) in all 
obstacles that we have considered are sums of three 
terms, two of which involve the square of bessel func 
tions. The third term (see, for example, Eq. (22)) com 
tains a product of two bessel functions each of whichis 
a function of a different variable. Since we have not 
succeeded in integrating this mixed term, we shal 
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jimit ourselves to cases in which k;/p is close to unity. 
Then the mixed term is unimportant, and one can ex- 
pect to obtain good results for the total scattering cross 
section without taking it into account. It was shown in 
S-I that in spheres with k,/o as large as 1.5 one ob- 
tains accurate values of ¢, in a large important range of 
x= 2mna/Xo 

Van de Hulst” has pointed out that as m= k,/ky—1, 
the appropriate parameter for the expression of the 
results is y= x(m—1). We shall show that as m—1, the 
total scattering cross section, o,, has a very simple form 
in several of the cases that we have considered in this 
paper. 

(i) Spheres 


The integral over all angles of the differential scatter- 
ing cross section of a sphere of radius a (see Eq. (8)) is, 
as k= (ki— ko) /(Ritko)—0, 


2r 7 
02k? ki— ko)?xra® f J 
0 0 


where w, is given by the equation above (7a). If one in- 
troduces the new variable u=aw, one obtains, after 
the integration over ¢, 


J ?(a1) 





sinOdQd®, 


z(m+1) 
o [429 0?x?(m—1)?/m ] uJ °(u)du. 
z(m—1) 


Now in Watson, p. 417, Eq. (9), it is shown that if 
Re(v)>3, 


” I'(v—3) 
J J(u) u (42 — 2°)"Idu —_ 


z qigrt 





~H.(2z), (38) 


where H,(2z) is the Struve function defined on page 
328 in Watson. The integral in o, corresponds to 


H,(2z). Then 


6 ~2'a?x?(m—1)?m- 


= 1}) Hy(2x[m+1)) 





(39) 
(x{m—1])*? = (x[m+1])>” 
where (Watson, p. 333, Eq. (3)) 
H,(z) = (2/2x)*(1+ 22%) — (2/2z)*(sinz+2— cosz). (40) 


To obtain the limiting expression for ¢, as m—1 we 
let y=x(m—1) and note that x(m+1)=~y(m+1)/ 
(m—1)— for any finite fixed y. Since Hj(2z)—(z/z)! 
as 2-0, we have 

o,/wa>—>2(m/y)*H3(2y) — 2[ (m—1)/(m+1) P 
or as m—1, 
o,/na?= 2(m/y)Hy(2y). (41) 
This is exactly the limiting expression obtained by 
Van de Hulst” and the authors.’ As x (and therefore as 
y) approaches infinity, o,/2a? approaches" 2. 


”H. C. Van de Hulst, Recherches A stronomiques de I’ observatoire 
@Utrecht 1 (1946). 
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(ii) Infinite Circular Cylinders 


To find the total scattering cross section per unit 
length of an infinite circular cylinder when the direction 
of the incident wave is normal to the cylinder axis, we 
integrate (22) over all values of @. Then when m is 
close to 1 (or x to 0) (say, m<1.5), 


oc 2nathythy(m—1)? i€ T2(aus)ea-*d 
wo? = k?+ k— 2Rok cos@. 


An alternative expression can be obtained for oa, if one 
introduces the new variable u=aw,. Then 


o 8k ke? (m— 1)*a* 

7 J?(u)du 
x , 
z(m—1) Uf [u?—x?(m—1)? ][x2(m+1)?—u?]}4 





42) 


where x=akp=27a/Xo. This integral does not seem to 
be expressible as a single function of any tabulated func- 
tions. However, if we let y= «(m—1), Eq. (42) simplifies 
considerably as m—1. For fixed y, x(m—1)=y(m+-1)/ 
(m—1)— and 





8ary’ma f J?(u)du 
(m+1)J, ulw—e} 
Then Eq. (38) implies that as m—1, 
o,/2a=nH,(2y), (43) 
where H,(2y) is the Struve function of 2y of order 1. 
Since H,(2y)-2/r as yoo (ie., as 2ra/Ay>~), 
a,/2a—2 as is required. This result has been stated by 


Van de Hulst." The Struve functions have been tabu- 
lated in Watson, p. 667. 


(iii) Finite Circular Cylinder 
The total scattering cross section of a finite circular 
cylinder can be approximated by integrating (27) over 
all angles and © when the incident wave is propagated 


at right angles to the cylinder axis: hence (neglecting 
terms of order x and x*) we have 


sin(}.LRp cos®) )? 
sin®d®, 
Lo cos® 


2m Ji 2 
ro [ESI 


and w3 is given by (25a): 
= k?+ ko? sin’?— 2kok, sin® cos@. 


It is convenient to introduce the new variable 
u=cos®. Then 


o,= La'ksko}(m—1)? f flu eran “ ba (44a) 


"1H. C. Van de Hulst, Astrophys. J. 112, 1 (1950). 





™= Latkyko?(m— vf Fe 


where 











1286 E. W. MONTROLL 


where 

f(u)= 

af *J P(a[ ko?(1—u?) +k? — 2k ko(1— U2)! nO), 
0 «6 @L Ro? (1—w?) + ky? — 2kiko(1—u*)! cosO] 





(44b) 


When Lp is very large, the main contribution to o, 
comes from values of u in the neighborhood of u=0. 
In this neighborhood f() is a slowly varying function 
of «u while the term in the curly brackets of (44a) 
decreases very rapidly. Hence as Lkyp—>~, the first 
term in the series, 


du 


o.= Lak yko(m— v4 / of { og 


. 1 sin?($ Lkou)du 
+" f +f 
- 421k? 





should give a good approximation to ¢, for all @ when 
k,/ko is close to unity. Since 


f eneoenee 4 (sits 2 =) 
———-} du=—[_ SiLky -———— }, 
il AL Rou Lk i 


' sin?( Luko) 
— u= 


- Lh 


Qa 
Sia- f z~! sinzdz, 
0 


we have (since += ako) 





2 
(Lko—sinLko), 
3,3 


where 


Os 2 sin*4 Lo 
——= m(m— | o siczts) 0 
2La Lko 
1 . 
+—(Lko—sinLky) f’’(0)+--- | 
Lk? 


By direct differentiation of (44b) we find 





ad J2(alke+ ke 2kiko cos@ }!)dO 
fo=2f 
0 a*[k+ ko? — 2kiky cosO | 
z(m+1) 


4 J ?(v)dv 
7 i of v? — 2?(m— 1)? ]§[x2(1-+m)?— ey 


z(m+l)r J 42(y) 9 
roma 








[x?(m?—1)—v* ]dv 





[o*— (m—1)*x? h[a*(1+-m)— 2} 


AND R. 





W. HART 


The integral f(0) is that which appears in the expres. 
sion for the total cross section of an infinite cylinder 
(Eq. (42)). It can be expressed simply in terms of taby. 
lated functions only when m—1: Then if y=x(m—1), 


f0)~L4(m— 1)? "Hi, (2y). (46) 


As m—1, we have the following expression for f’’(0): 


hover, (PROV _& 
ro=»f || 


u 


— f J 1(v)J2(v)dv oa 
Os, eee . 





We can make a rough estimate for the conditions 
under which the f’’(0) term in (45) can be neglected 
when compared with the /(0) term by examining (46) 
and (47). For sufficiently large v, |J1(v)Jo(v)| <2/xn. 
Hence for sufficiently large y we have 


1 
if’ (0)| <-f a : 


When Lf is large, SiLko~~}m; and as y>~, H,(2y) 
~2/n. Therefore, for large Lko and y, the first term in 
the bracket of (45) is ~2/xy*. On the other hand, the 
absolute value of the f’(0) term under these conditions 
is £1/Lkoy*. This means that as long as L>Sa, the 
second term of (45) is less than 10 percent of the first. 
Our initial hypothesis that the internal field of a cyl- 
inder can be approximated by that of an infinite 
cylinder of the same radius is probably not valid when 
L<S5a. On this basis, when L/a becomes sufficiently 
small so that the second term in (43) is important, little 
confidence can be put in (43). In the range of validity 
of (43) we have for the ratio of the scattering cross 
section of a finite cylinder to that of a section of an 
infinite cylinder of the same length 


o.(L)/o.( 0) = (2/m){ SiLko— (Lko)(1—cosLko)}. (48) 


As the relative index of the scatterer to that of its 
surrounding medium approaches unity, we have 


o,/2La=2H,(2y)[Si(Lko) 
—(Lko)“(1—cosLko) } (48a) 





where 2a is the geometrical cross section of the cyl- 
inder and y= 27a(m—1)/Ao. Equation (48) is probably 
correct for all super-soft cylinders with L> Sa. 


(iv) Thin Disks 


The total scattering cross section of a thin circular 
disk can be approximated by integrating (36) over all 
angles (6, @) when the direction of propagation of the 
incident wave is normal to the plane of the disk. If in 
the integration terms of order « and x’ in (36) are neg 
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lected we obtain 





* J?(x sin®) sin*éko(m— cos) 
c= 8e'(m—1)* f ie 


sin® ko?(m—cos®) 





i J °(x sin®) ( sin?6ko(m—cos®) 
= 8ra?(m— yf | 
(m—cos®)? 
sin?dko(m-+ cos®) 
(m+ cos®)? 


sin® 





| d@, (49) 


where x= 27a/Ao= ako. When ~x is very large and 6 is 
small, the main contribution to o, comes from values 
of ® in the neighborhood of 6=0. 

To approximate o, we expand the terms in the bracket 
of (49) in a power series in (cos#— 1) and neglect terms 
of order (cosh—1)? and higher. We let 6= dk» and find 


sinA(m—z) _ 
A(m—z) ~ A(m— 1) 





| sina(m—1) 


sinA(m—1 


+A(z— o| —coss(n—1) Of 6 pat], 


(m—1 
sinA(m-+-z) - 
A(m+z) A(m+1) 


sinA(m-+ 1) 
—A(z— | cosstmt » ow 1)?A*) . 
A(m+1) 





{sina (m+ 1) 


These expansions can be expressed in a slightly more 
compact form by introducing bessel functions of half 
integral argument: 








a : ) t7satm—1 

A(m—z) 2A(m—1) 
+A(2—1)J;(A[m—1])+---}, (50a) 

sinA(m-+z) ; 

A(m+z) toad Filet 1) 


—A(z—1)J,(A[m+1])+---}. (S0b) 
If we let z=cos and substitute (50a) and (50b) into 
(49), we find integrals 


¢,=82a?(m— ys | seat —1))s: 


2A(m—1) 
+2J,(ALm—1)Jy(ALm—1]) (52-51) + | 


+ [ere ta. 207 (alm 1D 


XJi(A[m+ 1])(se—si)+ -- +, (51) 
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where 
i J,2(x sind) 
7 i} hh tins (52a) 
sin® 
ix J°(x sin®) 
am f —_———-cos#d®, (52b) 
0 sin® 


Terms of order A‘ and higher have been omitted in 
(51). ; 

To integrate s2 we let u=sin® and then apply Eq. 
(3) on p. 136 in Watson to obtain 


52= f [JP (xu)/uldu=—3[JP(x)+Je(x)]J. (53a) 


We calculate s; by substituting the integral repre- 
sentation of J,?(x sin®), 


J ?(x sin®) , 


J ,(2x sin® sinu)sinu cos*udu, 
0 


= (4x/m)sin® 


which has been derived by Chaundy™ into (52b) and 
interchanging the order of integration. Then 


4x 
J (2x sin® sinu)d®. 


0 


4x 
Si:= (4x/7) f sinu cos*udu 
0 


But according to Watson, p. 374, Eq. (3) and p. 333, 
Eq. (3), 
hr 
J (2x sinu sin’)d&= (2x sinu)—'[[1—cos(2x sin) ]. 
0 


Hence (see Watson, p. 48, Eq. (1)), 


ir 
$1= (2/7) [1—cos(2x sin) |cos*udu 
0 


2 (/=){ («/4)— [1-01 costa 


=$[1—a7J,(2x)]. (53b) 


If we substitute (53a) and (53b) into (51), we obtain 
our final estimate to the total scattering cross section 
of a thin disk. As m—1, the first bracketed term (which 
is a function of A(m—1) gives the dominant contribu- 
tion to o,. For very thin disks (A=2r5/Ay—0 and 
x= 2ma/Xo) 


o,/ma?~4(m—1)?(m—1)*A*[1—(1/x)J1(2x)]. (54) 
(v) Total Cross-Section Curves 


In Fig. 1 we have plotted the total scattering cross 
section per unit area of geometrical cross section of an 


2 T. W. Chaundy, Quart. J. Math. (Oxford Series) 2, 144 (1931). 
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Fic. 1. Total scattering cross section per unit area of geo- 
metrical cross section, ¢,/2a of an infinite circular cylinder of 
radius a as a function of y=2(m—1)xa/Xo. Here m=Xo/A1, Ao 
= wavelength of incident plane wave, and \,= wavelength of in- 
cident plane wave in scattering material. 


infinite circular cylinder as a function of y= 2(m—1)a/ 
Ao as m—1. Notice that as y increases (that is, as the 
cylinder radius increases), o,/2a becomes a damped 
oscillating function of y which approaches 2 as (see 
Eq. (43)) yoo. 

In Fig. 2 we have plotted the ratio of the total scat- 
tering cross section of a finite cylinder of length L to 
that of a section of an infinite cylinder of the same 
length as a function of 27L/Xo (Eq. (48)). Notice that 
the end effect in a finite cylinder is less than 10 percent 
when 27L/Xo> 10. 

A quantity proportional to the total scattering cross 
section of a thin circular disk of radius a is plotted in 
Fig. 3 as a function of 27a/Xpo. 


APPENDIX I. SCATTERING BY SOFT INFINITE 
CYLINDER, ACOUSTICAL CASE 


In Sec. 2 the following expression was obtained for 
the scattered wave function for an infinite cylinder 
(see Eqs. (12b) and those below (15)). 


Ve= DY Ao, nt"H,™ (rho sinao)e*™*+éz*o cosa0, = (AT—1) 


n=—o 
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Fic. 2. Ratio of total scattering cross section of a finite cylinder 
of length L to that of a section of an infinite cylinder of the 
same length, o,(L)/o,(), as a function of Lko=2xL/Xo. 


where 
—m*J | (aky*)J .(ako*)+ gJ n(aki*)J n’ (ako*) 


~ m*H, (akc*) Jn! (aki*)— gJ n(aks*)[H a (akg*)] 





and the various symbols are defined in Sec. 2. Using the 
approximation (16) for the denominator of A,,n, we 
obtain 


wiag(ko*k:*)he'a(h*—¥0") [1 — ix(—1)etiehi*] 
(m*+ g)[1+« exp4iak,*] 
X {In (aki*)(m*/g)J n(ako*) 
—J (aki*)J n'(ako*)}. (AI-2) 





In the text we have summed (AI-1) only in the case 
of g=1 (that is, when the density of the acoustical 
scatterer is the same as that of the surrounding medium), 
By employing the method developed for scattering by 
spheres in S-I, we can sum (AI-1) in the general case of 
g~1. At distances sufficiently far from the scatterer s 
that the asymptotic expression for the hankel function 
is valid, we have 


A 


27 
v| : , exp[i(rko*— 4a+-zky cosa) | 
TT Ro 





XD Az, ne. (AI-3) 


n=>—@® 


The series in y, can be summed through the identities 


io) 


Jo(u)= > e'”?J,(y)J2(2), 


Jo(v)= > (—1)"e'"*J,.(y)J2(2), 
where 
u?= y+ 2°—2yz cos8, 
v= y’+2°+ 2yz cos. 
These series imply 
—Ji(u)u“"(y—z cosb)= > e'*J,'(y)Jn(z), (AT-4a) 


4) 


~J(o)v-(y-+-2 cosd—)= ¥ ei™Ty'(y)In(2). (ALD 


By combining (AI-4) with (Al-3) and (AI-2) and 
letting 
wr= kit kyo? — 2ki*ko* cosé, 
Vor = kt ko?+ 2ki*ko* cos6, 
we finally obtain the expression 
2rk,* (a?ko*)? 
R_ (m*+ g) 
+m*(1—g)cosO ]—ixuge?***1"[(g— (m*)?) 





\y.|?= wo J 1(awo)[(g—(m*)*) 





2 
— m*(1— g)cos@ ] / (1+ 2x cos4ak;*+F 





| 





Th 
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for the square of the modulus of the scattered wave 
function. The differential scattering cross section is 


o,(0)=R| ¥.(0) | . 
APPENDIX II. SOME MULTIPLE INTEGRALS 


Several multiple integrals are needed in various 
parts of this paper. The first is 


a Qa 
J(a’, b’) -{ f r exp{ir(a’ cosé 
0 0 


—b’ cosl@—¢]})}drd@. (A II-1) 


Let y and 7 be two new parameters defined by 
a’—b’ cosg=y siny, 
b’ sing=y cosy. 
Then 
y= a"+b”"—2a’d’ cose, 
a’ cos6—b’ cos(6@— v) = —y sin(6—y), 
so that 


J(a’, v= f rer sin(@—y) rdédr 
0 0 


= dn f rJulrydr= (2ra/y)Ji(ay). (A II-3) 


(A II-2) 


The second is somewhat more complicated: 


rT 2r 
ra)= f f einlar sinft cosé@—az2 sin{ cos(6—@ )+a3 cos{] 
aii -sintd¢d@. (A II-4) 


Let us define the three new parameters y, z, and 6 by 
a3=2 cosy, 


a\— a2 cosO =z sinB-siny, 
a2 sinO =z cos@ siny. 
Then 
2? sin*~= a;?+ a2?— 2a a2 cosO, 
2? cos*y=a;’, 
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Fic. 3. Total scattering cross section, o,, of a thin disk of 
radius a and thickness 6 as a function of x=2xa/Xo. Here 2 =a,/ 
4ra*A?(m—1)? and A=2765/Xo. 


so that 


2=([a+a2?+a3’— 20142 cos }}. 


We see that 


(A II-S) 


a; sing cos@— a2 sint cos(@— @) 
= —sin{[cos0(a2 cosO—a;)+ a2 sin® sind ] 
=z sint siny[_cos@ sin8— cos sin@] 
= —zsin¢ siny sin(@—8). 


Hence 


tg Qn 
Fa)= f f ein sing siny sin (6—8) giz cosy cost sintd¢d@ 
0 0 


rg 2x—8 
a f ei27 cosy cost singas f e~inz sing siny sind’ 79’ 
0 —B 
=2n f ei27 cosy cost J) (nz siny sint)sin¢dt 
. =2n(2e/nz)'Jy(nz). (A IL-6) 
The integral over 6’ is given in Watson, p. 20, Eq. (5). 


That over ¢ is given on p. 379, Eq. (1) with »=4 and 
r=0. 
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AC Characteristics of Si p-n Junction 
Sakujt KomaGata, Mitio Hatoyama, Mortoicu! SHipuyaA, WATARU 
SASAKI, TAKAAKI YAMAMOTO, AND MAKOTO KIKUCHI 
Electrotechnical Laboratory, Tokyo, Japan 
(Received June 4, 1951) 


OME ingots prepared from silicon, containing small amounts 

of unknown impurities, were revealed to be partly of p- and 
partly of n-type. Several samples of p-m junctions were cut from 
these ingots. Rectification ratio, taken for +1 volt, was from 80 
to 100, while a, the slope of the curve of logi versus » was as small 
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Fic. 1. Conductance, G, and susceptance, S, of a p-n junction 
as functions of frequency, f. 


as 15, in contrast to the value, 40, observed for Ge p-n junctions 
by Shockley e al.! The ac impedance of the samples was meas- 
ured for frequency range of 200 c/sec to 150 kc/sec. 

Figure 1 shows the dependence of the real part, G, and the 
imaginary part, S, of the admittance of a sample at 0°C, 0-bias 
voltage. For comparison, theoretical values by Shockley? are repro- 
duced with dotted lines in the same figure, but in arbitrary units. 
The result that the susceptance, S, approaches a constant value 
for higher frequencies seems to be in agreement with the result 
given by Odell and Fan* that the capacity of a p-m junction in 
Ge*exhibits a frequency dependence of approximately 1/f for a 
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Fic. 2. Resistance, R, 
and capacity, C, of a 
p-m junction as func- 
tions of bias voltage. 
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frequency range from 50 kc to 4 Mc/sec. It seems to be interesting 
too, that deviations from constant values of the resistance, R, and 
the capacity, C, of a p-m junction occur already at frequencies as 
low as 200 c/sec. 

Figure 2 shows the dependence of R and C on dc bias voltage 
for various frequencies at 0°C. The static differential resistance 
Ro, is also shown. The maximum slope of the curves of logR versus > 
decreases for increasing frequency, from 15 for Ro to 5.5 for 
150 kc/sec. 

1 Shockley, Pearson, and Sparks, Phys. Rev. 76, 180 (1949). 


2 W. Shockley, Bell System Tech. J. 28, 435 (1949). 
3N. H. Odell and H. Y. Fan, Phys. Rev. 82, 763 (1951). 





The Yielding Behavior of Iron Single Crystals 
A. N. HOLDEN 


Knolls Atomic Power Laboratory, Schenectady, New York 
(Received June 14, 1951) 


HERE have been three experiments recently’ in which 
small amounts of carbon or nitrogen were deliberately added 
to otherwise pure single crystals of iron, and the effect of these 
solute elements on the yielding of iron crystals was studied. In one 
experiment! no initial yield point was observed in the stress-strain 
curves of single crystals. In the other experiments*® initial yield 
points of small magnitude were observed. In all experiments strain- 
aging enhanced the yield point effect in the crystals tested. 

This seeming discrepancy in the behavior of iron single crystals, 
of course, has made the formulation of an adequate theory of the 
yield point more difficult. Therefore, the study of the behavior of 
these carburized crystals has been continued with the following 
results: 

1. No yield point has been found with completely annealed iron 
crystals containing carbon under conditions of very accurate align- 
ment during testing. 

2. Very small amounts of plastic strain introduced intentionally 
into a carburized crystal followed by 100°C annealing for as short 
a period as one hour caused the crystal to exhibit a yield point on 
testing. For example, 0.25 percent strain followed by the 100°C 
anneal caused a crystal to have sharp yield point and 0.4 percent 
yield point elongation. After 7.5 percent strain and one hour at 
100°C the yield point elongation was increased 0.75 percent. 

3. There is a size effect in the yield point elongation for strain- 
aged crystals. If two identically oriented strain-aged crystals are 
chosen that differ appreciably in cross section, the larger crystal 
will have the more pronounced yield point and the greater yield 
point elongation. 

4. If a carburized single crystal is strained by light rolling, aged, 
and tested in tension, it shows no yield point. If an identical crystal 
is strained in tension, aged, and retested in tension, it exhibits a 
marked yield point. 

A satisfactory explanation of these experimental results requires 
that the Cottrell theory* of solute-anchored dislocations be ex- 
tended to describe in detail the release of an avalanche of disloca- 
tions at the yield point. Cottrell’s most recent extension of his 
original hypothesis’ requires that the dislocations which break 
free of their solute atmospheres at the upper yield stress become 
piled up behind polygonization boundaries until the stresses on 
these boundaries are high enough to release an avalanche. There 
is no evidence of polygonization boundaries in undeformed an- 
nealed single crystals. Furthermore, it seems unlikely that a crys 
tal, deformed and aged at room temperature, would polygonize, 
yet such low temperature strain-aging will cause a yield point. 

An alternative to Cottrell’s boundary interaction mechanism is 
the interaction of dislocations by collision. The yield point in 
dislocation collision process may be described as follows: At equi- 
librium, carbon atoms have anchored all the randomly distributed 
dislocations (boundaries are not required). When the crystal is 
then stressed in tension, several dislocations may be torn free of 
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their carbon atmospheres at the upper yield point. Each of these 
free dislocations accelerates to a velocity approaching the velocity 
of sound. After moving freely some distance they pass close to 
other anchored dislocations. The interaction forces between the 
moving and anchored dislocations as they pass are sufficient to 
knock the latter free of their solute atmospheres. A repetition of 
this process frees an avalanche of dislocations. 

Simple calculations have shown that the interaction forces are 
sufficient to knock anchored dislocations free if the moving disloca- 
tions pass within three or four atom distances of the anchored 
dislocations; suck events are considered successful collisions. 
Furthermore, the moving dislocation has enough kinetic energy 
to impart the required interaction force to the anchored dislocation 
if it moves at an appreciable fraction of the velocity of sound, a 
not improbable velocity.* A complete quantitative treatment of 
the dislocation collision theory has been submitted for publication 
elsewhere. 

As a consequence of the required distance (three or four atom 
distances) for successful collision, it can be shown that an annealed 
crystal containing 10* dislocations? per cm? would have a mean 
free path between collisions of 10! cm, whereas a strain-aged 
crystal containing 10" or 10” dislocations’ per cm? would have a 
mean free path between collisions of 10~* or 10-* cm. 

One may conclude that the yield point is small in an annealed 
crystal because the mean free path between dislocation collisions 
approaches the specimen size. One may conversely explain a more 
extensive avalanche of dislocations and larger yield point effect 
in a strain-aged crystal. Furthermore, the existence of a crystal 
size effect on the magnitude of the yield point can be readily 
explained by the mean free path concept, since the longer the slip 
path through the crystal the greater the number of successful 
collision events. 

The experiments in which rolled crystals did not have yield 
points when aged and tested in tension (while crystals deformed 
small amounts in tension) aged, and retested in tension did have 
yield points, indicate that dislocations added along slip planes 
operative during rolling do not affect the mean free path for dis- 
location collision on slip planes operating during tension. 


1A. N. Holden and J. H. Hollomon, Trans. Am. Inst. Mining Met. Engrs. 
185, 179 (1949). 


2A. H. Cottrell and A. T. Churchman, Trans. Am. Inst. Mining Met. 
Engrs. 185, 877 (1949). 

3H. Schwartzbart and J. R. Low, Trans. Am. Inst. Mining Met. Engrs. 
185, 637 (1949). 

4A.H. Cottrell, Report of Conference on Strength of Solids, Bristol, Proc. 
Phys. Soc. (London) (1948) 30. 

5A. H. Cottrell, Sympostum on Plastic Deformation of Crystalline Solids, 
Mellon Institute, Pittsburgh (1950) 60. 

*F. C. Frank, Symposium on Plastic Deformation of Crystalline Solids, 
Mellon Institute, Pittsburgh (1950) 89. 

7 There are many estimates of each value. See A. H. Cottrell, Progress in 
Metal Physics, Vol. I, 77 (1949). 





Catalysis by Metal Powders 


DonaLp P. SMITH 
Montrose, Pennsylvania 
(Received June 18, 1951) 


ROM many x-ray studies, and from behavior with hydro- 

gen,'? it appears that massive metals in their usual some- 
what strained (cold-worked) condition, owe important properties 
to widened lattice intervals and the stressed fields of interatomic 
forces within these intervals. It may not be amiss to call attention 
to reasons for regarding finely divided metals, obtained in the 
ordinary way, as merely strained metals in small masses, which 
also owe certain characteristics to their departure from the state 
of the ideal lattice, their disstructure or rifted condition. 

Metallic powders, blacks, or sponges are in general prepared by 
chemically withdrawing atoms from compounds (or alloys) of the 
metal desired, at temperatures below thosé at which the formation 
of the new crystal readily proceeds, its recrystallization range. 
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Metal thus set free, and not then or later subjected to higher 
temperature, forms its own proper lattice incompletely and has in 
part a skeletal arrangement derived from the oxide (etc.) from 
which it came. The presence of a certain amount of perfected 
metallic lattice is shown by corresponding lines in the x-ray diffrac- 
tion spectrum; while the marked broadening of these lines may 
result either from fineness of particles, or from local widenings of 
spacing within the continuous portions of lattice, “parameter 
variation.” As in the much-studied case of massive metal after 
cold-working, the x-rays fail to distinguish clearly between these 
two causes.’ 

Heated to somewhat higher temperatures, either during or after 
their formation, such powders show much sharpened diffractions 
and are less active in all processes involving hydrogen—changes 
often attributed solely to coalescence of minute particles to larger 
ones, with accompanying loss of surface area, but which must be 
due, in part at least, to the closing of widened intervals resulting 
from formation of sound lattice. The fact is not to be overlooked 
that the narrowing of lines and loss of activity take place in the 
same range of temperatures in which recrystallization and the 
disappearance of the cold-worked and rifted state occur in strained 
massive metal. 

During atom-withdrawal, a gradual building-up of metal results 
which is the reverse of the breaking-down of lattice in the course 
of plastic deformation. The metal goes through the open-lattice 
stages of mechanical deformation, but in the opposite order; and 
although the metal is itself obtained in highly subdivided state, 
since the starting material is a finely-divided substance, yet condi- 
tions within the individual particles are similar to those which 
exist in cold-worked massive metal. Each particle is filled with 
widened intervals. 

In the strained massive metal, mechanical work has forced some 
of the atoms from their positions of equilibrium in the normal 
lattice and has formed rifts, or locally widened spacings; while in 
the powder which results from atom-withdrawal a chemical process 
has left gaps where the new lattice has not yet taken form. The 
gaps may well differ from the rifts in arrangement and shape, but 
the two types of opened spacing are essentially alike as to basic 
properties. Both are intervals only slightly wider than the normal 
lattice interstices; and both are places of intensified and inhomo- 
geneous force fields, since in each instance work (mechanical or 
chemical) has been done in what amounts to a removal of atoms 
to a distance from their neighbors greater than that at which they 
would rest in perfect lattice. 

If now we imagine a hydrogen atom introduced into such a field, 
it follows that the forces exerted upon electron and proton in 
opposed directions are great in certain places, but fall to very small 
values at some other points within the rift. At the one point, then, 
a hydrogen atom will be ionized, while at the second electron and 
proton combine, so that hydrogen exists in different (simultaneous) 
degrees of ionization. The widened lattice interval is hence an 
environment in which processes of electron transfer are induced 
and hydrogen is activated. 

Finely-divided metal displays activity as catalyser, particularly 
in regard to reactions involving hydrogen. Since in this state the 
metal has a very large surface, it is commonly assumed that reac- 
tion takes place on the surface; but it is to be noted that if the 
solid particles contain widened intervals, the accessibility of these 
rifts, to an outer phase, also increases with surface area. Mere 
proportionality of reaction to surface, even if definitely proved, is 
therefore ambiguous evidence as to the seat of activity. 

From examination of copious empirical information,? it is evi- 
dent that in iron and the other metals of its (endothermically 
absorbing) type, in massive form, absorbed hydrogen is contained 
only in the widened intervals and does not penetrate the normal 
(un-work-altered) lattice; while the investigations of Sieverts* 
make it equally clear that the solubility of hydrogen in iron is the 
same whether the metal is in the form of wire, foil or powder. With 
regard to the higher ‘temperature range (recrystallization range), - 
where alone a solubility equilibrium can be reached,? the conclusion 








1292 LETTERS TO 


is hence unavoidable that, in the finely-divided metal also, all 
hydrogen taken up from the gas phase is contained in the widened 
intervals of the crystal particles. Nothing remains for surface 
condensation, since otherwise the solubility observed for powder 
would be greater than that found for foil. 

In the temperature range below recrystallization, evidence of 
just this kind is not applicable; yet the facts are highly instructive. 
Thus, in the case of the well-studied iron, much-strained massive 
metal in the form of foil has been found to take up 5.7 ml/gram, 
at 25°C and one atmosphere.’ The illuminating thing is that this 
large supersaturation of massive metal in rifted condition is many 
times greater than the amount adsorbed by finely divided iron.® 
Although, as we have seen, the capacities of the two forms of iron 
are the same at higher temperatures, passage to the lower range 
results in a greater increase of capacity, not for the powder but for 
the cold-worked massive metal! The absorptive capacity of the 
widened lattice intervals is evidently much larger in severely 
strained massive metal than in powder; yet the capacity of the 
deeper-lying rifts is realized only slowly, as hydrogen diffuses 
inward by comparatively long paths; whereas the smaller rift- 
capacity of the powder, being all of it relatively superficial—near, 
though not on, the surface—is more quickly attained. Hence the 
activating (ionizing) action which has effect only within the rifts, 
is far more quickly available to an outer phase, in the case of the 
finely divided metal; and the extended surface of the powder 
obtained by atom-withdrawal is essential to the high velocity of 
interaction between phases which is necessary for contact catalysis. 

If the contact activity of finely divided metal is dependent on 
the stressed and unstable state of the rift-fields, it is clear that 
the energy stored during their formation must contribute a part 
to any contact reaction, with in consequence a gradual relief of 
the unstable condition and progressive decline of activity. The 
well known decay of such catalysts may accordingly be due to the 
fact that their action is not wholly “catalytic,” rather than to 
accumulation of impurities from the reaction mixtures! 

The prevailing conception of contact catalysis by finely divided 
metals is based almost exclusively upon study of such reactions. 
By regarding the metal powder as strained metal in small masses, 
and attributing its activity to the same characteristics of widened 
lattice intervals which are evident in larger masses, instead of to 
special properties of surfaces and adsorbed layers, we should gain 
the advantages of resting our considerations upon a much wider 
and more varied basis of facts, and of bringing outwardly diverse 
processes under a common point of view. 

The considerations so briefly outlined above seem therefore to 
merit more careful examination. 

1D. P. Smith, Hydrogen In Metals (Chicago University Press, 1948). 

2D. P. Smith, ~~ 5 113, 348-52 (1951). 

* Stokes, Pascoe, and Lipson, Nature ist. 137 (1943). 

4 Sieverts and Hagen, Z. physik. Chem. A155, 314- vy (1931). 

’G. A. Moore, Trans. Am. Inst. Mining Met. Engrs. 135, 255-91 (1939). 

*W. A. Dew and H. S. Taylor, J. Phys. Chem. 31, 277-90 (1927). A. F 


Benton and T. A. White, J. Am. Chem. Soc. 52, 2325-36; 53, 3301-14; 54, 
1820-30 (1930-31). A. F. Benton, Trans. Faraday Soc. 28, 202-18 (1934). 





Erratum: On the Spectrum of Isotropic 
Temperature Fluctuations in an 
Isotropic Turbulence 
[J. Appl. Phys. 22, 469-473 (1950)] 


STANLEY CORRSIN 
Department of Aeronautics, The Johns Hopkins University, 
Baltimore, Maryland 


T has been called to my attention by Professor M. Z. v. 
Krzywoblocki that for low speed flow of an ideal liquid 


(p=constant) ¢, is replaced by c, (specific heat at constant vol-. 


ume) in the approximation to the energy equation. The analysis 
and results as given apply to very low Mach number flow (M—0) 
of a perfect gas.' 


1S. Corrsin and L. S. G. Kovasznay, J. Aeronaut. Sci. (to be published). 
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Tropospheric Propagation beyond the Horizon 
J. FEINSTEIN 
National Bureau of Standards, Washington, D. C. 
(Received June 23, 1951) 


EASUREMENTS on vhf transmissions in recent years have 
indicated that field strengths deep in the diffraction region 
are many orders of magnitude greater than values calculated on a 
smooth earth diffraction theory, even when allowance is made for 
the refraction produced by the standard atmosphere.' Of the 
theories advanced to date in explanation of this effect, the 
mechanism of scattering by inhomogeneities resulting from atmos. 
pheric turbulence has perhaps been the most successful.? However, 
certain discrepancies in wavelength and angle dependence, to. 
gether with the fact that the abnormal fields apparently exist 
everywhere at all times, lead one to seek an explanation solely in 
terms of the model represented by the standard atmosphere. 

It is well known that a discontinuity in the characteristics of a 
medium traversed by a wave gives rise to a reflected as well as a 
transmitted wave. Such reflections have been invoked in explana- 
tion of certain tropospheric effects.’ It is usually overlooked that a 
gradual change in medium properties produces the same effect in 
principle, although the magnitude of the continuously reflected 
components is then generally quite small. Since the atmosphere 
possesses a gradient of refractive index, such partial reflection 
must occur. By this mechanism wave energy is delivered to areas 
far beyond the line of sight of the transmitter. Since the diffracted 
field decays exponentially with distance in this region, while the 
reflections go down at first only as an inverse power of the distance, 
the latter dominate rather rapidly. 

It is perhaps surprising in view of the many wave treatments of 
this problem extant that these contributions have not been taken 
into account. The answer appears to lie in the approximations 
employed in evaluating the results of rigorous formulations of the 
problem. The WKB method, or one equivalent to it, is generally 
employed to approximate the wave functions characteristic of the 
inhomogeneous atmosphere when evaluation of the contour inte- 
gral or other type of representation of the field is undertaken. It is 
well known that WKB solutions represent traveling waves for 
which the loss in intensity occasioned by the continuous partial 
reflections is taken into account, while the new waves arising from 
these reflections are entirely neglected.‘ 

Rigorously, a true wave solution is required for an exact answer, 
but it appears feasible to estimate the magnitude of the field 
produced by these reflections from a hybrid ray and wave theory, 
along the following lines: The contributions to the field at any 
receiving point will come from those partial reflections originating 
in areas for which the equal angle law of specular reflection is 
satisfied. When the transmitter 7, and the receiver R, are on the 
surface of the earth this locus is the perpendicular bisector of the 
path TR. Across any differential distance dz along this locus we 
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Fic. 1. Dependence of signal produced by atmospheric gradient 
reflection upon distance and antenna height. 
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consider the beam incident upon the troposphere to be a plane 
wave whose angle of incidence a with respect to the spherical 
stratification is a function of the geometry and of atmospheric 
refraction; if, further, we replace the spherical surfaces of constant 
refractive index by plane stratification, and take account of 
the convergent effect upon the wave front through a coefficient 
C(a), then it may be shown that the differential reflection is 
(Kdz/2a*)C(a) for small (but not too small) a,5 where K =dn/dz, 
the gradient of refractive index. Adding the contributions all along 
: with their proper phase: 


2(a9) ae exp(ikza) - C(a) - fr(a)- fr(a)dz, (1) 





where the f functions are the effective antennae patterns. 

The validity of this treatment is contingent upon its develop- 
ment from a rigorous wave formulation. 

For a “high” antenna, the pattern function may be replaced by 
its average value over a lobe, while for a low antenna, one may 
employ k(a—ao) as pattern function. Consequently curves (a) 
(Fig. 1) are calculated from 


| R| ~LK (eo) ]/ (Akexs*)- (2/m)?-C, 


where C is the mean value of the convergence coefficient; for 
curves (b), 


|R| YK /(4ke*)+(2/m)-(n/D)-C. 


1 Proc. Inst. Radio Engrs. 39, 391 (1951). 

2 Booker and Gordon, Proc. Inst. Radio Engrs. 38, 401 (90%. 

3K. A. Norton, Advances in Electronics, I, (1948) p. 402 ff 

4H. Bremmer, Terrestrial Radio Waves (1949) p. 158. 

’ By a modification of the work by H. Bremmer, Physica 15, 593 (1949). 





The Electrical Resistivity of Cold-Worked 
Copper-Iron Alloys 


F. E, HETHERINGTON AND JAMES REEKIE 


Department of Physics, Royal Military College of Canada, 
Kingston, Ontario 


(Received July 6, 1951) 


HILE an annealed metal of high purity has, as a rule, a 

well-defined value of electrical resistivity, this resistivity 
may be profoundly affected by the addition of even minute 
amounts of impurities or by internal strains present in the pure 
metal. In general these factors result in an increased resistivity, 
and it seems probable that impurities and internal strains both 
give rise to randomly distributed scattering centers which are the 
basic cause of the increase in resistivity. 

In this paper we wish to report briefly on measurements of 
electrical resistivity made on a series of alloys of copper with small 
amounts of iron, the alloys being used both in the annealed state 
and after various degrees of cold-working carried out by drawing 
through dies. From the results some significant conclusions can be 
drawn regarding the behavior of the iron impurity when subjected 
to stress in the metal. The experimental method was somewhat 
similar to that used by Rutter and Reekie! in previous work on 
pure metals, but the accuracy of measurement has been consider- 
ably improved by modifications to the galvanometer and potenti- 
ometer circuits. 

The lowest curve of Fig. 1 shows how the electrical resistivity p 
of the annealed alloys varies with the amount of iron present. p 
increases rapidly from a value of 1.75X10-* ohm-cm for pure 
copper (Johnson, Matthey, Ltd., ““H-S” brand, iron content less 
than 0.001 percent) to 3.49X 10-8 ohm-cm when 0.26 percent of 
iron is present in the copper. For alloys with greater amounts of 
iron there is a sharp decrease of resistivity followed by a further 
slow rise. The general shape of this curve is maintained even when 
the alloys are cold-worked, as is shown by the results obtained for 
36, 61, 82, and 95 percent reduction of area. There is, however, a 
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noticeable accentuation of the peak at an iron content of 0.26 
percent. All these curves were measured at a temperature of 
27.0°C. 

During annealing the alloys were slowly cooled in the furnace 
from 850°C down to a room temperature of 25°C, and we can 
interpret the curve for the annealed metal as indicating a solu- 
bility for iron in copper of approximately 0.26 percent at 25°C. 
This result is in reasonable agreement with the early work of 
Hanson and Ford? With higher iron content the iron is precipi- 
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Fic. 1. Electrical resistivity of copper-iron alloys for various degrees of 
cold-working. (Cold-working shown as percent reduction of area.) 


tated on slow cooling, and the peak in the resistivity curve prob- 
ably comes about in part as a result of stress relief in the copper 
lattice when the iron begins to precipitate. 

By taking the value of resistivity observed for each alloy at a 
given degree of cold-working, and subtracting from this the value 
observed for the resis tivity of pure copper at the same degree of 
cold-working, we can obtain the quantity [pcu-re—pcu]. If we 
now plot this quantity against the degree of cold work, we will 
obtain a measure of the effect of cold-working on the iron in the 
state in which it is present in the alloy. Carrying out this process 
for several of the alloys used we obtain the results shown in Fig. 2. 
It will be at once evident that there is a significant change in the 
character of the curves as we proceed from low to higher iron 
conténts. For the alloys containing less than 0.15 percent iron the 
resistivity is increased by a constant amount, irrespective of the 
degree of cold work. The alloys with 0.15 and 0.23 percent iron 
show constant values of [pcu-re—pcu] until the cold-working has 
reached 80 percent and 30 percent reduction of area, respectively. 
Alloys with iron contents of 0.26 percent or more show a continu- 
ous increase in [pcu-Fe— cu]. From these curves we conclude that, 
if the iron is in solution, cold-working does not change that portion 
of the resistivity contributed by the iron, the effect of the iron 
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Fic. 2. Effect of cold-working on the component of resistivity because of 
the iron in copper-iron alloys. (Weight percent iron as shown.) 


merely being to increase the number of randomly distributed 
scattering centers, or, possibly, the number of dislocations. If, 
however, the iron is precipitated, then cold-working does increase 
the resistivity due to the iron. Furthermore, we can see that in the 
case of the alloys with 0.15 and 0.23 percent iron sufficient cold- 
working results in a partial precipitation of the iron. This would 
be expected to occur when the strain energy due to the deformation 
exceeded the activation energy necessary to allow migration of 
the iron atoms. 


1J. W. Rutter and J. Reekie, Phys. Rev. 78, 70 (1950). 
2 D. Hanson and Grace W. Ford, J. Inst. Metals 32, 335 (1924). 





A Note on a Simple Technique for Obtaining 
Absorption Coefficients of Metals and Alloys 
for Complex Radioactive Spectra* 


FRANK E. Jaumot, Jr.,f AND Foster C. Nix 


Randal Morgan Laboratory of Physics, University of Pennsylvania, 
Philadelphia, Pennsylvania 


(Received July 25, 1951) 


CCASIONALLY it is desirable to know the absorption 

coefficient of a metal or alloy whose physical properties 
make it inconvenient or impossible to measure this quantity by 
the usual method of using thin foils. This is especially true when 
the absorption coefficient is desired for a complex spectrum, so 
that calculation of the quantity is difficult and uncertain. As an 
example, in the course of doing some work on diffusion of cobalt 
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Fic. 1. A plot of activity as a function of the thickness 
of absorbing layers of cobalt-aluminum alloy. 


in cobalt-aluminum alloys, the writers found it necessary to 
know the value of the absorption coefficient of very thin cobalt. 
aluminum alloys for the complex spectrum! of Co®. These alloys 
are so brittle that the use of foils was prohibited. 

It was found possible to measure the absorption coefficient by 
dispersing powdered alloy in a solution of polyvinyl alcohol. The 
polyvinyl alcohol was obtained as Elvanol, the Du Pont trade 
name for the product, and a saturated water solution prepared 
according to the manufacturer’s directions. This solution was 
diluted with an equal volume of distilled water before the powdered 
alloy was added in order to insure a sufficiently thin mixture. The 
mixture of polyvinyl alcohol and metal was poured onto a flat 
metal plate and when dry could be peeled off with ease, thus pro- 
ducing a foil composed of polyvinyl alcohol and powdered metal 
for which the absorption could be measured. The mass per sq cm 
of the foil was determined. At the same time a foil of pure poly- 
vinyl alcohol was made and the mass absorption coefficient deter- 
mined. Polyvinyl alcohol decomposes into its monomer, a gas, at 
temperatures above 150°C, thus allowing us to bake the alloy 
pigmented foil in a furnace until only the powdered alloy residue 
remained. From this we determined the mass per sq cm of the 
powdered alloy. 

Making the assumption that the absorption mechanism in the 
pigmented foil was identical to absorption through two layers, one 
of polyvinyl alcohol, the other of powdered metal, we corrected 
for the absorption of the polyvinyl alcohol and obtained the ab- 
sorption coefficient for the alloy by plotting N; against x, where 
N,, the calculated activity through the powdered alloy is given by 


Ni=Ne, 


and x; is the mass per sq cm of the powdered alloy. WN is the ac- 
tivity observed through the pigmented foil, uz is the absorption 
coefficient, and x2 the mass per sq cm of polyvinyl alcohol. N, of 
course, is given by 


N=Nee**, 


where Np is the activity of the source without absorbers, and 
and x are the absorption coefficient and mass per sq cm of the 
pigmented foil. 

The metal should be powdered as finely as possible, not larger 
than 300 mesh, and the mixture should be thoroughly stirred so 
that reasonably uniform foils are obtained. 

This method was checked by determining the absorption coefi- 
cient of powdered cobalt for the complex spectrum of Co™, which 
had been previously determined for cobalt foils,? and by comparing 
the results obtained for powdered aluminum using the aforemen- 
tioned technique, with the absorption coefficient obtained using 
thin aluminum foils. The values obtained for cobalt agreed to 
within three percent, and those for aluminum to within less than 
one percent, both being well within the experimental error involved 
in absorption measurements. 

As an illustration of the results obtainable by this method, we 
give in Fig. 1 the plot of the activity, Ni, in counts per minute, 
as a function of the thickness in g per cc, of a cobalt-aluminum 
alloy near the 50-50 atomic percent composition. 

* This work was supported by the ONR. 

+ This note is part of a dissertation by Frank E. Jaumot, Jr., presented 
to the faculty of the Graduate School of the University of Pennsylvania 1 
partial fulfillment of the requirements for the Ph.D. de 


gree. 
1 AEC, Isotopes Branch: Catalog and Price List No. 3 (July, 1949). 
2F. C. Nix and F. E. Jaumot, Jr., Phys. Rev. 82, 72 (1951). 
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Edge Effects in a Circular Aperture 


RoBERT E. HOUSTON AND ROBERT H. NOBLE 


Department of Physics, Michigan State College, 
East Lansing, Michigan 


(Received July 2, 1951) 


N a program of electromagnetic field measurement in small 
apertures, a set of data was taken with the same diameter- 
wavelength ratio (d/X) =2 as one of Andrews.' The actual wave- 
length (50 cm) was somewhat larger than Andrews’ wavelength 
(8 cm) so that we were able to investigate the field closer to the 
edge of the aperture. The rest of our experimental procedure was 
similar to his except that we used point by point readings of a 
galvanometer in place of his recording system. 
Our data, as shown in Fig. 1, were found to be in very good agree- 
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Fic. 1. Field measurements along the magnetic diameter and along the 
electric diameter of a circular aperture (d/A) =2. Ordinates are in arbitrary 
units proportional to galvanometer deflections. 


ment with his curves except at the edge of the aperture. Here the 
field rises sharply and approximates the function (a?—7?)~$ which 
would be predicted in this region from the work of Bethe.” 

A more complete description of these results for a series of aper- 
tures will be published in the near future. 

We are grateful to Dr. R. D. Spence for suggesting this problem 
and for many helpful discussions. 


1C. L. Andrews, J. Appl. Phys. 21, 761 (1950). 
2H. A. Bethe, Phys. Rev. 66, 163 (1944). 





Structure of Oxide Replica for Electron Microscopy 


SHIGETO YAMAGUCHI 
Scientific Research Institute, Bunkyo-ku, Tokyo, Japan 
(Received July 9, 1951) 


T is known that the oxide replica of metal gives often a more 
resolved image than Formvar replica in electron microscopy.! 
The oxide films formed on metals and able to reprint faithfully 
their surface forms are called “structureless.” It has been shown 
in the present study that the structureless oxide film formed on 
the etched surface of iron single crystal is composed of oriented 
crystals of small size (50-100A). 
The specimen of iron single crystal used here was a piece of wire 
(diameter: 2 mm, length: 10 mm). This was etched with an ethanol- 
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Fic. 1. Oxide film formed on etched surface of iron single crystal. 


bromine (10:1 by vol.) solution for 10-30 seconds. The oxide 
replica was formed on this etched surface by the method of the 
author? and stripped off the substrate by Mahla-Nielsen’s method.’ 
The micrograph of this replica film is given in Fig. 1. In Fig. 1 
there are the octahedral crystals oriented parallel to each other 
and inclined by 25°-30° regarding the (001) plane to the replica 
film. The electron diffraction pattern from the etched surface of 
the specimen showed the same inclination of (001) plane to the 
macroscopic surface as that in microscopic observation. This coin- 
cidence verified the faithfulness of oxide reprint performed here. 

The electron beam running perpendicular to the replica film of 
Fig. 1 gave a diffraction pattern of Fig. 2, which verified the 
existence of regularly oriented Fe;O, crystals. The oxide crystals 
formed on metal single crystal orient themselves according to the 
orientation of the substrate.* The orientation of the oxide particles 
formed on one boundary face (e.g., (111) plane) of single crystal 





Fic. 2. Diffraction pattern obtained with the beam perpendicular 
to the film of Fig. 1. 
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should be distinguishable from that of oxide particles of another 
face (e.g., (111) plane). Therefore, the diffraction spots obtained 
from the oxide replica might appear in the positions in a focal 
plane of objective lens, which are determined by the orientation 
of the oxide particles to the incident beam. Since the diffraction 
pattern obtained on a focal plane of microscope behaves as a new 
diffraction lattice for the observable image (Abbe’s principle),® the 
oxide replica of Fig. 1 giving separated diffraction spots of Fig. 2 
on a focal plane should give a more resolved image than Formvar 
replica. 

1 Metallurgical Application of the Electron Microscope, Institute of 
Metals, Monograph and Reports Series No. 8 (1950). 

2S. Yamaguchi, J. Appl. Phys. 22, 680 (1951). 

3 E, M. Mahila and N. A. Nielsen, J. Appl. Phys. 19, 378 (1948). 

4H. R. Nelson, J. Chem. Phys. 5, 252 (1937). 


Svan Dorsten, Nieuwdorp, and Verhoeff, Philips Tech. Rev. 12, 33 
(1950); Haine, Page, and Garfitt, J. Appl. Phys. 21, 173 (1950). 





Fracture and Yield Stress of 1020 Steel 
at Low Temperatures 


A. S. ELp1n* anp S. C. COLLIns 


Laboratory of Cryogenic Engineering, Massachusetts Institute 
of Technology, Cambridge, Massachusetts 


(Received July 10, 1951) 


HE large helium refrigerator recently completed at the 
Massachusetts Institute of Technology has been equipped 
with a tensile testing machine capable of applying loads up to 
60,000 pounds. A series of investigations of some of the mechanical 
properties of metals at low temperatures has been planned. The 
first project undertaken and herewith reported was the result of a 
suggestion by Professor E. Orowan of the need for measurements 
of brittle strength as a function of temperature. 

All of the specimens were cut from a single bar of hot-rolled 
1020 steel just as it came from the mill. At first, standard one-half 
inch samples were used. At reasonably low temperatures, however, 
breakage always occurred in the threaded section. Thereafter the 
diameter was reduced to 0.320 inch, a figure which was found by 
experience to be the largest section which could be depended upon 
to break in the test section rather than in the region of the three- 
quarter inch threads. 

Several samples were tested at various temperatures in the 
interval 12°K to 61.5°K. All exhibited typically brittle fracture 
with no reduction in area. At 75°K, however, a definite yield point 
was observed and the increase in length before rupture was un- 
expectedly large. Several additional specimens were broken at 
temperatures lying between 61.5° and 75°K. The results indicate 
a sharp transition near 61.5°K below which no reduction of area 
occurs and above which the reduction of area increases rapidly 
as the temperature at which the sample is broken rises above 61.5°. 
These results are shown in Figs. 1 and 2. 
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The fracture stress was computed from the breaking load and 
the cross-sectional area at the zone of fracture. This is plotted 
against temperature in Fig. 3. A maximum at somewhat above 
100°K is observed. The magnitude of the strain before breakage, 
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a negative temperature coefficient of the effect of strain hardening 
upon strength, and a decreasing brittle strength with temperature 
could account for the shape of the curve between 61.5° and 100°K. 
Throughout this temperature range the fracture appeared to be 
wholly brittle. At 104° the broken section contained at its center a 
dark, fibrous dot about two millimeters in diameter, the remainder 
of the area having a brilliant crystalline appearance. As the tem- 
perature at which the specimen was broken was increased, the size 
of the darkened circle increased until it covered the whole area at 
185°K and gave the appearance of a fracture that is character- 
istically ductile. 

In order to avoid bending moments in the specimen resulting 
from lack of symmetry in the assembly when force is initially 
applied, three strain gauges were mounted on each specimen and 
by repeated trials an adjustment was found which gave sym- 
metrical loading. A small load of a few hundred pounds was kept 
on the sample while the tensile machine was lowered into the cold 
well and during the cool-down period so that the adjustment 
would not be lost. This precaution undoubtedly helped to reduce 
scattering of the results. 
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The load is applied to the specimen by a hydraulic cylinder at 
the warm end of the assembly. A load cell is interposed between 
the hydraulic cylinder and the pull rod which is attached to the 
specimen. The load cell is equipped with strain gauges which 
have been calibrated. 

The temperature down to 20°K is measured by a copper- 
constantan thermocouple. The cold junction is soldered to a 
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Fic. 4. General view of tensile testing machine and helium cryostat. 


massive part of the testing machine near the center of the speci- 
men. It is possible that the temperature of the specimen was 
slightly different from that recorded. Below 20°K a hydrogen vapor 
pressure thermometer is used. At the 60° level the sensitivity of 
the thermocouple circuit used is 0.2°, but there could be an error 
of 0.5° in estimating the temperature. At the 30° level the error 
could be as high as 3°. 

Figure 4 shows the tensile apparatus ready. to be lowered into 
the 8-foot well of the refrigerator. 

Acknowledgment of the valued assistance of Professor R. B. 
Green is cordially made. 

Facilities and funds for the work reported were provided by the 
ONR. 


* Submitted in partial fulfillment of the requirements for the D.Sc. degree. 





A New Method for the Production of Air Shocks 
in a Tube* 
A. B. LAPONSKY AND R. J. EMRICH 


Lehigh University, Bethlehem, Pennsylvania 
(Received July 19, 1951) 


N a conventional shock tube,’ shocks are produced by bursting 

a diaphragm which separates a high from a low pressure sec- 
tion of the tube. This arrangement produces almost immediately 
an essentially uniform flow in the tube behind a shock of constant 
strength. 

An interesting aspect of the flow is the formation and growth of 
the shock. The bursting diaphragm is not well adapted to the 
study of these processes because the diaphragm bulges, breaks in 
an erratic manner, and permits the shock to form at nearly full 
strength almost immediately. A means for initiating a transient 
flow in which flow variables change over measurable distances and 
times is required. The acceleration of a relatively massive piston 
at one end of a gas-filled tube produces such a flow. 
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A piston accelerated from rest at 10* m/sec? for a distance of 1 
meter initiates in air a flow in which a shock will appear within 10 
meters from the initial piston position. Such a flow is being pro- 
duced by pneumatically impelling a 50-gram piston in a 3.5-cm- 
diameter tube. Piston and shock positions as functions of time 
are being measured. The predictions of one-dimensional flow 
theory? are, to a first approximation, confirmed by the results. 


* This work is being supported in part by the ONR under contract 
N7onr-39302. 


1 Bleakney, Weimer, and Fletcher, Rev. Sci. Instr. 20, 807 (1949). 
2R. Courant and K. O. Friedrichs, Supersonic Flow and Shock Waves 
(Interscience Publishers, Inc., New York, 1948), p. 168. 





A Method for Determining the Energy Stored 
During Cold-Working of Metals 


M. B. BEVER AND L. B. TICKNOR 


Department of Metallurgy, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 


(Received August 11, 1951) 


HE energy which metals retain as the result of cold-working 
is an important factor in the deformation process as well 
as in some processes, such as recrystallization, which may occur 
subsequently to cold-working. Methods used for measuring this 
energy have included comparisons of heating curves of a worked 
and an annealed sample and determinations of the heat of solution 
of such samples in aqueous media; the stored energy has also been 
found as the difference between the total work of plastic deforma- 
tion and the heat evolved during cold-working. While the spread 
of published values'—ranging from less than 0.1 to more than 10 
cal per gram of metal—may be partly due to differences in the 
manner and degree of cold work and in the kind of metal, there 
exists a real need for more convenient methods of satisfactory 
accuracy. 

A new calorimetric method employs the fact that the heats of 
solution are small in some metallic systems, especially compared 
to heats of solution of metals in aqueous media. The heat of solu- 
tion is measured in a calorimeter which contains the liquid solvent 
metal and which consists of a vacuum-walled glass flask in a con- 
stant-temperature salt bath regulated to +0.02°C. A long neck 
admits thermocouple leads and a stirrer and allows introduction 
under vacuum of solute samples. These are held at 0°C until 
dropped into the metal bath. The calorimeter is calibrated with a 
sample of tungsten of known heat capacity. A temperature change 
of 0.01°C in the metal bath can be observed. 

Initial work was carried out with tin as solvent and additions 
of silver and gold-silver alloy in the annealed and cold-worked 
state. An alloy of 75 wt. percent gold and 25 wt. percent silver 
has been found to have a heat of solution in tin at 240°C which is 
just sufficient to heat the alloy from 0°C to 240°C; consequently, 
this alloy can be added with an over-all thermal effect near zero. 
Thus, the stored energy, found by taking the difference between 
the total heat effects resulting from adding worked and annealed 
samples, is almost the whole of the effect observed. Accidental 
variations may be largely eliminated by alternately adding an- 
nealed and worked samples in a single run. 

The alloy sample had been cold-worked by rolling wire of 0.032- 
in. diameter to ribbon of 0.0048 in.-thickness. With this alloy two 
runs were made with four additions each. The second and fourth 
samples of Run A and the first and third samples of Run B were 
cold-worked; the other samples were annealed. Since the heat of 
solution depends on the concentration, the energy of cold work is 
found by taking differences between worked and annealed samples 
at the same composition. Table I summarizes the results from 
which an average value of 0.28+0.1 calories per gram is obtained. 

Earlier measurements with silver cold-rolled from wire of 0.032- 
in. diameter to ribbon 0.0065 in. thick gave a value of 1.30.5 
calories per gram, but the device of alternating annealed and cold- 
worked samples in a single run was not used. 
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TABLE I, 








Addition No. 
2 3 





Composition of resulting 0. 1.75 
solution (wt. percent) 
Run A—heat effect (cal/g) 
Run B—heat effect (cal/g) 
Energy of cold work (cal/g) 


0.59* 
0.78 
0.19 


2 
0.64 
0.33* 
0.31 








* Cold-worked samples. 


Details of experimental procedure and heats of solution of group 
IB metals in tin will be published elsewhere. Further investigations 
of the stored energy of cold work and especially its dependence on 
such variables as the extent of working are planned. 

The writers are greatly indebted to Dr. John Chipman for 
suggesting the use of a metallic solvent and also acknowledge 
much helpful advice given by Dr. Carl Wagner. The work was 
performed under Contract No. AT-30-GEN-1002 with the AEC. 


1 For a review of the literature see W. G. Burgers in G. Masing, Handbuch 
der os (Akademische Verlagsgesellschaft, Leipzig, 1941), III, 2. 
p. 96 ff. 





Erratum: The Emission of Radiation from Diatomic 
Gases. I. Approximate Calculations 
[J. Appl. Phys. 21, 685-695 (1950) ] 


S. S. PENNER 


Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, California 


HE quantity am appearing in the second equation from the 
bottom of p. 692 should be replaced by a@n_i,n. Equation 
(19) should read 
an = rer N7/3yc= 2 ant, n- (19) 
Thus the integrated absorption for a given number of molecules 
Nr and for all vibrational transitions of the form n—1-—n is, 
approximately, independent of the degree of vibrational excitation 
and hence of the temperature. At moderate temperatures this 
statement is seen to be in agreement with the result of a more 
rigorous analysis by B. L. Crawford, Jr., and H. L. Dinsmore, 
[J. Chem. Phys. 18, 983, 1682 (1950) ]. 
Recent integrated.absorption measurements on CO [J. Chem. 
Phys. 19, 807, 817, 974 (1951) ] have shown that ao,:=237 cm 
atmos and ao2= 1.64 cm™ atmos™ for this molecule. 





Column Behavior under Conditions of 
Compressive Stress Wave Propagation 
GEORGE GERARD, College of Engineering, New York University, 
New York, New York 
AND 


HERBERT BECKER, Combustion Engineering-Superheater, Inc., 
New York, New York 


(Received July 26, 1951) 


T is the purpose of this note to show that under certain condi- 
tions of impact loading, the compressive stress supported by 
a perfect column at buckling may be of any magnitude in excess 
of its static euler load. This phenomenon is due to the fact that 
the magnitude of the stress at the end of the bar which is intro- 
duced by velocity impact depends only upon the velocity of impact 
and the density and elastic modulus of the material. Initially, the 
stress is only at the struck end. Since there is no stress in the 
remainder of the bar, it is only after the stress has been trans- 
mitted over a “critical” length of the bar that buckling can 
possibly occur. 





In the following analysis, it is assumed that the stresses are” 
elastic and, therefore, the stress wave is propagated as a single 
discontinuity and is constant over the length traversed by the 
wave. The equilibrium differential equation for the column is 


P'Werzz + Wrz = 0, (1) 


where: p is radius of gyration, ¢ is strain, w is lateral displacement, 
and x is coordinate along length. 

The assumed boundary conditions are clamping at the fixed end _ 
of the bar (x=0) from which point the compressive stress wave 
originates. The additional boundary conditions must be prescribed 
not at the other end of the bar but at the stress wave front where 
the first influence of the disturbance is felt. Since there can be no 
lateral deflection in front of the stress wave, continuity conditions 
prescribe the boundary conditions at the wave front. Thus, 


w(0, t) =wz(0, t) =0 
w(L., t) =w,(L,, t) =(), 


where L,=ct and is the distance the stress wave front has traveled 
along the bar at velocity c. A selution which satisfies Eq. (1) as 
well as the boundary conditions of Eq. (2) is given by 


w=a(1—cos2rx/L.) for O<x<I,. (3) 


By substituting the appropriate derivatives of Eq. (3) into Eq. (1), 
the following nontrivial solution is obtained 


Ler=2rpe. (4) 


Figure 1 is a photograph of a 0.010 in. by 0.50 in. aluminum 
alloy column which was buckled by the propagation of a straim 
wave of 0.0064 in./in. magnitude. It can be observed that a serie 
of buckles of approximately } in. wavelength formed above the 
end which was clamped in a jig. From Eq. (4), the estimate 
critical length based on the assumption of elastic buckling is 0 
in. which is in excellent agreement with the observed value. | 

It is concluded that both experiment and theory indicate thal 
under conditions of velocity impact in which the magnitude of th 
compressive stress wave may be far in excess of the static eulé 
load, buckling occurs in a critical length less than the actual lengt 
of the bar. This phenomenon indicates that a column can momeé 
tarily support a compressive stress of any magnitude which m 
be introduced by velocity impact. Therefore buckling un¢ 
velocity impact does not place an upper limit to the compressi 
stress which can be carried by a slender bar in the sense of § 
buckling. 


(2) 





